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Создан в целях информирования читательской аудитории о новейших достижениях и перспективах в 
области механики, машиностроения, информатики и вычислительной техники. Издание является форумом 
для сотрудничества российских и иностранных ученых, способствует сближению российского и мирового 
научно-информационного пространства. 


Журнал включен в перечень рецензируемых научных изданий, в котором должны быть 
опубликованы основные научные результаты диссертаций на соискание ученой степени кандидата 
наук, на соискание ученой степени доктора наук (Перечень ВАК) по следующим научным 
специальностям: 


1.1.7 — Теоретическая механика, динамика машин (технические науки) 

1.1.8 — Механика деформируемого твердого тела (технические, физико-математические науки) 

1.1.9 — Механика жидкости, газа и плазмы (технические науки) 

1.2.2 — Математическое моделирование, численные методы и комплексы программ (технические науки) 

2.3.1 - Системный анализ, управление и обработка информации, статистика (технические науки) 

2.3.3 - Автоматизация и управление технологическими процессами и производствами (технические науки) 

2.3.5 — Математическое и программное обеспечение вычислительных систем, комплексов и компьютерных сетей (технические 
науки) 

2.3.7 - Компьютерное моделирование и автоматизация проектирования (технические, физико-математические науки) 
2.3.8 — Информатика и информационные процессы (технические науки) 

2.5.2 — Машиноведение (технические науки) 

2.5.3 — Трение и износ в машинах (технические науки) 

2.5.5 — Технология и оборудование механической и физико-технической обработки (технические науки) 

2.5.6 — Технология машиностроения (технические науки) 

2.5.8 — Сварка, родственные процессы и технологии (технические науки) 

2.5.9 — Методы и приборы контроля и диагностики материалов, изделий, веществ и природной среды (технические 
науки) 

2.5.10 — Гидравлические машины, вакуумная, компрессорная техника, гидро- и пневмосистемы (технические науки) 


Индексация: РИНЦ, СуБе епшка, СтоззКеё, Оипепз1оп$, РОАТ, ЕВ$СО 

Наименование Свидетельство о регистрации средства массовой информации ЭЛ № ФС 77 - 78854 
органа, от 07 августа 2020 г., выдано Федеральной службой по надзору в сфере связи, 
зарегистрировавшего информационных технологий и массовых коммуникаций 
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Построение диаграммы предельных деформаций формоизменения 
листовых заготовок из авиационных алюминиевых сплавов 


С.И. Феоктистов, И.К. Андрианов 
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Аннотация 

Введение. Современное развитие штамповочного авиастроительного производства неразрывно связано с 
оценкой предельных возможностей листовых заготовок. Однако малоизученным является вопрос 
бездефектного формоизменения заготовок из авиационных алюминиевых сплавов. Важность данного вопроса 
связана с тем, что алюминиевые сплавы достаточно часто используются при изготовлении тонкостенных 
изделий авиационного назначения. При реализации процессов формообразования возможно появление 
различных дефектов — гофрообразования или недопустимого утонения. В связи с этим целью работы являлось 
построение диаграммы предельных деформаций основных авиационных сплавов и проведение сравнительного 
анализа кривых предельного деформирования для данных материалов. 

Материалы и методы. Для учета больших деформаций были использованы логарифмические деформации, 
обладающие свойством аддитивности. Построение диаграммы предельных деформаций формоизменения 
проводилось в постановке деформационной теории пластичности. Вопрос построения диаграммы предельных 
деформаций рассмотрен на основании критерия положительности производной силы нагружения. В области 
отрицательных значений наименьших главных деформаций для построения кривой предельного 
деформирования использовался критерий Хилла, а в зоне положительных значений главных наименьших 
логарифмических деформаций — критерий Свифта. При построении диаграммы предельного деформирования 
использовалась степенная аппроксимация закона упрочнения. 

Результаты исследования. Получены кривые предельных деформаций для авиационных сплавов: АМг-6, 
Д16АТ, АМГг2М, 1201-Т, АМцМ. Согласно проведенному сравнительному анализу областей безопасного 
формоизменения, сопоставлены значения деформаций начала шейкообразования и их влияние на изменение 
положения кривой предельного деформирования заготовок: чем больше деформация шейкообразования, тем 
выше положение кривой предельных деформаций. Описана концепция диаграммы предельных деформаций 
Килера. Представлены подходы к построению критериев Хилла и Свифта, используемых по результатам 
испытания листовых образцов на разрыв. 

Обсуждение и заключения. На основании построенных кривых предельных деформаций для авиационных 
сплавов АМг-6, Д16АТ, АМг2М, 1201-Т, АМцМ выяснили, что наибольшую область безопасного 
формоизменения имеет сплав АМг2М, наименьшую — сплав 1201-Т, что объясняется отличием относительных 
деформаций начала шейкообразования. Проведенное исследование позволило оценить возможности 
бездефектного формоизменения тонкостенных заготовок из основных авиационных алюминиевых сплавов. 
Применение построенных диаграмм предельного деформирования позволит прогнозировать появление 
разрывов в процессе формообразования листовых заготовок. 


Ключевые слова: листовая штамповка, диаграмма предельных деформаций, логарифмические деформации, 
диаграмма Хилла-Свифта. 
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Введение. Современное развитие авиационного производства неразрывно связано с исследованием 
процессов листовой штамповки. Одной из ОСНОВНЫХ проблем в задаче формообразования является 
прогнозирование дефектов, в частности, утонений, разрывов, гофрообразования. Данные вопросы связаны с 
оценкой предельных возможностей заготовки. Под предельным формообразованием листовой заготовки будем 
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понимать способность материала деформироваться до требуемой геометрии без шейкообразования или 
разрушения. 

На сегодняшний день проблема прогнозирования дефектов листовых заготовок в процессе штамповки 
решается с помощью следующих методов: 

— эмпирических, основанных на механических испытаниях на простое растяжение металлических образцов, 
тонких листов и лент, испытаний на изгиб, а также методов испытаний труб на раздачу и протяжку; 

— теоретико-эмпирических, которые основаны на использовании и распространении результатов испытаний 
образцов на одноосное растяжение на другие схемы деформирования заготовок; 

— теоретических, в основе которых лежит использование критериев предельного деформирования, в 
частности, при изготовлении тонкостенных изделий. Основоположниками этих методов были Дж. Закс, 
Р. Хилл, А. Д. Томленов, В. Д. Головлев, Г. Д. Дель, 3. Марчиняк, А. Д. Матвеев, Дж. Д. Любан [1-20]. 

Следует отметить, что недостатком эмпирических и теоретико-эмпирических методов является 
ограниченный круг использования результатов. 

Наиболее важным шагом в решении вопроса прогнозирования дефектов тонкостенных изделий стала 
разработка концепции диаграммы предельных деформаций (Е.О), предложенная Килером [16-20], которая на 
сегодняшний день является общепринятой в решении задач листовой штамповки. ЕГО-диаграммы широко 
используются в САЕ-программных системах АОТОЕОКМ и РАМ-5ТАМР 2С. 

Экспериментальные методы построения диаграмм деформирования основываются на методике испытаний, 
представленных в трудах Марчиняка и Наказимы. Следует также отметить, что вопросы пластического 
разрушения листовых заготовок рассматривались в [17-19]. В последнее десятилетие интерес к построению 
ЕГО-диаграмм существенно вырос, большинство зарубежных исследований направлено на экспериментальное 
построение диаграмм для конкретных материалов, а также численное моделирование процессов 
формоизменения с помощью конечно-элементных методов, что отражено в трудах [20-30]. Теоретические 
аспекты построения ЕГО-диаграмм и диаграмм деформирования Ш рода представлены в работах [31, 32]. 

Важно отметить, что диаграмма предельных деформаций позволяет оценить начало шейкообразования, 
которое в процессе деформирования заканчивается разрушением образца. ЕГО-диаграмма связывает значения 
главных логарифмических деформаций, действующих в плоскости листа. Диаграммы предельных деформаций 
позволяют не только прогнозировать разрушение заготовки, но и оценивать наличие других дефектов, в 
частности, складкообразование, утонения, которые, в свою очередь, понижают качество штампуемой детали. 
Основными зонами ЕГО-диаграммы являются зоны разрушений, возможных разрывов, безопасного 
формоизменения, вероятного образования складок и складкообразования (рис. 1) [31]. 


Рис. 1. Диаграмма предельных деформаций: 1 — зона разрушений; 2 — зона возможных разрывов; 
3 — зона безопасного формообразования; 4 — зона вероятного образования складок; 5 — зона складкообразования [31] 


Материалы и методы. При описании процессов формообразования тонкостенных заготовок момент 
перехода в стадию пластичности определяется в соответствии с критерием Губера-фон Мизеса [33]: 


9; = \/0? + 02 — 0102 = от, 
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где 01, 02) — главные истинные напряжения, при 03 = 0 ввиду малой толщин заготовки; 0; — интенсивность 


истинных напряжений; от — предел текучести. 
В соответствии с деформационной теорией пластичности зависимость между интенсивностью напряжений и 
интенсивностью логарифмических деформаций определяется как: 


ый 


ия (и-) | 
о (0 
е2 = (92 —тоз), 
где е,, е›, — главные деформации. 


Поскольку в задачах листовой штамповки процессы формоизменения могут происходить в несколько 
переходов, Поэтому рассматриваются большие деформации. Применение относительных деформаций 
недопустимо. В связи с этим деформированное состояние в (1) представлено в истинных логарифмических 


2 
е; = — [е? +е2 + е\е. 
МЕ] 


Отношение главных деформаций и главных истинных напряжений: 


Ще в — 92 
=, В =. 0) 


деформациях. 
Интенсивность главных деформаций: 


На основании соотношения (1) и (2), связь между а и В определяется: 


В — а. (3) 


Согласно выражению (2), критерий Губера-Мизеса имеет вид: 


в; = 0/1 - В+ В? = от, 


а интенсивность деформаций: 


де +а +а?. (4) 


Для оценки начала шейкообразования на сегодняшний день используется критерий положительности 
производной силы нагружения, основоположниками которого являются Дж. Закс и Дж. Д. Любан [17]. 
Согласно описанному критерию, деформирование образца является устойчивым при положительном 
приращении растягивающего усилия. Момент неустойчивого деформирования с последующим растяжением 


е; = 


начинается при ДР = 0 и продолжается при ДР < 0 (рис. 2). 


Рис. 2. Индикаторная диаграмма одноосного растяжения образца [33] 


Диаграмму деформирования алюминиевых сплавов, получаемую по результатам испытания на одноосное 
растяжение, при теоретическом анализе диаграмм предельных деформаций формоизменения аппроксимируют 
степенной функцией [18, 19, 31, 32]: 

05 = Ае" или д; = Ае!, 
где ох = Р/Е — истинное напряжение; Р — растягивающее усилие; Е— текущая площадь сечения образца; 
е = ш(1 + ЛИ/Ё) ,аАит — коэффициенты степенной аппроксимации. 
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Используя критерий положительности производной силы нагружения и степенную аппроксимацию 
диаграммы деформирования Ш рода, можно получить связь между предельной деформацией образца в момент 
возникновения рассеянной шейки и коэффициентом степенной аппроксимации п при условии 
ДР = 0 (рис. 3) [32]: 

@ш = Фи = П. (5) 


Важно отметить, что соотношение (1) выполняется при использовании логарифмических деформаций. В 
случае использования относительных деформаций ВОЗМОЖНО ТОЛЬКО приближенное равенство. 


Ртах 


01 


а) 6) в) 


Рис. 3. Возникновение рассеянной шейки при одноосном растяжении плоского образца: а — образец; 
б — напряженное состояние; в — деформированное состояние [33] 


Рассмотрим напряженное состояние пластины, к которой по краям приложены два растягивающих усилия, 
т.е. заготовка испытывает двухосное растяжение (рис. 4). 


02 


6) в) 


Рис. 4. Напряженно-деформированное состояние пластины при двухосном растяжении: 
а — пластина; б — напряженное состояние; в — деформированное состояние [33] 


Используя критерий положительности производной силы нагружения, определим величину предельной 
деформации е; и в момент, когда силы Р, или Р› максимальны. Для случая, когда Е›/Ё, = соп$%, Р./Р, = соп$, 
в момент максимума растягивающих усилий АР, = АР. = 0. Г. В. Свифтом [19] предложено соотношение для 
предельной деформации: 

(1 -В+ 823/72 

Согласно выражениям (2)—(4), соотношение для описания кривой предельных деформаций при различных 

показателях степенной аппроксимации в зоне упрочнения имеет вид: 


4(е, —п)(е› + 2е,)3 — З(е, — 2п)(е, + 2е)(е› + 2е.)? — 
—3(е, + 2п)(е› + 2е,)(е. + 2е>)? + 2(2е. +п)(е, + 2е>)3 = 0. (6) 


еш = 4" 


Из экспериментальных и теоретических работ известно, что после возникновения рассеянной шейки 
пластическая деформация образца продолжается. В дальнейшем может возникнуть локализованная шейка, 
которая отличается от рассеянной не только размерами, но и тем, что её возникновение и развитие 
осуществляются в условиях плоской деформации при интенсивном утонении образца по толщине. 

Согласно исследованиям Р. Хилла [20], критерий предельного деформирования определяется моментом 
образования локальной шейки, при котором приращение суммарного усилия равно нулю. В этом случае 
соотношение для предельной деформации определяется выражением: 

ен = 2т 622422 
т 1+В 


(т) 
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С учетом выражений (2), (4), выражение (7) для построения диаграммы предельных деформации по 

критерию Хилла имеет вид: 
е. +е› — п = 0. (8) 

Следует также отметить, что описанный подход согласуется с методом конечных элементов, который 
широко применяется для построения ЕГ.О-диаграмм различных материалов в последние годы [21-29]. 

Результаты исследования. На практике построение ЕГ.О-диаграммы будем осуществлять с использованием 
двух критериев [31], а именно: критерия Хилла, который используется при е› < 0 согласно (8); критерия 
Свифта, который используется при е> > 0 согласно (6). Рассмотрим применение данных соотношений для 
построения кривой предельных деформаций широко используемых в авиационной промышленности 
алюминиевых сплавов при известных значениях деформации шейкообразования: #„ = 0,18 (АМг-6), =. = 0,16 
(Д16АТ), =„ = 0,2 (АМг2М), =; = 0,06 (1201-Т), =. = 0,1 (АМЦМ) [34]. 

Определяя е„ по формуле е„ = ш(1+=):) и используя выражение (5), получим значение показателей 
деформационного упрочнения: п = 0,17 (АМг-6), п = 0,15 (Д16АТ), п = 0,18 (АМг2М), п = 0,06 (1201-Т), 
п = 0,09 (АМцМ). Тогда диаграмма Хилла-Свифта, согласно (6), (8), примет вид, изображенный на рис. 5. 
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Рис. 5. Н.О-диаграммы Хилла-Свифта для различных авиационных сплавов: 
1 — АМг-6; 2 — Д16АТ; 3 — АМг2М; 4 — 1201-Т; 5 — АМцМ 


Обсуждение и заключения. Согласно построенным кривым предельных деформаций наибольшая область 
безопасного формоизменения из пяти исследуемых сплавов у сплава АМг2М, наименьшая — у сплава 1201-Т, 
что обусловлено различиями в деформационных свойствах материалов, в частности, отличием деформаций 
начала шейкообразования. У сплава АМг2М относительная деформация начала шейкообразования 
составляет 20 %, ау сплава 1201-Т — 6 %. 

Таким образом, на основании данных о кривых упрочнения, степенной аппроксимации и деформации 
шейкообразования для авиационных алюминиевых сплавов АМг-6, Д16АТ, АМг2М, 1201-Т, АМЦМ построены 
кривые предельных деформаций формоизменения, позволяющие определить зону безопасного деформирования 
листовых заготовок. Результаты исследования имеют практическую значимость при решении задач листовой 
штамповки для данных материалов с целью прогнозирования недопустимых утонений, разрывов и 
складкообразований тонкостенных заготовок. 
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К вопросу о построении математических моделей 
мембранной теории выпуклых оболочек 
Е.В. Тюриков © 
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Аннотация 

Введение. В работе рассмотрены вопросы построения математических моделей безмоментного состояния 
напряженного равновесия упругих выпуклых оболочек с использованием методов комплексного анализа. При 
этом впервые рассмотрены оболочки с кусочно-гладкой (ребристой) боковой поверхностью. Целью работы 
являлось отыскание классов оболочек, для которых возможно построение содержательных математических 
моделей. 

Материалы и методы. С помощью методов теории разрывной задачи Римана-Гильберта для обобщённых 
аналитических функций получен критерий безусловной разрешимости соответствующей статической задачи 
для уравнения равновесия выпуклой оболочки с ребристой боковой поверхностью. Этот критерий в сочетании с 
методами теории обобщённых аналитических функций представляет собой инструмент построения 
математических моделей состояния безмоментного напряжённого равновесия упругих выпуклых оболочек. 
Результаты исследования. Разработан метод построения математических моделей безмоментного состояния 
напряжённого равновесия выпуклой оболочки при действии переменной внешней нагрузки и условии 
концентрации напряжений в угловых точках срединной поверхности. Введение в граничное условие векторного 
параметра, а также понятий «порядок квазикорректности» и «квазиустойчивость» позволяют провести как 
количественное, так и качественное сравнение математических моделей. Найдены классы оболочек, для 
которых описание математических моделей дается в терминах геометрии границы в окрестности угловых точек 
срединной поверхности. Полученный результат в применении к пологим выпуклым оболочкам позволяет дать 
геометрический критерий квазиустойчивости. Установлено, что для пологой оболочки, не являющейся 
квазиустойчивой, единственной адекватной математической моделью является вероятностная. 

Обсуждение и заключения. Предлагаемый метод построения двухпараметрического семейства задач с 
модифицированным граничным условием позволяет моделировать состояние безмоментного напряженного 
равновесия для достаточно широких классов выпуклых оболочек с кусочно-гладкой боковой поверхностью при 
условии втулочной связи. При этом разработанный алгоритм вычисления индекса граничного условия 
позволяет ответить на вопрос о существовании адекватной математической модели для оболочки с боковой 
поверхностью произвольной конфигурации, а для оболочек специального вида (например, пологих или 
оболочек вращения) — сформулировать геометрический критерий существования математической модели. 


Ключевые слова: тонкая упругая оболочка, обобщенная аналитическая функция, задача Римана-Гильберта, 
индекс граничного условия, математическая модель. 
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Введение. Рассмотренные в работе вопросы изучались и описывались еще в работах И. Н. Векуа [1,2] и 
А. Л. Гольденвезера [3, 4], положивших начало применению методов теории обобщенных аналитических 
функций к безмоментной (мембранной) теории тонких упругих оболочек и теории изгибаний поверхностей. К 
настоящему времени окончательные результаты в этом направлении получены для выпуклых оболочек с 
гладким краем (т.е. с гладкой границей её серединной поверхности). Наиболее значимые из них — 
корректность и квазикорректность основной задачи со статическим граничным условием с односвязной и 
многосвязной серединными поверхностями — есть следствия того факта, что индекс соответствующей задачи 
Римана-Гильберта является инвариантом связности поверхности. Применение автором методов И. Н. Векуа к 


Е.В. Тюриков. К вопросу о построении математических моделей мембранной теории выпуклых оболочек 


задачам теории выпуклых оболочек с кусочно-гладким краем’ [5,6] выявило связь между «геометрией» 
срединной поверхности в окрестности её угловой точки и картиной разрешимости соответствующих задач 
Римана-Гильберта с разрывным коэффициентом граничного условия. Использование этих методов в 
работах [7-9] позволило получить эффективную формулу для индекса при некоторых дополнительных 
геометрических условиях на угловые точки поверхности, и, как следствие, геометрический критерий 
квазикорректности основной граничной задачи. 

Цель настоящей работы — построение математических моделей безмоментного состояния напряжённого 
равновесия выпуклых оболочек с ребристыми боковыми поверхностями на основе геометрического критерия 
квазикорректности основной граничной задачи. 

Материалы и методы. Пусть 5 — односвязная поверхность заданного класса регулярности [7] с кусочно- 


п . 
гладким краем Г.=| | Г. и угловыми точками р. (1=1,...п). Зададим на 5 вдоль Г, кусочно-непрерывное 
ТБ р; 
т 


векторное поле г= {(5),В(5)} ‚ допускающее разрывы первого рода в точках р,, где а (5) ; В(5) (о +В? =1, 
В>0) — касательная и нормальная составляющие, 5 —щ натуральный параметр, и обозначим через Л 
заданный в[9] гомеоморфизм поверхности 5, на комплексную плоскость &=х-й{у. Пусть область 
р=/(5) — образ поверхности при отображении на плоскость < с границей Г и угловыми точками 
4; =У (р,). Рассмотрим следующую задачу (задача К): найти в области р комплекснозначное решение и (=) 


уравнения: 
%-(=)-В(2)*(5)=Е(<з), ЕР, (1) 
удовлетворяющее условию Римана-Гильберта 


Ве{^(6),*(5)} =1(5), СеГ, 2) 


где 


м6) = (6) [В© (6) -а(6) 56] (3) 
5(6)=5(6)+8. (6), (б)=а(б)+и, (6), Г=-—, 5, Е (1=1,2) — вещественнозначные функции, 
комплекснозначные функции 7(6) | ^(5) гёльдеровы на каждой из дуг Г, = Иь) ; 5 (>, +, ) ‚В(=), 


Е(=) — функции класса Г, (2), г>2. При этом предполагается, что решения класса М!” в точках разрыва 


4, допускают «интегрируемую бесконечность» т.е. допускают оценку |№( <] < с0п$ |= — 4; 


| 


‚ |: 


Следуя [10], класс таких решений обозначим через Н’. 
Как известно [11], статическая граничная задача безмоментной теории для упругой выпуклой оболочки с 


ребристой боковой поверхностью в математической постановке есть задача К, где и (=) —щ комплексная 


2 & * 
функция напряжений, ЁР(=) — комплекснозначная функция внешней нагрузки. При этом условие ЕН 


равносильно условию концентрации напряжений в угловых точках серединной поверхности. 

Построение математической модели состояния равновесия оболочки проведём на основе результатов о 
разрешимости задачи К для поверхности специального вида (канонического купола [9]). Для упрощения 
изложения будем считать, что в каждой угловой точке направление одной из дуг совпадает с одним из главных 
направлений К, (К, ) и называть купол 2-каноническим (1-каноническим). Задачу К для канонического купола 
К назовём канонической, если направление поля г в каждой угловой точке р есть направление обобщённой 

з “ . &2 
касательной [7, с. 46]. Введём обозначения: д; — отношение соответствующих главных кривизн К, К, в точке 
р (0<5, <1), р(\,) — угловая точка р с внутренним углом у,, Т(у,} — множество (сектор) направлений 


обобщённой касательной в этой точке, Т — множество непрерывных на Г, векторных полей г, задающих в 


каждой угловой точке Р(у й направление обобщённой касательной. Как установлено в [7], каноническая 


задача К квазикорректна для любого поля гЕТ ‚, если п>2. Задача К есть семейство К" задач (1)-(3), каждая 


' Тюриков Е.В. Обобщенная граничная задача Гольденвейзера для безмоментных сферических куполов. В: Современные проблемы 
механики сплошной среды: тр. ХУ межд. конф. Ростов-на-Дону; 2010. С. 290-293. 
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из которых задаётся выбором векторного поля г. Следуя И. Н. Векуа [1], задачу К" назовём 5 -квазикорректной 
в классе Н’, если она безусловно разрешима в этом классе, а её решение зависит от 5$ вещественных 
произвольных постоянных (5 — порядок квазикорректности). 

Определение 1. Каноническая задача ЕЁ называется квазиустойчивой относительно поля направлений 
обобщённой касательной, если задача В" 5-квазикорректна для любого поля ТЕТ. 

Замечание [. В силу теоремы о разрешимости задачи Римана-Гильберта для обобщённых аналитических 
функций [11] задача К квазиустойчива тогда и только тогда, когда индекс к задачи К есть инвариант поля 
ГЕТ. При этом индексом задачи назовём индекс соответствующей задачи сопряжения с коэффиицентом 


^(6) = (5) (5) . В случае к> -1, порядок квазикорректности 5 =к+1. 
Ниже используется техника [6, 8] вычисления индекса граничного условия вида (2), а также понятие [10] 


особенного узла р’ задачи (1), (2) или особенной точки 4 =У ( р,) разрыва граничного условия (2). Следуя [6], 


направление обобщённой касательной в угловой точке назовём особенным, если соотвествующая точка разрыва 
граничного условия (2), (3) есть особенный узел задачи К. 
Определение 2. Угловую точку р(\) назовём точкой неустойчивости задачи В, если сектор Т (у) 


содержит особенное направление. 
Введём обозначение: У, 6 — односторонние пределы в угловой точке р(у) касательного к Г, единичного 


вектора, причём вектор сб задаёт главное направление К, на поверхности в точке р, а внутренний угол У 


задаётся парой (-у,6). 


Утверждение 1. Если направление вектора г в точке р(\) совпадает с направлением вектора У, то 


точка а=У ( р) есть особенный узел граничного условия (2) тогда и только тогда, когда: 
1 
у = агссо$ ——. (4) 
1+5 


Если же вектор У заменить на вектор б, то: 
у= агссе Е я (5) 


где + — единственный положительный корень уравнения: 


1+64 1+6% Е 1 
2 + — (6) 

9+1 9+ к(1+г) +48 1 

Вы 1 
Здесь Е = 5 ‚ К=КК, [12, с.164]. При этом агссоз—_— <агсска у. 
Обозначим 
1 

6 = агссо$ ——, ц = агссме5\/Е . 7 
5’ га (7) 


Следствием утверждения 1 является утверждение 2. 


Утверждение 2. Угловая точка р (у) есть точка неустойчивости задачи К тогда и только тогда, когда: 


9<у<и. (8) 
Угловая точка р (у) есть точка 1-типа (2-типа), если выполнено условие 0<у<0 и<у< т) . Как 


установлено в [8], точка А-тина (К =1,2) есть точка устойчивости. Если р(У) — точка неустойчивости, то 
единственное особенное направление г, обобщённой касательной разбивает сектор Т(у) на два связных 
множества 7" (у) и 7 (у). 


* 
Утверждение 3. Индекс к задачи К в классе Н вычисляется по формуле: 


к=-4+У (4-к)), (9) 


11 
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где п — число угловых точек границы, к, =т для точки устойчивости р(\) т-типа (т = 12), ик, =5 в 


случае гЕТ® (У) ($=12) для точки р(\) неустойчивости. 


Замечание 2. Формула (9) для индекса задачи К в классе Н, ограниченных решений согласно [10] 


принимает вид: 
к=-4+\ (3-к). (10) 
1 


Отметим, что условие принадлежности решения задачи Ю классу Н’ есть условие ограниченности 
интеграла энергии растяжения оболочки [13, с. 83] в окрестности угловой точки. 

Из формул (9), (10) следует утверждение 4. 

Утверждение 4. Если граница Г. содержит единственную угловую точку [-типа, то задача К безусловно 
разрешима в классе Н’ и имеет единственное решение УтЕТ ; если п>2, п=п +п,, где п, — число угловых 
точек К- типа (К=1,2), то задача Е квазиустойчива относительно гЕТ с порядком квазикорректности 
5=2и, +п, -3. 

Если граница Г, содержит точки неустойчивости, то согласно (9) задача К в классе Н’ не является 
квазиустойчивым относительно поля геЕТ. Однако можно выделить такие классы полей, относительно 
которых задача К является квазиустойчивой с различными порядками квазикоретности. Классы таких полей 


можно задавать, выбирая в каждой точке неустойчивости р(у) направление обобщённой касательной только в 
1 2 
одном из секторов Т® (), т } (у). 
Замечание 3. Как нетрудно убедиться, формулы (9), (10) вместе с утверждением 4 справедливы и в случае: 
п 
и<У< +0”, (11) 
2 
где © — достаточно малая величина, заданная поверхностью 5. При этом условие малости @ не является 
обязательным. Например, для омбилической (к, =А,) точки р(У) достаточно положить ®’` = г 


Рассмотрим теперь 1-канонический купол К. В этом случае описание особенных узлов граничного 
условия (2) также даётся утверждением 2 с той лишь разницей, что в равенстве (4) и уравнении (6) величину 
2 К, 2 К, > > 
5`=-- <1 необходимо заменить на 5” = >1. При этом для угловой точки р(У) неустойчивости задачи К 
1 2 


неравенство (8) принимает вид: 


и акк =. (12) 


На основании очевидного графического анализа уравнения (6) заключаем: 
1° величина и= агссёе Е есть функция двух параметров, а именно: иц= и (5, К) ‚ КЕ (0; +00) ‚ если 


К | 
5° ыы ш= (5, К, е (0;-00) ‚ если 5* = 


1 2 


= 


УК Е (0;-+0о) ‚если 6’ = — ; 
1 


5—1-0 


2° Шт (6%) =5- Нтр(5,*,)= 


[а 


3° Нт н(6.6.)=5. Шт (5, =0 УК, с (0;-0) , если 5* 1; 


6—1+0 В 
4° для любого фиксированного дб из правой (левой) полуокрестности единицы функция И как функция 
д 

аргумента К, (К) есть медленно меняющаяся функция [14] в том смысле, что и(5,, к г Е (=) для любого 
фиксированного 5 = 5%, К, ® (0:40) (1=1,2). 

Результаты исследования. Введём в рассмотрение семейство поверхностей 5, те [0,=) ‚ где т — малый 
параметр, каждый их которых есть срединная поверхность тонкой упругой оболочки У. из некоторого 
семейства {У.}, причём У, и 5, совпадают с оболочкой У и поверхностью 5 соответственно. Будем 


полагать, что для Уте [0,=) выполнены следующие условия: 
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1) поверхность 5. есть канонический купол класса регулярности №”, г>2, с внутренними углами 
величины у, в угловых точках р, (1=1,...п) границы Г,, причём поверхность 5, есть 2-канонический (1- 
канонический) купол 5 с границей Г; 

2) отношения главных кривизн К“, К“ поверхности 5, в угловых точках совпадают с величинами 6” в 
угловых точках поверхности 5 ; 

3) ЕС ЕА, +5, (<), где те, (т)=0, К, = И 

Пусть / — отображение поверхности 5. на комплексную плоскость 5=х+й), заданное выше, 
р. =/1(5.) — семейство ограниченных в плоскости < односвязных областей с соответствующими Г, = / (1. } 


и угловыми точками а, = Л ( Р; ) Рассмотрим семейство задач К., те [0, =) ь 


%-(=)-В. (2)*(&)=Е (2), <ЕБ,, (13) 


Ве(^. (0) »(5)}=5.(6), СеГ. (14) 


в котором функции В. (=), ^.(6), 7.(6) определены серединной поверхностью 5, оболочки У. согласно [8], 
причём ^.. (6) имеет вид (3). 

Предполагается, что боковые поверхности оболочек У. имеют общие рёбра, содержащие точки р, 
(1=1,....п)). Отметим, что задачу (13), (14) Ух = [0,=) можно рассматривать как семейство задач во ‚ каждая из 
которых задаётся выбором векторного поля на 5. вдоль [.. 


Рассмотрим 1-канонический купол 5. Справедливо утверждение 4. 
Утверждение 4. Если в каждой из угловых точек р(у) купола 5 выполняется условие: 


о (15) 


а * [а .. а х 
то Уте [0,2=) задача К. квазиустоичива Н относительно поля направлении обобщённой касательнои с 


порядком квазикорректности 5 =Зп-3. Если же в каждой точке р(\) выполнено условие: 
Л 2 п 2 
агссе ру В (16) 


2 - 
где @’ определено замечанием 3, величина в, =®,(=) задана семейством 5,, те[0,=), то задача К, 


квазиустойчива в классе Н` с порядком квазикорректности 5=2п-3. 
Для доказательства достаточно воспользоваться утверждением 3, условиями 1-3 и известными 


свойствами [1, с. 97] отображения Л. Утверждение 4 остаётся в силе для 2-канонического купола 5, если 


условия (15), (16) заменить условиями: 


у < агссю УР — в, (17) 
о еЫИ. (18) 
1+5 2 


соответственно. 
Уточним теперь понятие угловой точки К -типа (к = 1.2) поверхности 5. Будем говорить, что точка р(\) 


границы Г есть точка 1-типа (2-типа) относительно семейства 5., те [0,=) ‚ если выполняются условия (15) 


и (17) для 1-канонической и 2-канонической точек соответственно (условия (16) и (18) для 1-канонической и 
2-канонической точек соответственно). Рассмотрим канонический купол 5, каждая угловая точка которого 
есть точка 1-типа или 2-типа относительно указанного семейства 5. Тогда следствием утверждения 4 является 


Утверждение 5. Если р и 4 — число точек [-типа и 2-типа соответственно (р +а= п), то УтеЕ [0,=) 
задача К. безусловно разрешима в классе Н” и квазиустойчива относительно направлений г обобщённой 
касательной с порядком квазикорректности 5 =3р+24-3. В частности, если р=0, то 4>2. 

Утверждение остаётся в силе, если класс Н” заменить классом Н, ограниченных в О. решений, а порядок 
5=3р+24-3 порядком 5=2р+а-3 при условии 2р+4>3. Таким образом, задача К. не является 


безусловно разрешимой в случаях: р=0, 432 ир=1, 4=1. 
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Случай омбилической (5=1) угловой точки Р(\) ‚ в которой любое направление является главным, 


р д аа 
требует специального рассмотрения. В этом случае согласно ВНЕ. а любое направление обобщённой 
касательной есть особенное направление задачи В. Однако и в этом случае утверждение 4 остаётся 
ь . т 
справедливым, если для семеиства {5} считать = 0 ‚ авеличины У, Ц заменить величинои З . 


Обсуждение и заключения. Полученные результаты могут быть использованы для построения 
математических моделей тонких и пологих оболочек положительной гауссовой кривизны с ребристыми 
боковыми поверхностями. Наиболее полные и продвинутые результаты как линейной, так и нелинейной теории 
упругих оболочек получены для тонких пологих оболочек. Детальное обсуждение понятия «пологая оболочка», 
а также описание различных вариантов теории приводится в [15, с. 29]. Линейная теория пологих выпуклых 
разработана И. Н. Векуа [2, 16]. В рамках этой теории вопрос о реализации состояния напряжённого равновесия 
пологой оболочки с ребристой боковой поверхностью при выполнении статического граничного условия 
общего вида сводится к задаче К , рассмотренной выше. 

Пусть Р — пологая оболочка [2, с. 164] с ребристой боковой поверхностью, 5 — её срединная 
поверхность с кусочно-гладким краем. Будем полагать, что в каждой угловой точке поверхности 5 


выполняется условие усиленной пологости К, = К, , которое равносильно следующему: 
0=1, (19) 


© К, К, 
гле б — любое из соотношений —, 


а. Условие (19) означает, что любую угловую точку Р(\) следует 
2 1 


считать как 1-каноничекой, так и 2-канонической точкой. Рассмотрим угловую точку р(\) ‚ в которой У=—,и 


о|а 


семейство поверхностей 5., те [0,=), заданное условиями (1)—(3). Согласно свойствам 1°—4° величин 0, ци 


утверждению 2 соответствующая задача К, Уте (0,=) не является квазиустойчивой в классе Н”, то есть точка 


п м“ К 
р(\) при УЯ 3 —ЩЩ точка неустоичивости задачи К. . В этом случае естественно предполагать, что в правои 


части формулы (10) соответств ющая точке У) величина к, (1<1<и) есть дискретная случайная величина 
Р 1 


с возможными значениями 1 и 2. Таким образом, если в каждой угловой точке р(у) серединной поверхности 


п 
выполнено условие У то в силу формулы (9) порядок квазикорректности задачи К есть дискретная 


случайная величина, принимающая значения Ш, т+1,..., т+п, где И — число угловых точек, т=2и-3 для 
класса решений Н’, и т=и-3 для класса Ну. 

Предлагаемый выше метод может быть использован для построения математических моделей теории тонких 
пологих оболочек с ребристыми боковыми поверхностями любой конфигурации. Для этого достаточно 
воспользоваться результатами” о разрешимости задачи К для сферических куполов с кусочно-гладким краем. 
Рассмотрим для определённости срединную поверхность 5 в предположении, что все угловые точки р(у) 


границы — «выходящие», то есть у<л. В этом случае угловая точка на сферической поверхности есть 
е И 
особенный узел граничного условия тогда и только тогда, когда -5 (^ а Отсюда следует, что 


«выходящая» угловая точка срединной поверхности пологой оболочки есть точка неустойчивости, если 
. дк 
выполнено одно из следующих условий: = З (^ =1.2 Таким образом, формула (9) для индекса может 


служить обоснованием следующей гипотезы: 

если задача К для пологой выпуклой оболочки безусловно разрешима в заданном классе решений, то её 
порядок квазикорректности есть дискретная случайная величина, принимающая целые значения К, К+1,..., 
К-М, где М — число точек неустойчивости; К — число, заданное набором угловых точек и выбором 
непрерывного векторного параметра г. 


? Тюриков Е. В. Обобщенная граничная задача Гольденвейзера для безмоментных сферических куполов. 


3 
То же. 
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В заключении отметим, что обоснованием этой гипотезы могут служить те же рассуждения, но проведённые 
для регулярных выпуклых поверхностей, удовлетворяющих условию локальной симметрии [17] в угловых 
точках. 
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Аннотация 

Введение. Статья посвящена проблеме подтверждения заданных уровней надежности при экспериментальной 
отработке сложной технической системы с последовательным соединением элементов. Такие задачи 
возникают, когда требуется принять решение об испытании системы в составе более крупной или об окончании 
экспериментальной отработки и запуске серийного производства. Цель исследования — обосновать 
сокращение сроков экспериментальной отработки. Задача — определить, принимается или отклоняется 
гипотеза Н.. 

Материалы и методы. Для реализации цели и задачи работы по результатам испытаний строится критическая 
область, описываемая неравенством. Формулировка задачи подтверждения требований базируется на 
известных подходах к проверке статистических гипотез. Задействуется понятийный аппарат теории 
информации, вероятности и статистики. Изучена теоретическая и прикладная литература о математических 
методах в теории надежности. Частные задачи работы решены известными способами. Так, вероятность 
получения точного числа успешных исходов в определенном количестве экспериментов определена по схеме 
Бернулли. Точный доверительный интервал, основанный на биномиальном распределении, получен из 
соотношения Клоппера — Пирсона. Теорема А. Д. Соловьева и Р. А. Мирного позволила оценить надежность 
системы по результатам испытаний ее компонент. 

Результаты исследования. Математически определены правила контроля, адекватные этапу 
экспериментальной отработки (при недостаточности данных о технической системе) и этапу серийного 
производства. Вероятность успешного исхода при испытании технических систем представлена через: 

— вероятность события для элемента системы; 

— значение доверительной вероятности; 

— требуемый объем испытаний. 

С этих позиций исследованы нулевая и альтернативная гипотезы и соответствующие им процедуры контроля 
надежности. Рассмотрены два положения. Первое допускает использование нулевой гипотезы доверия 
Н. => В. с альтернативой Н Ре Р: для подтверждения требований (Р,, у) к показателю 


надежности одного параметра при любых (Ру, у) 2 При этом достаточно одного безотказного испытания. 


Второе положение рассматривает последовательную техническую систему с № независимыми элементами, 
которые испытываются отдельно от системы по схеме Бернулли для одного параметра. Рассмотрим требования 
к системе в виде совокупности величин (Р., ) и требования к любому ее элементу (Р;,, у). Они совпадают, 


если планируемый исход испытаний соответствует случаям выполнения соотношения Р= Пт :Р, =Р 


нулевая альтернативная гипотеза выбирается из теории проверки статистических гипотез. 
Обсуждение и заключения. Стратегию экспериментальной отработки следует реализовать в два этап: поиск и 
подтверждение надежности элементов серией безотказных испытаний. В этом случае планируемый объем 
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испытаний каждого элемента определяется с учетом доверительной вероятности, нижней границы 
доверительного интервала и требований к показателям надежности одного параметра технической системы. 
Если допустимо использование нулевой гипотезы доверия, для подтверждения требований к показателю 
надежности достаточно одного безотказного испытания. 


Ключевые слова: экспериментальная отработка, проверка статистических гипотез, надежность технической 
системы, нулевая гипотеза, альтернативная гипотеза, гипотеза недоверия, гипотеза доверия, доверительная 
вероятность, объем безотказных испытаний, модель испытаний биномиального типа, схема Бернулли, 
уравнение Клоппера — Пирсона, теорема А. Д. Соловьева и Р. А. Мирного. 
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Введение. Рациональные методы подтверждения заданных уровней надежности актуальны для 
экспериментальной отработки сложной технической системы, когда принимается решение о возможности ее 
испытания в составе более крупной структуры или о завершении экспериментальной отработки и начале 
серийного производства. Такие же задачи возникают и при серийном производстве, если нужно: 

— оценить готовность предприятия к выпуску серийной продукции на основании данных испытаний 
установочной партии; 

— сделать вывод о соответствии выпускаемой продукции требованиям технической документации с учетом 
данных по эксплуатации. 

Цель работы — получить приемлемое решение для планирования и сокращения объема испытаний с 
помошью методов интервального оценивания показателей надежности последовательных технических систем. 
Для реализации заявленной цели необходимо определить, принимается или отклоняется гипотеза Н.. 

Материалы и методы. Задачу подтверждения требований целесообразно формулировать в терминах 
теории проверки статистических гипотез [1-4]. Пусть Р — надежность технической системы, Рг — некоторый 
фиксируемый (требуемый) уровень для Р. До проведения испытаний относительно Р можно выдвинуть три 
исходных предположения: Р = Р,, Р < Р,,Р > Р.. 

Каждое из них называется нулевой гипотезой, если записано виде: 

Пе РР = ре РВ 

Множество Н, = {Р =Р,} содержит лишь один элемент, поэтому гипотеза М, = {Р = Р,} называется 

простой. Гипотезы вида Н,= {Р <Р,} и Н, = {Р>Р,} называются сложными. Наряду с нулевой 


гипотезой, выражающей заранее сформулированную точку зрения, задают альтернативную гипотезу Н, 

выражающую противоположное высказывание Н,(Н,^ Н =О). Задействуем понятийный аппарат теории 

информации [1], вероятности и статистики, применимый к решению подобных задач. Рассмотрим две 
совокупности множеств Нь и Н: 

РЕ РР (р 

Н=РОРЬН = РВ (2) 

Гипотезу Н. в (1) будем называть жесткой, или гипотезой недоверия. Действительно, в случае (1) 


первоначально (до испытания) исходят из позиции недоверия к уровню качества системы. Показатель 
надежности Р принимается не выше некоторого фиксируемого уровня Р”,. Гипотезу Н, в (2) будем называть 


гипотезой доверия, так как в этом случае первоначально исходят из положения, что показатель надежности Р не 
меньше некоторого фиксируемого значения РУ. 

Смысл величин Р’, и РГ различен. В (1) Р’, — такое браковочное значение, что при Р< Р’, система 
считается неприемлемой. В (2) Рг — такое значение, что при Р > Р, система считается приемлемой для 
использования. Очевидно, что Р, > Р’,. 


Результаты исследования. Итак, необходимо определить, принимается или отклоняется гипотеза Но. В 
теории статистических гипотез для этого по результатам испытаний строится критическая область. Она 
описывается некоторым неравенством. Причем нулевой гипотезы (вследствие исходного Доверия к ней) 
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придерживаются до тех пор, пока это разумно с точки зрения принятого уровня значимости а. Поэтому заранее 
ясно, что процедура контроля надежности в случае (1) будет существенно отличаться по сравнению со 
случаем (2). 

Действительно, далее убедимся, что для отклонения гипотезы Н, в (1) и принятия гипотезы 
Н = {Р > В о выполнении требования к показателям надежности следует использовать критическую 


область (или условие): 


Р>Ру, (3) 
где Р — нижняя граница доверительного интервала для Р при значении доверительной вероятности 
у =1-а; Р’, — нижняя граница браковочного интервала для Р при значении доверительной вероятности 


у=Т-а. 
В случае (2) для принятия гипотезы Н, о соответствии значения параметра Р предъявляемому требованию 
следует использовать условие: 

РВ, (4) 
где Р — верхняя граница доверительного интервала для Р при значении доверительной вероятности 
у=1- о. 

В (3) и (4) под требованием к показателю надежности одного параметра Р понимается совокупность величин 
(Р’,, у) или (Р+т, у) , задаваемых до проведения испытаний. 


Пусть в условиях схемы Бернулли проведено одно успешное испытание. Тогда, используя соотношения 
Клоппера — Пирсона, находим нижнюю и верхнюю границы одного параметра при значении, например, 
у= 0,95: 


Р=(1- у) =1-т= 0,05; РЕ. 
При этом даже для весьма умеренных значений Р’, [0, 05; 0, 95] условие (3) не выполняется, в то время 


как (4) выполняется при любом Ру. Покажем справедливость принятого положения. 


Первое положение. Если допустимо использование нулевой гипотезы доверия Н, = |Р > В] при 
альтернативе Н = {Р < Вы то для подтверждения требований (Р,, у) к показателю надежности одного 
параметра при любых (Р., у) достаточно одного безотказного испытания. 

Если исходная гипотеза Н, — это гипотеза недоверия из (1), то необходимо существенно большее число 
испытаний. Так, при т = 0, получаем п> [ов (1 _ у)/ ю8Р; >>1. Это вполне справедливо, т. к. при проверке 


гипотез первоначально исходят из справедливости нулевой гипотезы Но. 

На этапе экспериментальной отработки нет достаточно полных данных о технической системе, поэтому 
целесообразно использовать правило контроля (3). На этапе серийного производства можно исходить из 
гипотезы доверия и использовать существенно более легкое правило контроля (4). Это допустимо, если, по 
данным экспериментальной отработки, условие (3) оказалось выполненным. 

Рассмотрим систему, состоящую из № независимых последовательно соединенных элементов, которые 
могут испытываться отдельно. Тогда вероятность успешного исхода при испытании технических систем: 


№ 
Р=ПР. (5) 
1 
Здесь Р; —Ш_ вероятность того же события ДлЯ 1-го элемента. К величине Р заданы требования в виде 


совокупности значений (Р,,у). Нужно спланировать процедуру контроля надежности одного параметра для 


каждого элемента системы, т. е. задать У! е [1]. п] пару (Р„., 7). 
№ 
Вследствие перемножения Р; в формуле (5) должно выполняться соотношение [П]Р,=Р. Кроме 
7 


того, /;, = У. 
В результате требуемый объем испытаний и; каждого элемента резко возрастает и даже при безотказных 
исходах всех испытаний становится неприемлемым. При РГ = 0,9 иу = 0,95, М = 100 ит=0, = М: 


Р-Р, = 0,999, п, =108(1 — 7) /1юеР,, = 3000. 


Такому способу планирования логически противоречит неравенство и, >п,, следующее из Р,, >Р; при 


т; = 0. Понятно, что при безотказных исходах требуемый объем испытаний п»; 1-го элемента, проводимых 


Механика 


29 


Бер://уези-Чопза.ги 


А4уапсеа Епетеетпх Кезеагсй (Ко5оу-оп-Ооп). 2023;23(1):26-33. е155М 2687—1653 


отдельно от системы, должен быть равен требуемому объему испытаний системы п.. Во избежание этого 
противоречия следует использовать теоремы А. Д. Соловьева и Р. А. Мирного [5-7]. Итак, при т = 0, У! = 1,М№: 
РЕШШ:Р.=/ (п, 0,5) = (1-1)”. (6) 
Здесь Р — нижняя граница доверительного интервала для показателя надежности технической системы по 
одному параметру, при значении ‘7 доверительной вероятности; Р; — значение нижней границы доверительного 
интервала для показателя надежности 1-го элемента системы при той же доверительной вероятности; и — 
минимальное число испытаний элементов системы; /(п,0,у) — корень уравнения Клоппера — Пирсона: 


т 


1 ==», м Вы РТ, (7) 
К=0 
Согласно [8—12]: 
Г(п,п4,у) < Р< (в, [па]. (8) 
— М = = — 
Здесь а =1-Р; Р =П (1 _ т, /п,); п = пт п‚; [14] — целая часть произведения па; }(и,па,у) — 


корень уравнения Л»(п,Р,п4 +1) == 


Из (8) следует: 
Р=шш:Р, =Рт, (9) 


1<1<М№ 
где Р „ — минимальное из Р; при значении доверительной вероятности у. 
Р, =шш:Р, > Р, <> (Р/> Р,)с(Р, > Р,)...С(Ру> Р;). 


т . 
1<1<М№ 


Это справедливо не только для случая (1), когда т=0, У! = 1, М, но и для случая (2), т.е. для исхода 
п= (п,п,, ...пу), т= (т, 0, 0, ..., 0) испытаний, где п и т — вектор испытаний и вектор отказов, если 
п =п. 

При этом отказывает только один элемент, подвергшийся испытаниям минимальное число раз. 


Действительно, в этом случае пд =т, — целое число. Расчеты позволили установить, что (9) приближенно 


выполняется Также при исходе испытаний пт, т.е. в случае (3), если для пары (п, ‚т, ) возможно 


Во всех упомянутых случаях (9) не перемножаются нижние границы, а система вырождается в один 
слабейший элемент. Это позволяет подтвердить следующее положение. 

Второе положение. Рассмотрим последовательную систему с № независимыми элементами, которые 
испытываются отдельно от системы по схеме Бернулли для одного параметра. Требования, задаваемые к 
системе в виде совокупности величин (Р,, 7) ‚ и требования к любому ее элементу (Р,., у,) совпадают, если 


планируемый исход испытаний соответствует упомянутым случаям выполнения соотношения (9), а нулевая 
альтернативная гипотеза выбирается, исходя из (1) и (2). 

Следствие. В случае (3) выполнение (9) — планируемый объем безотказных испытаний (№М- 1) элементов 
определяется из: 


п, > п, = [08 (1 -%) / 08Р;.. (10) 
Объем отказов условно первого элемента при ти находится из соотношения: 
Р, =Л(и,т,у) =Рг. (р 


Доказательство (11) основано на том, что условие Р =Р,> Р, выполняется, если Р, =Р, = Р, и 
Р, <М![2, №] . Последнее же соотношение удовлетворяется, если выполняется (10), т. к.: 
— т 
п, >2п >Р, = (1 -у)">РВ,. 
Во всех рассмотренных случаях предполагалось, что до проведения испытаний нет никакой информации 
относительно Р, если не считать очевидный факт Ре[0, 1]. Однако может быть известно, что Р > Р,, где 


Рн= 0. Значит, Р = Ре [ Ри, 1] . Величину Рн можно найти по данным испытаний или расчетов к моменту 


планирования испытаний на надежность. Еще нет методики определения Ри, и ее разработка —Щ задача будущих 
изысканий. Но если известна величина Ри, по формуле полной вероятности можно найти: 


Р =Р+ (Р. 


О.Ю. Царев и др. Подтверждение показателей надежности при экспериментальной отработке сложной технической 


Отсюда: 
Р =(Р-Р,)! (1 -Р,), Р= (Р-Р,)!А -Р,). (12) 


Последнее соотношение выполняется вследствие монотонности по Р, зависимости Р = (Р - Р,). 


С учетом (12) соотношение (2) примет вид: 


Р‚+(Г+Р, ) Г (п,па,у) <Р5Р, + (1 -Р,) (па, у) . (13) 


Пусть условием принятия решения о соответствии технической системы требованиям (Р,, у) по-прежнему 


служит (3), где нижняя граница доверительного интервала определяется из (13) с учетом Р= РЕ[Р., : 


Тогда из соотношения 


РВ, (19) 
находим планируемый объем безотказных испытаний по одному параметру для каждого из № элементов: 
п, > п, =108 (1 — 7) ю8(Р,‚-Р,) / а-Р,). (15) 
Величина и ’ убывает по Рин. Значит, пу (10) при Рн = бит» = 0 при Ри = Р.. 


№ 
Пример. Требования к показателям Р= ПВ надежности системы заданы для одного параметра в виде 
7 


совокупности величин (Рут = 0,90; у = 0,95). Число элементов технической системы М = 100. По имеющимся 
данным, Рг> Рн= 0,70. Требуется найти планируемый объем испытаний каждого из № элементов, если 
выполнение требований надежности проверяется условием (14). Из (15) находим: 

п 21% = [08 (0,05) / 108 (0,90 — 0,70) / (1-0,70) =7. 

Заметим, что при Рин = Оп; >п* = 29. 

Соотношения (10) и (15) позволяют планировать необходимый объем испытаний 1-го элемента технической 
системы при определенной последовательности экспериментальной отработки системы. Весь процесс 
экспериментальной отработки делится условно на два периода: поиск и подтверждение требований надежности 
для решения о переходе к следующему этапу испытаний или о принятии технической системы к серийному 
производству. В первом периоде возможны доработки, и целесообразно использовать модели с переменной 
вероятностью Р успешного исхода испытания системы. 

Данные, полученные в первом периоде, могут использоваться для расчета величины Ри. 

Во втором периоде имеем дело с установившимся вариантом конструкции технической системы и 
технологического процесса. Это позволяет использовать рассмотренные выше модели испытаний 
биномиального типа с постоянной вероятностью Р. 

Пусть первый период отработки № элементов технической системы закончен и ставится вопрос о 
подтверждении требований к показателю надежности системы по одному параметру. Целесообразно 
подтвердить надежность, если положительное решение принимается только в случае безотказного исхода 
испытаний последней серии для каждого из М элементов. Эта стратегия удобна тем, что она основана на 
минимально возможном числе испытаний элементов технической системы и допускает простые аналитические 
решения (10) и (15). 

В общем случае есть смысл исследовать стратегию, допускающую отказы элементов при проведении 
испытаний и основанную на оптимизации некоторой целевой функции. Но здесь ограничимся рассмотрением 
лишь упомянутой стратегии с безотказными заключительными сериями. 

Обсуждение и заключения. Результаты научных изысканий позволили сформулировать приведенные ниже 
ВЫВОДЫ. 

1. Даже при большом числе М элементов системы возможно планирование объема их испытаний. В этом 
случае методы интервального оценивания показателей надежности последовательных технических систем 
позволяют получить приемлемое решение, однако только по одному параметру. 

2. Стратегия экспериментальной отработки технических систем тесно связана с методом подтверждения 
надежности элементов. Рациональная стратегия экспериментальной отработки предусматривает после 
поискового периода подтверждение надежности элементов заключительной серией безотказных испытаний. В 
этом случае планируемый объем испытаний каждого из М элементов не зависит от М№ и определяется по 
соотношению (15), в которое входят требования (Рт,/) к показателям надежности одного параметра 
технической системы в целом и Рн. Получаемый по (15) объем испытаний каждого элемента технической 
системы при любом их числе невелик, если на систему из М№ элементов заданы умеренные требования 
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(Рг= 0,80 ...0,95; у = 0,90... 0,95), однако только по одному параметру. При этом объем убывает с ростом 
величины Ри. 

3. В серийном производстве при испытаниях модернизированных технических систем можно использовать 
метод контроля со сменой нулевой и альтернативной гипотез. Если допустимо использование нулевой гипотезы 
доверия, то для подтверждения требований к показателю надежности достаточно одного безотказного 
испытания. 
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Метод решения задачи о движении нагрузки по ледяному покрову — 
водоема по сложной траектории : 
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А.В. Галабурдин © 
Донской государственный технический университет, Российская Федерация, г. Ростов-на-Дону, пл. Гагарина, 1 


Ч сааБагат @ та|.га 


Аннотация 

Введение. Освоение полярных районов Мирового океана способствовало повышению интереса к изучению 
волновых процессов в водоемах с ледяным покровом, обусловленных действием подвижной нагрузки. В 
большинстве работ отечественных и зарубежных ученых ледяной покров рассматривался как упругая или 
вязкоупругая пластина, нагруженная прямолинейно движущейся вертикальной силой. Однако при 
моделировании воздействия транспортных средств на ледяной покров представляет интерес рассмотрение 
задач, в которых сила движется по более сложной траектории. Поэтому целью данного исследования является 
разработка метода исследования поведения ледяного покрова под действием силы, движущейся по траектории 
сложной формы, подчиняясь произвольному закону движения. 

Материалы и методы. Предложен метод решения задач о действии движущейся произвольным образом силы по 
ледяному покрову водоема конечной глубины, заполненного невязкой несжимаемой жидкостью. Ледяной покров 
рассматривался как вязкоупругая пластина, лежащая на поверхности жидкости, находящейся в состоянии 
потенциального течения. На верхнюю поверхность пластины приложена сосредоточенная сила, движущаяся по 
произвольной замкнутой траектории и являющаяся периодической по времени. На нижнюю поверхность пластины 
со стороны жидкости действует гидродинамическое давление. В силу периодичности приложенной к пластине 
нагрузки при решении дифференциальных уравнений, описывающих поведение данной системы, использовали 
интегральное преобразование по времени. Далее, применяя традиционные методы, получили формулы для 
вычисления напряжений и перемещений в пластине и компонент вектора скорости частиц жидкости. Эти формулы 
представили в виде повторного интеграла. Для вычисления интегралов использовали численные методы. 
Результаты исследования. Расчеты показали, что прогиб ледяного покрова заметно увеличивается с ростом 
скорости и касательного ускорения движения нагрузки. К росту прогиба также приводят увеличение времени 
релаксации льда и уменьшение радиуса траектории движения нагрузки. Рассчитано распределение вектора 
скорости частиц жидкости по глубине водоема. 

Обсуждение и заключения. Предложенный метод показал свою эффективность при решении задач о действии 
подвижной нагрузки на ледяной покров водоема. С его помощью исследовано влияние на прогиб ледяного 
покрова кривизны траектории движения и механических свойств льда, кинематических характеристик движения 
нагрузки. Метод моделирует воздействие транспортных средств на ледяной покров водоема. Результаты его 
применения могут быть использованы при строительстве ледовых дорог или аэродромов на льду. 


Ключевые слова: бесконечный ледяной покров, движущаяся нагрузка, произвольная замкнутая траектория, 
переменная скорость. 
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Введение. В последнее время, в связи с освоением северных территорий и районов Мирового океана, огромное 
внимание уделяется проблемам изучения поведения ледяного покрова под действием внешней нагрузки. Этим 
вопросам посвящено большое количество работ отечественных и зарубежных ученых. Воздействие на 
вязкоупругую плавающую пластину импульсной движущейся нагрузки рассматривалось в [1]. Воздействие 
подвижной нагрузки на ледяной покров, примороженный к стенкам канала, представлено в [2, 3]. Колебания 
ледяного покрова, обусловленные движущейся с постоянной скоростью нагрузкой, рассматривались в [4]. 
Поведение полубесконечного ледяного покрова при действии на него равномерно движущейся нагрузки 
изучалось в [5]. Распространения волн, возбуждаемых вдоль канала с ледяным покровом, рассмотрены в 
работе [6]. В [7] при моделировании применялись нелинейные модели. В монографии [8] приведены результаты 
исследований поверхностных волн в море с плавающим битым и сплошным ледяным покровом. Разрушение 
льда под действием движущейся нагрузки рассматривалось в [9]. В приведенных выше работах в основном 
рассматривалось прямолинейное движение нагрузки. Целью данного исследования являлась разработка метода 
решения задачи о воздействии силы, движущейся по ледяному покрову произвольным образом. Достижение 
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поставленной цели позволяет более точно моделировать воздействие на ледяной покров транспортных средств, 
движение которых происходит часто по достаточно сложным траекториям и по сложному закону. 

Постановка задачи. Рассматриваются колебания бесконечного ледяного покрова, лежащего на поверхности 
водоема конечной глубины, под действием силы, движущейся произвольно по замкнутой траектории. Ледяной 
покров моделируется тонкой вязкоупругой пластиной, механические свойства которой описываются моделью 
Кельвина-Фойгта. Водоем наполнен несжимаемой жидкостью. 

Материалы и методы. Изгиб ледяного покрова описывается дифференциальным уравнением [7]: 

(1+ тои, + ЗОРИ + КИ! + Б9, Фо = 9, (1) 
где И/(х,у,ё) — прогиб ледяного покрова; д=Ей/(12(1—и?)); Е — модуль Юнга; и — коэффициент Пуассона; й — 
толщина ледяного покрова; т, — время релаксации деформаций; Д’=0,*+20,? 9. +д,*; р„ — плотность льда; р, — 
плотность воды; с? = р.й/О, К=р,э/О, = р»/О; Р(х,у,ё) — действующая на поверхности льда нагрузка; Ф(х,у, 5,1) 
— потенциал движения жидкости. Поведение жидкости подчинено уравнению: 


ДАФ=0. 
Граничные условия на границе лед-вода при <=0 и на дне водоема <=-Н (Н— глубина водоема) имели вид: 
ОИ/=0.Ф 5=0, 0.Ф &=-Н =0. 


Предполагалось, что сила Р (х, у, #) перемещается по произвольной замкнутой траектории у произвольным 
образом. При этом полагалось, что Р=Р(5(1)), где 5— дуговая координата, отсчитываемая от некоторой 
фиксированной точки кривой 7. Если функция 5(1) является периодической с периодом Т, то вполне очевидно, 
что периодической будет и величина Р(5(!)), причем также с периодом Т. Параметрическое задание траектории 
принималось в виде: 

[ = хо(®) 
у = У@)’ 
где — время. 

Рассматривая установившийся процесс и применяя интегральное преобразование по переменным х иу, а 
также конечное интегральное преобразование по { на интервале [0; 7], получим: 


(1- Фо) ро с + ИНЬ Фо. -0=Рор, 
д2Фо - р? Фо=0, 
где \» Фо — изображения неизвестных функций УМ и Ф, р’=7+ а’; ^; а — параметры интегральных 
преобразований, соответственно, по переменным х и у; ®„=2ли/Г, п=0, 1,2, 3,.... Решая второе уравнение с 


учетом граничных условий, получим: 
Фо=-ю„Тосй(р(<=Н))/р5й(РН). 

Тогда ИЗ первого уравнения следует: 

М’, =Ро/(с?+ЬсИ(рН р) в" -р*-К+Е ,т,р*). 
Движущаяся сосредоточенная сила аппроксимировалась функцией: 

Р(ху)==”ехр(-=”((х-хь (и) +(у-у,()) Ид, 
где = — параметр. 

Далее, выполняя вполне очевидные преобразования, можно получить: 


О О: п„со8 (в, (#—т)) „эт (о, (#—т)) а 
вт Г р -+К ар. я + — я 7, (рК)е р атар 


п п 


00 
"и=(с?+Ьсй (РНУ/р) в?-р*-К , г=оитор*, В? =6?+ В, 
б=х,(т)-х, В= уо(т)У. 

Зная М, можно по известным соотношениям определить формулы для вычисления компонентов вектора 
перемещения и тензора напряжений в любой точке ледяного покрова. 

При численной реализации предлагаемого метода возникает проблема, связанная с необходимостью 
вычисления интегралов, один из которых в качестве подынтегральной функции имеет сильно осциллирующую 
функцию при достаточно больших и, а второй является несобственным с бесконечным верхним пределом. 

Для вычисления интегралов от сильно осциллирующих функций применялась квадратурная формула, 
полученная на основе использования метода кубических сплайнов [11]: 


Ь М- ®х,1 ох 
Ре 1 е Н-е ' 
[его =-—, У, - (м,„-М,), 
а © 1= 7 
где А; — длины элементарных отрезков, на которые разбивается интервал [а; 5]; $(%) — аппроксимация 


Хх)кубическим сплайном, М; = $ (х). 
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со 
При вычислении несобственного интеграла использовалось приближенное соотношение т (фар = 

А о 
у ЕФ)ар при достаточно большом А. Значение А выбиралось таким, чтобы оценка допускаемой ошибки 


со < 
| ] ы Г (р)ар| была достаточно малой. 
Используя традиционные методы, можно получить необходимые оценки для вычисляемых величин. 
Например, для величины прогиба ледяного покрова данная оценка имеет вид: 


(с? +ьет( АН)! А) + А* (1+ ,л,) + 
(с? чьем (АН) А) — А*- К] +[в,л,А* ] 


А со 
\ = (фар + |, м(ф)ар. 
При вычислении суммы ряда ДЛЯ ускорения его сходимости применялся метод сигма-множителей Ланцоша. 


[ее] 


|№(2)4Ф 


А 


< 45?е^ = 


Результаты исследования. Расчеты проводились для случая действия единичной сосредоточенной силы, 
которая двигалась по замкнутой кривой, представленной на рис. 1: 


Г,:51 =. 
5(0 = ее 7» +— 
251(9 1) 2 


При этом полагалось, что толщина ледяного покрова В=0,25 м, модуль Юнга материала пластины 
Е=500 000 000 Н/м*, коэффициент Пуассона и=1/3, плотность льда р=900кг/м*, плотность воды р=1000 кг/м*, 


, “=т, 1=[0:Т]. 


глубина водоема Н=5м, &=2,5. Радиусы, определяющие форму траектории движения силы, принимались 
равными А,=15 м, А2=9 м, Ю;=3 м (рис. 1). 


У, м 


-4 -3 -2 -1 0 1 2 3 4 Х,м 
Рис. 1. Траектория движения сосредоточенной силы 


На рис. 2а приведено изменение прогиба ледяного покрова, а также максимальных значений нормальных 


напряжений (при 2=+5) 5х и 5у (рис. 6, в) и при скорости движения единичной сосредоточенной силы 
у =7,4022 м/с и касательном ускорении \=0. 


М, мх 108 2! 9х, НИМ? у, Нар 
5 | | 200 | 
200 | | | | 
0 о 2857 | 
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- -200 1 
ь ‚м -50 1 | 
/ Ра | У, м % 
маш и’ о 270. и. ЗВ 
й Хм. Хм -50 Х,м `-50 0 | 
а) 6) в) 


Рис. 2. Изменение прогиба И напряжений ледяного покрова: а — изменение прогиба; б — изменение напряжения 5х; 
в — изменение напряжения 5, 
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На указанных графиках красной точкой обозначено место приложения силы, а синим цветом — траектория 
движения. 

На рис. 3 представлена зависимость изменения максимального значения прогиба ледяного покрова от 
величины скорости движения силы (а) в момент {=Т/2. Положение силы на траектории в данный момент 
отмечено красной точкой (рис. 1), при этом касательное ускорение \=0. Приведено два графика, один из 
которых соответствует радиусу траектории В ‚=15 м (сплошная линия), а второй — радиусу В =5 м (пунктирная 
линия). На графике б на рис. 3 изображена зависимость максимального значения прогиба ледяного покрова от 
величины касательного ускорения в момент 1=Т в точке траектории (2;0) на рис.1 при скорости у=0. На 
графике на рис. 3 изображена зависимость максимального значения прогиба от величины времени релаксации 


деформации. 


У, мх 10“ УМ, мх 10* У, мх10? 

-1,0 

-0,5 
9 

0 -1,4 
-1,6 

=1,5 
-1,3 

5 10 У, м/с 2 4 УМ! м/с? 4 6 8 то, С 
а) 6) в) 


Рис. 3. Зависимости прогиба ледяного покрова от: а — от скорости; б — касательного ускорения; в — времени релаксации 


Поведение жидкости представлено на рис. 4, где изображено распределение вектора скоростей частиц 
жидкости по глубине водоема при {=Т/2, скорости движения силы у=3,7011 м/с и ускорении \=0. 


0 
50 -20 


Х, м 


Рис. 4. Распределение вектора скоростей частиц жидкости по глубине водоема 


Обсуждение и заключения. Исследовано влияние кривизны траектории движения и механических свойств 
льда на прогиб ледяного покрова скорости и ускорения движения нагрузки. Расчеты показали, что прогиб 
ледяного покрова существенно зависит от скорости движения силы и ускорения движения. 

Однако характер качественного изменения перемещений и напряжений в ледяном покрове, обусловленных 
действием движущейся силы при изменении скорости и ускорения движения, менялся слабо. 

Заметное влияние на прогиб ледяного покрова оказывают механические свойства льда, в частности, время 
релаксации. 

Полученные результаты и предложенный метод решения подобных задач могут быть использованы при 
строительстве ледовых дорог или аэродромов на льду. Кроме того, предложенный метод решения показал свою 
эффективность и может быть использован при решении других подобных задач. 


А.В. Галабурдин. Метод решения задачи о движении нагрузки по ледяному покрову водоема по сложной траектории 
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Аннотация 

Введение. Данные о возникновении начальных разрушений, полученные испытаниями на стандартных 
образцах, не всегда можно экстраполировать на реальные сварные соединения и конструкции. Это обусловлено 
отличиями концентраторов в соединениях, т. к. после сварки возникает значительная структурно-механическая 
неоднородность зоны термического влияния и концентраторов напряжений. Протяженные, глубокие 
концентраторы рассматриваются как трещиноподобные дефекты, в вершинах которых образуется объемное, 
сложное напряженное состояние. Решается вопрос построения критических диаграмм начала возникновения 
предельного состояния в вершине концентратора, которое зависит от уровня внешней нагрузки и 
теоретического коэффициента концентрации. 

Материалы и методы. Для исследования напряженного состояния задействовали аналитические методы. 
Проанализирована литература по теме. Учтены особенности проверенных физических моделей и 
закономерности поведения материалов. Характеристики сплавов стали взяты из открытых источников и 
обобщены в виде таблицы. Нелинейные уравнения решались в прикладных программах МаЧафб. Построенные 
авторами диаграммы позволяют отследить корреляцию опасного уровня теоретического коэффициента 
концентрации напряжений и уровня внешней нагрузки. Для оформления графической части работы 
использовали Сигуе ЕИйп» Тооох МаЧаБ. 

Результаты исследования. Дана характеристика разрушений от концентраторов напряжений в сварных 
соединениях. Наглядно показано развитие трещин в зоне сплавления. Указаны условия, стимулирующие и 
тормозящие разрушение. Определен теоретический коэффициент концентрации напряжений @т. Показано, 
каким образом данный показатель зависит от ширины, высоты шва и от толщины свариваемой детали. 
Рассмотрены острые концентраторы напряжений с теоретическим коэффициентом концентрации @.. = 5...14 и 
более. Для этого случая приводится аппроксимирующая формула, которая учитывает максимальное 
напряжение в концентраторе в первом полуцикле, исходную деформацию и коэффициент асимметрии цикла 
нагружения. Через эти элементы задается показатель повышения максимальных напряжений в зависимости от 
числа циклов нагружения. Аналитически показаны условие текучести, напряженное состояние и коэффициент 
перенапряжения, учитывающий повышение первого главного напряжения для сложного напряженного 
состояния. Описана модель критического состояния в вершине острого макроконцентратора напряжений. Она 
представлена как зависимость относительных напряжений зарождения разрушения 02°/о%2 от концентратора. 
Проанализированы возможные вариации этой модели. Представлены зависимости относительных значений 
07° /оо от теоретического коэффициента концентрации @2° = ат. Для проверки физической адекватности 
данной модели построены графики, которые отражают изменения относительного напряжения внешней 
нагрузки при критическом состоянии в вершине концентратора напряжений. Обоснована неизбежность 
бифуркации как результата исследуемых процессов. Указаны два направления дальнейшего развития событий: 
хрупкое разрушение и потеря устойчивости напряженного состояния с переходом к росту пластических 
деформаций. Момент бифуркации определен как критическое состояние в очаге концентратора. 


СКА. Молоков, В.В. Новиков, М. Дабалез, 2023 
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Обсуждение и заключения. Анализ и расчеты, выполненные в рамках представленной научной работы, 
позволили, в частности, сделать выводы о роли ключевых факторов исследуемых процессов. Установлено, 
например, что работа сплава стали при высоком теоретическом коэффициенте концентрации напряжений 
зависит от характеристик напряженного состояния. При жестком состоянии возможно сдерживание сдвиговой 
деформации и наступление предельного состояния при меньшем значении теоретического коэффициента 
концентрации напряжений. При обычной прочности стали (в сравнении с высокой) фиксируется большее 
влияние объемности напряженного состояния на значение теоретического коэффициента концентрации 
напряжений. Вероятность разрушения зависит от сопротивляемости материала росту макротрещины. В 
будущих изысканиях возможно уточнение аналитических моделей и результатов, оценка эффективных 
коэффициентов концентрации напряжений. 


Ключевые слова: сварное соединение, теоретический коэффициент концентрации, дефекты сварных 
соединений, объемное напряженное состояние, концентрация напряжений, предел текучести, макротрещина. 


Благодарности. Авторы выражают признательность сотрудникам департамента морской техники и транспорта 
ДВФУ за консультативную помощь в разработке темы, а также рецензентам за ценные замечания по структуре 
рукописи. 


Для цитирования. Молоков К.А., Новиков В.В., Дабалез М. Оценка появления начальных разрушений от 
концентраторов напряжений в сварных соединениях и элементах конструкций. Айдуапсеа Епетеетия 
Кезеагсй (Кобюу-оп-Ооп). 2023;23(1):41-54. В рз://401.0то/10.23947/2687-1653-2023-23-1-41-54 


Опвта[ агие 


ЕуашаНоп оЁ Фе Оссоггепсе о? Пийа ЕаЙоге$ гот 59`е5$ Сопсепй`аюг$ ш У/@4ед ош 
апд Зетисига Еетен$ 


Коп$апеп А Моюкоу' ‚ Умегу У Муко" ©, Мопаттад рарае 
'Раг Еауеги Ке4ега! Ошхуегзну, 8, биКВапоуа 51., У|ауо$юкК, Ваз ап Недеганоп 
2\аауозюк Зе Ошщхуегзцу, 41, Сого]еуа 51., У1ауозюК, Ваз$1ап Недеганоп 

Ч Зрекиит01 1277 @ ятай.сот 


АЪ5егасЕ 

Гитодисйоп. Оаа оп Фе оссштепсе оЁ па аИагез оате4 @гоцэВ {е5Яп® оп зап4агА запр!ез саппо! а!\ауз Бе 
ехнаро!ае4 © геа| \е!4е4 ]о1{з ап згас@агез. ТЫ$ 15 дие ю Ше 41Шегепсе Бебмееп Ше сопсепгаог$ ш Ше 0115, 
Бесаизе аНег \е]ше Шеге 15 а уетНсапе зисв@ага| ап@ тпесвашса| Беегосепейу оЁ фе Веа-аНесе4 ап зше$$ 
сопсепгаюг топе. Ежепае4, Аеер сопсепга{юог$ аге сопз1Аеге4 аз сгасК-ПКе деЁес($, аё \Возе уегисез а уоитейлс, 
шшвах1а| $е5$ 5ае 15 Югтед. Те рарег аЧЧгеззез Ше 155ие оЁ сопзгасипе спйса| Фаотата$ оЁ Фе опзеё оЁ Фе 
Шило $бае аё Фе сопсепгаог уецех, \мсН 4ерепа$ оп Ше 1еуе] оЁ ежегпа! 1оа ап4 Ше Феогейса| сопсепганоп 
соеЁйслепи. 

Маепа5 апа Мейо4у. Апа!уйса| тефо4$ \еге изе4 ю за4у Фе з9е55$ 5ае. Тве Шегавте оп Фе юр1с \аз 
апа]у7е4. ТБе Геагез оЁ ргоуеп рБуз1са! плоде!5 ап4 райегиз оЁ Берау1ог оЁ таепа[5 \еге {аКеп шо ассоппе. ТВе 
спагас(ег1$Ис$ оЁ ее] аПоуз \еге {аКеп Нот ореп зоигсез ап4 зиплитанте ш а ‘абае4 Ююгт. МопПпеаг едиаНоп$ 
\еге зо]уе4 ш МаНаБ аррНсаНоп$. ТВе Фазтата$ сопзисе4 Бу Фе аи@ог5 епа Ме ю НасК Фе согеаНоп оЁ Фе 
Чапоегои$ 1еуе] оЁ Ше Шеогенса! $ие5$ сопсепигаНоп Гасог ап Ше 1еуе] оЁ ежегпа| 1оа4. Сигуе Еййпе Тооох 
Майаь уаз изе ю 4ез1еп Ше отарШс рагё оЁ Те \огК. 

Вези5. Те сБагасет$Нс оЁ датазе Нот $4е5$ сопсепгаюг$ ш \е!4е4 ]оз уаз г1уеп. ТБе сгасК ргораганоп т 
Фе аз1оп 7опе \уаз зво\уп. Те соп41оп$ зитшайпе ап шЫ119п5 аезбгасвоп \уеге ш4сае4. ТВе Феогейса| зе5$ 
сопсепёгайоп Гасог ат \аз зресеч. [ уаз зВо\п Ном 1$ ша1саюг 4ерепае4 оп Ше ул, Фе Вет? оЁ Фе 5еат 
ап Фе ШсКпез$ оЁ Ше \е!4е4 раг. Асие 31е5$ сопсешгаюгз ИВ Шеогейса| сопсепгаНоп #ас!юог &. = 5... 14 апа 
тоге \еге зи Чеа. Рог 115 сазе, ап арргохипаНпе юга \аз э1уеп Фа юоК шю ассоцпЕ Фе тахипит $е$$ ш фе 
сопсепгаюг ш Ше ВгзЕ Ва суе, фе пиа1 деРогтайоп, ап4 е Тоа4 гайо. ТЬгоиеН позе @етеп, ап шФсаюог оРап 
шсгеазе ш тахипит $е5е5 \аз зе! аереп4те оп Фе питфег оЁ Тюа4тз сусез. ТВе Но\ соп оп, Ше згез$ зае, 
ап Ше оуегуоЦасе Расог, \усВ юоК шо ассоцие Фе шсгеазе ш Ше Вт рипс1ра! уо[азе Рог а соб те4 5#гез$ зе, 
уеге апа1!уйсаПу зВо\уп. А то4е| оЁ Ше си@са| 5ае ай Фе арех оЁап асще $е5$ тасго сопсепгаюг \'аз дезсгтфеа. Ё 
\’аз ргезете4 аз Ше Череп4епсе оЁ е гейануе зеззез оЁ Ше пиНаНоп оЁ дезгасНоп 02° / 00,2 оп Ше сопсешгаюг. 
Розе уамаНопз оЁ 15 то4е| уеге апа!уте4. ТЬе дереп4епсез$ оЁ тейануе уашез 00° / 002 оп Ше Шеогейса1 


сопсептгайоп {аспюг &?° = ат. \еге ргезете4. То спесК Ше рвузса! адеапасу оЁ #15 то4е|, огарбз \еге сопзнисе4 


К.А. Молоков и др. Оценка появления начальных разрушений от концентраторов напряжений в сварных соединениях 


Фа теЙесе4 сВапсез ш Фе ге]аНуе $е5$ оЁ Ше ежегпа]| 1оа4 аё а сийса| зайе аё Фе зез$ сопсешгаюг арех. ТВе 
шеупа Шу оЁ Ыигсайоп а$ а гезаЕ оЁ Фе зи е4 ргосез5ез \’аз уа|ааеч. Т\ууо дтесйопз оЁ ег деуе!ортепи оЁ еуеп5 
у\еге шФсаеа: ие дезиасНоп апа 105$ оЁ уабПиу оЁ Ше зшеззе4 зе ул Фе шапяйоп ю ап шстеазе ш р!азйс 
ЧеРогтаНоп$. ТВе потлепе оЁ Магсаноп У’аз 4ейпе4 аз а сийса!| 5е ш Ше Ююсиз оЁ Ше сопсепгаог. 

Обсизяопт апа Сопсшяоп. 'ТВе апа!уз1$ ап4 сасшайопз рефогтед ут Фе Натеууо!К оЁ Фе ргезешеа зслепйЙс \уогК 
епаЫе4, ш рагасцаг, ю 4га\ сопса$1юп$ аБоцё Ше гое оЁ Кеу Расюгз оЁ Фе ргосеззез ипаег звиау. [ ууаз езбабИзВеа, Гог 
ехатре, фаЁ Ше орегайоп оЁ а 5ее! аПоу а а Шей ШФеогейса 5ез$ сопсепганоп Фасюг дереп4е4 оп Ше сНагасцен$Ис$ оЁ 
Фе $е5$ зе. ш а па зе, п ууаз роззе ю шШЫи 5пеаг аеогтпайоп ап4 Ше опзеё оЁ Ше Ниле зе аЁа 1о\уег уаше 
оЁРШе Шеогенса! зез$ сопсепагайоп Расюог. \/ Фе изиа] змепо@й о $ее] (ш сотрайзоп © ШВ), а отежег ппрасё оЁ Фе 
уоате оЁ Ше 5е55 зе оп Фе уаше оЁ Ше Шеогенса! 5'е5$ сопсепёганоп Ёас{юог \’аз гесог4еа. Те ргоБаБшу оЁ Пите 
Дерепае4 оп Фе гез1${апсе оЁ пе таепа] ю Ше ото\/В оЁ а тасгосгаск. ш Вибиге гезеагсВ, 1 1$ розз1Ые ю гейпе апа!уйса! 
п1о4е[5 ап4 гезиз, еуащае еНесНуе зге5$ сопсепганоп асюгз. 


Кеуумога$: мее4 ]ошё ШФеогенса! сопсептайоп Расюг, 4еесз оЁ \ее4 ]о5, уопитейлс $е$$ зе, 5е5$ 
сопсепгайоп, у зтепе, тастосгаск. 


Аскпо\еазетет 5. ТБе аифогз \ошША ПКе о ФапК Фе з{аЁ оЁ Фе ПерагитепЕ оЁ Магте Епотееппе апа Тгапзрог, 
ЕЕЕО, Юг Фе! ау1зогу а551$апсе 11 Ше гезеагсв шо Фе юр:с, аз \еП аз 10 Ше геулеуегз Гог уащае соплтеп$ оп 
фе зигасвеаге оЁ Фе тапизсре. 


Еог сНаНоп. МооКоу КА, МоущКоу УУ, Рабе7 М. ЕуашаНоп оЁ Фе Осслитепсе оЁ Пийа! РаПагез ош $ез$ 
Сопсепига юг ш \е4е4 Люшб ап Зиасига! Еетепв. Адуаисе4а Епвтеептиз Кезеагсй (Козюу-оп-Ооп). 
2023;23(1):41-54. В&рз://А01.0г®/10.23947/2687-1653-2023-23-1-41-54 


Введение. Статическая и усталостная прочность снижается из-за дефектов сварных соединений. Это могут быть: 

— трещины, образовавшиеся во время или после сварки; 

— концентраторы напряжений (подрезы, несплавления, поры, сварочные кратеры, высокая чешуйчатость, 
резкое изменение формы шва, наплывы и др.). 

В последнем случае определяющими будут: 

— форма, размер концентратора и положение в сварном соединении; 

— показатель напряженного состояния в вершине [1]. 

При циклических нагрузках в сварных соединениях часто возникают и развиваются трещины, 
спровоцированные концентраторами напряжений [2]. Риск зависит еще и от того, насколько напряженное 
состояние сварного соединения с концентратором близко к возникновению макротрещины и распространению 
ее при рабочих он и циклических нагрузках о_,. Очевидно, что последнее существенно зависит от размеров, 
формы макродефекта и напряженного состояния в его вершине. Например, продолговатый и узкий 
макроконцентратор напряжений опаснее, чем округлый. Вероятно, в его вершине сразу образуется трещина, и 
вместе с ней исходный концентратор будет представлять одну длинную макротрещину [3]. Величина такой 
«общей» макротрещины может оказаться критической, что снизит статическую прочность [4]. Это зависит от ее 
длины и механических характеристик материала в вершине исходного концентратора. 

Цель исследования —Щ разработка аналитической модели оценки опасного уровня концентрации 
напряжений. Речь идет об образовании трещины в вершине концентратора, резком снижении несущей 
способности сварных соединений и деталей конструкций. 

Материалы и методы. Научные изыскания в рамках заявленной темы базировались на известных 
физических моделях и закономерностях поведения материалов. Проанализирована теоретическая и прикладная 
литература. Из источников извлечены иллюстративные и справочные материалы. 

Напряженное состояние исследовалось аналитическими методами. Выявлялась взаимосвязь жесткости 
напряженного состояния и величины теоретического коэффициента концентрации напряжений. Результаты 
представлены в виде диаграмм. Такой способ визуализации дает возможность отследить корреляцию опасного 
уровня внешней нагрузки и теоретического коэффициента концентрации напряжений. 

Для проверки результатов расчетов по физической адекватности авторской модели строились графики 
зависимостей для критического состояния. В работе использованы сведения о широко распространенных 
конструкционных сталях феррито-перлитного класса в состоянии поставки (сталь 10, 22К, 50, СтЗсп, 37ХНЗА, 
ЗОХГСА и др.). Их механические характеристики получены из открытых источников' и обобщены в таблице 1. 
Данные визуализировали в Сигуе ЕИйпе ТооФох Майа5. 


' Сергеев Н.Н., Сергеев А.Н. Механические свойства и внутреннее трение высокопрочных сталей в коррозионных средах. М., Вологда, 
2020. 431 с. 
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Характеристики сталей 


Марка стали о, МПа о., МПа т фк 
10 320 190 0,17 0,73 
15Г 410 245 0,148 0,55 
Ст3Зсп 450 270 0,16 0,71 
22К 540 310 0,16 0,69 
50 680 350 0,16 0,62 
10ХСНД 540 390 0,132 0,71 
37ХНЗА 1014 743 0,12, 0,6 
ЗОХГСА 1750 1360 0,09 0,44 


Таблица 1 


Для вывода формул задействовали математический аппарат. При решении нелинейных уравнений 
использовали пакет прикладных программ МаЧаБ. 


Результаты исследования 
Характеристика разрушений от концентраторов напряжений в сварных соединениях. С начала 2000-х 


годов исследуются концентрации напряжений в сварных соединениях и конструкциях применительно к 
промышленным задачам [2-3, 5-9]. 
Примеры разрушений сварных соединений из-за концентрации напряжений приведены на рис. 1. 


9) е) 


Рис. 1. Разрушения крестовых и стыковых соединений” от концентраторов напряжений: а, е — от зоны непроплавления, без 
разделки кромок; 6 — от подреза, с разделкой кромок; в — от резкого перехода в зоне сплавления; г — подрез; 
д — от концентратора, образованного ступенчатым переходом в стыке кромок [10] 
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? Молоков К.А., Новиков В.В., Турмов Г.П. Основы расчетного проектирования сварных конструкций. Владивосток, 2019. Т. 1. 204 с. 
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Разрушение крестовых соединений с большим непроплавлением кромок (обозначено стрелками на рис. 1 а) 
возможно при нагрузках: 

— статических (в результате охрупчивания зоны сплавления); 

— циклических (вызывают образование трещин в вершинах концентратора). 

В [5] анализируется распределение напряжений для стыковых, угловых, крестовых соединений и 
прерывистых связей судовых конструкций. В стыковых и крестовых сварных соединениях особое значение 
имеет область основного металла, примыкающая к сварному шву. Это самое слабое сечение, определяющее 
прочность соединения при переменных напряжениях. Следует отметить, что швы всегда достаточно длинны по 
сравнению с толщиной металла и концентратором напряжения (подрезом в стыковом или крестовом 
соединении). Это препятствует сдвигу в области концентратора и создает двухосное, а чаще многоосное 
напряженное состояние [11, 12]. При рассмотрении разрушения сварных швов с концентраторами (рис. 1 а, д, е) 
учтем концепцию возникновения пластической деформации в области вершины [5, 7]. Кроме того, отметим, 
что концентраторы примыкают к зоне сплавления, которая после сварки испытывает продольные остаточные 
напряжения растяжения. Они влияют на степень жесткости напряженного состояния в очаге 
концентратора [1, 13-14]. Все это препятствует прохождению сдвига в очаге концентратора и может 
способствовать возникновению минимальной хрупкой макротрещины при высоком теоретическом 
коэффициенте концентрации напряжений ©, даже от статической нагрузки. 

В [9] описаны показательные разрушения от сварных дефектов в трубопроводах большого диаметра. 
Установлено, что при длительной эксплуатации (более 20 лет) появляются трещины усталости в подрезах 
продольных сварных соединений труб. Трещины в подрезах в зоне сплавления, как правило, разветвляются и 
распространяются по двум направлениям (рис. 2): одна отдаляется от зоны сплавления, другая развивается по 
зоне сплавления (по шву). Это говорит о существенном влиянии механических, структурных характеристик 
материала” и механохимической неоднородности на процесс распространения начального разрушения [13, 14]. 
С другой стороны, разветвление трещины может затруднять процесс дальнейшего разрушения и переход 
трещины в магистральную. 


Рис. 2. Разрушение в сварном соединении трубы от концентрации напряжений в зоне сплавления. Сталь 14ХГС. Одинарной 
стрелкой показаны усталостные трещины; двойной — кромка разрыва трубы; штрихами — очертания сварного шва. 
Здесь Ь — ширина шва; й — высота шва; 6 — толщина свариваемой детали; 

у — угол перехода от основного металла к шву [9] 


Теоретический коэффициент концентрации напряжений в зоне сплавления можно определить по формуле, 
учитывающей эффект концентрации от неплавного перехода в этой зоне [5]: 


а, 1+ 1лА И 907, (1) 
где у = А[90 ехр(—р/ро) + В] — угол перехода от основного металла к шву; р — радиус перехода от 
основного металла к шву; ро =1мм; А = 0,94...0,17; В = 0,8; Б — ширина шва; № — высота шва; 
6 — толщина свариваемой детали (рис. 2). 

Однако данная зависимость не учитывает измененные механические характеристики материала. В этих 
зонах обычно присутствуют закалочные структуры, механические характеристики которых могут значительно 
(до 30 %) отличаться от характеристик шва и основного металла. 

Можно косвенно контролировать трещинообразование и напряженное состояние в области вершины 
концентратора напряжений через коэффициент концентрации их интенсивности. Так, при разрушении образца 
от усталости, в вершине концентратора образуется одна макротрещина, а в коррозионной среде — 
группа (рис. 3). 


3 Молоков К. А. Оценка поврежденности феррито-перлитных сталей в условиях малоциклового нагружения. В: Наука. Инновации. 
Техника и технологии: проблемы, достижения и перспективы: сб. тр. конф. Комсомольск-на-Амуре; 2015. С 125-128 ОВС: 
Ьрз://еНЬгагу га/{ет.азр?14=23752241&зе!4=46181945 (дата обращения: 31.10.2022). 
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а) 6) 


Рис. 3. Усталостное разрушение стали 45 от У-образного концентратора напряжений: а — на воздухе; 
б — вводе — коррозионной среде. Увеличение в 134 раза" 


После образования группы макротрещин (рис. 3 6) роль макроконцентратора снижается, и дальнейшее 
развитие разрушения может быть приостановлено либо заторможено. Это подтвердили эксперименты, которые 
свидетельствуют о росте усталостной прочности в образцах с концентрацией напряжений в коррозионной 
среде”. Такая группа трещин как бы разрыхляет материал, изменяет напряженное состояние над поверхностью 
и под ней в вершине макроконцентратора, что препятствует образованию лидирующей трещины, способной 
развиваться далее при меньшем уровне циклических нагрузок. 

Рассмотрим ситуацию на поверхности или в толще металла в массивных телах под напряжениями он на 
бесконечности при плосконапряженном состоянии. В этом случае целесообразно использовать решение 
Колосова и Инглиса для концентратора эллиптической формы: 


бр = 1+2 а/р, (2) 


где ри обтах — радиус кривизны и максимальная компонента напряжений на поверхности соответственно, в 
вершине выреза; а — большая полуось эллипса, или половина наибольшей длины области, перпендикулярной 
направлению поля внешней нагрузки. 

Зависимость (2) хорошо согласовывается с инженерной практикой. В соответствии с принятой концепцией 
теоретический коэффициент концентрации напряжений у выреза определяется глубиной выреза и радиусом 
кривизны контура в его вершине, но не зависит от формы контура [15]. 

Благодаря развитию информационных технологий оценка теоретического коэффициента концентрации 
напряжений не представляет больших трудностей. Результаты по расчетам концентрации напряжений и 
деформаций в упругой и упругопластической постановке задачи можно получить средствами 
конечноэлементного анализа САЕ [16-—18, 12]. 

Эксперименты показали, что при определенной концентрации напряжений, внешней статической нагрузке и 
упругопластической стадии материала в вершине концентратора на некотором расстоянии от поверхности 
возникала начальная трещина разрушения. Она быстро распространялась к оси, перпендикулярной 
направлению растяжения пластин с надрезами. Распространение видно на снимках (рис. 4). Процесс 
зафиксирован в пластинах при плоском напряженном состоянии и при плоской деформации. В последнем 
случае появление трещины в толще металла под поверхностью выявили с помощью ультразвукового 
дефектоскопа. 


а) 6) в) г) 


Рис. 4. Зарождение подповерхностной трещины в вершине концентратора при статическом растяжении: а — начало 
возникновения макротрещины; 6, в —увеличение размера вглубь и кповерхности; г — выход трещины на поверхность 09] 


* Молоков К. А. Оценка поврежденности феррито-перлитных сталей в условиях малоциклового нагружения. 
> Сергеев Н. Н., Сергеев А. Н. Механические свойства и внутреннее трение высокопрочных сталей в коррозионных средах. 
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Таким образом, при указанных состояниях и достаточно остром концентраторе напряжений образуется зона 
многоосного напряженного состояния, в которой возможно нарушение сплошности. Здесь будет выше уровень 
наступления текучести, а предел прочности в ослабленном месте из-за концентратора может увеличиться до 
20—60 % (это зависит от материала) [11]. 

Концентрация напряжений и условия напряженного состояния. Рассмотрим достаточно острые 
концентраторы напряжений с теоретическим коэффициентом концентрации @., = 5... 14 и более. Как правило, 
в роли таких концентраторов выступают упомянутые выше дефекты. В вершине концентратора напряженное 
состояние можно рассматривать таким же, как при плоской деформации. В [18] приводятся результаты 
аналитического решения и экспериментальные данные для острых концентраторов при плоской деформации, 
устанавливается область допустимых значений [ат] и радиус в вершине [р] для стали 09Г2С. Однако не 
описано влияние ее механических характеристик. К тому же не ясно, основаны ли аналитические решения на 
таких исходных параметрах 09Г2С, как коэффициент упрочнения т, предел текучести бу, предел прочности 
обв, критическое сужение при разрушении фх и др. Известно, что пластичность и прочность существенно 
влияют на уровень допустимой концентрации напряжений и возникновение хрупкого разрушения в вершине. 

Для упрочняющихся материалов циклические напряжения в вершине острого концентратора усугубляют 
возможность возникновения хрупкого разрушения, поэтому за критерий опасности примем такой ат, при 
котором вероятно возникновение хрупкого разрушения в условиях статической нагрузки. Повышение 
максимальных напряжений в зависимости от числа циклов нагружения М№М задается аппроксимирующей 
формулой: 

бтах = Ошах + А '18(М), (3) 
где отах — максимальное напряжение в концентраторе в первом полуцикле; А — коэффициент, зависящий от 
уровня исходной деформации и коэффициента асимметрии цикла нагружения в вершине концентратора. 

Принимаем очевидное допущение, что появление хрупкого разрушения может быть обусловлено 
следующими причинами: 

— усиление интенсивности трехмерного растяжения в непосредственной близости к вершине концентратора; 

— увеличение сопротивления пластической деформации материала вблизи вершины как последствие 
сопротивления сдвигу в острых концентраторах напряжения. 

Последнее обстоятельство оказывается тем сильнее, чем выше скорость нагружения и острее концентратор 
напряжений. 

Рассмотрим упругопластическое плоское напряжение в вершине концентратора и упругое состояние в 
сечении брутто: бтах = @тбн » 602, бн ® 002 Ибн < 002. Здесь бу2 — условный предел текучести материала 
при одноосном растяжении. При &.. = 7,3 жесткость напряженного состояния близка к состоянию для трещины. 
Коэффициент безопасности ТЕдр для трещиноподобного дефекта при плоской деформации незначительно (в 
переделах 10 %) отличается от Трдр для трещины, т. е. для такого состояния, когда имеет место приближенное 
равенство ТЕдр & ТрЕдр: 


1-(во/от ела. |1 (воз/лыт/во)* [1 (боз/пт/во)* 
{1 ботт/оту дай [1 (в/о) [1-(вол/оо 


Здесь Пт — коэффициент запаса по пределу текучести, оо — локальная прочность материала у вершины 


(4) 


трещины или острого концентратора. 

В ПП термин «локальный предел текучести» применяется к острым концентраторам напряжения, когда 
примыкающая к вершине область концентратора испытывает сложное напряженное состояние (СНС). 
Очевидно, что этим термином обозначается напряжение течения материала на контуре концентратора. На 
поверхности в его вершине будет плоское напряженное состояние, под поверхностью возникнет трехосное 
СНС. Там резко возрастет жесткость напряженного состояния, поэтому начальное разрушение сплошности 
образуется на некотором расстоянии от поверхности (рис. 4). Установлено, что начало наступления локальной 
текучести в зоне концентратора напряжения не совпадает с уровнем, определенным расчетами по критериям 
прочности, в частности по Мизесу — Губеру — Генки. При увеличении концентрации напряжений возрастает 
разница между: 

— величиной предела текучести обо, теоретически рассчитанной по этому критерию; 

— экспериментально определенной величиной напряжения локального течения материала 00 в области, 
примыкающей к вершине концентратора. 

Результаты экспериментов на плоских образцах разных сталей и сплавов [11,16] показывают, что 
отношение 052/002 хорошо аппроксимируется линейной зависимостью от ат. Тогда условие текучести можно 
записать в виде: 
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0; = 002 (0,9 + 0,14), (5) 
где о; — интенсивность напряжений по Мизесу. 

Авторы ограничились экспериментальными исследованиями до значений ат = 10, и испытания 
проводились только на плоских образцах с центрально расположенными концентраторами. При плоской 
деформации жесткость напряженного состояния оказывается несколько выше. Поэтому логично предположить, 
что текучесть при некоторых значениях ©.т наступала бы не раньше, чем при плоском напряженном состоянии. 
И вполне вероятно, что наклон прямой по уравнению (5) мог быть несколько иным [11]. Тем не менее, 
остановимся на этой зависимости. 

Для момента начала текучести предположим, что с переходом от упругих деформаций к 
упругопластическим не меняется отношение второго и третьего главного напряжения к первому 01 = бщах. Это 
доказано экспериментально. 

Приравняем к (5) интенсивность напряжений начала текучести по Мизесу согласно четвертой 
(энергетической) теории прочности и введем коэффициенты отношений между главными напряжениями при 
плоской деформации: 0: = 0;/Р; в. = 901; 03 = ит(1 + ада. [15]. Здесь ит = 0,5 — коэффициент Пуассона в 


пластической области, а р = коэффициент перенапряжения, учитывающий повышение первого главного 
напряжения в случае СНС. У полученного относительно 4 уравнения два решения. Первое — для 
растягивающей о2 компоненты, второе — ДЛЯ сжимающей. Растягивающая компонента 02 увеличивает 
жесткость напряженного СОСТОЯНИЯ В случае СНС. После преобразований И сокращений можем выразить ее так: 

1-2 (1+2). (6) 


Для предельного случая он = 6,2 можно пренебречь изменением радиуса в вершине концентратора и 
наступлением некоторой начальной глобальной текучести для сечения нетто. Тогда единственное значение &т, 
удовлетворяющее жесткости эквивалентного трещине напряженного состояния, равно 7,3. Это подтверждает 
полученное ранее примерное равенство (4). Для других значений, как видно из (6), пропорциональность 4 будет 
меняться в зависимости от значения © т. И он. 

Для трещины коэффициент О определяется из условия равновесия в упругопластической области и 


вычисляется для плоской деформации по формуле: 
(1+т) (1-2) 


ПЕ 8. (1) 


2 
где т — коэффициент степенного упрочнения; и — коэффициент Пуассона в упругой области. 


Данный коэффициент для феррито-перлитных сталей принимает значения от 0,22 до 0,26. 

Множитель перед О в формуле а = 1- 20/\3 для трещины составляет 1,156, а а для феррито-перлитных 
сталей — 0,73. Действительно, при подстановке @г = 7,3 (6) при бн = бо, получаем 4 = 0,73. Далее найдем @т, 
если 0о2/сон > 1, при котором реализуется такое же напряженное состояние, как при @т. = 7,3 ибо2/он — 1. 

В [20] исследовано влияние ат на эффективный коэффициент концентрации напряжений К., а для 
типичного алюминиевого сплава показано, что в диапазоне ат от 7 до 13 наблюдается глобальный экстремум 
значения К.. Вполне вероятно, что такое существенное максимальное значение может свидетельствовать о 
хрупком разрушении на начальных этапах циклического нагружения или несплошности в вершине 
концентратора. В данном случае ресурс будет зависеть только от дальнейшей способности материала 
препятствовать распространению макротрещины. 

Модель критического состояния в вершине острого макроконцентратора напряжений. Практически 
все дефекты сварного соединения трубопроводов, корпусных судовых конструкций и др. создают 
концентрацию напряжений. В [9] представлена аналитическая модель, разработанная на базе понятия 
«локальный коэффициент концентрации напряжений». Кроме того, получена простая зависимость на основе 
формулы Нейбера, связывающей теоретический обобщенный коэффициент концентрации и коэффициенты 
концентрации интенсивности упругих напряжений и деформаций. В итоге для критической величины 
пластичности в очаге концентратора получена следующая зависимость: 


2/(1+т) 
бвас 
в: = 189 _—, (8) 
где ас — критическое значение теоретического коэффициента концентрации напряжений; Ов — условный 


предел прочности; ЕЁ — модуль упругости; 71 — показатель степенного упрочнения. 

Зависимость (8) дает такие высокие значения, которые редко встречаются на практике. Например, при 
критической логарифмической деформации для стали 50 критический коэффициент концентрации — 26, а для 
высокопрочной стали 37ХНЗА с ов = 1014 МПа — более 18. Можно заключить, что такая пластичность до 
нарушения сплошности в первых циклах нагружения недостижима. Это не удивительно, так как известно, что 
величина предельной пластичности для образцов с концентраторами напряжений существенно ниже, чем для 


К.А. Молоков и др. Оценка появления начальных разрушений от концентраторов напряжений в сварных соединениях 


образцов без концентрации напряжений. Кроме того, &‹ должно определяться не только свойствами материала, 
но и условиями развития пластической деформации перед разрушением. Похожие результаты дает 
зависимость, используемая в [9]: 


_ УЁ'Еразр"Оразр 
ат = У Внебьь, 6) 


С. А. Куркин использовал ее для сравнительной оценки чувствительности материала к концентрации 
напряжений. Здесь ат характеризует концентратор, в котором трещина возникает при конкретной нагрузке он. 
Далее будем использовать для этого случая понятие «наиболее вероятное зарождение макротрещины в 
концентраторе» (Боги сгасК) а?°. Истинные деформацию и напряжение в момент возникновения несплошности 
в вершине концентратора обозначим соответственно Еразр И бразр. Их можно определять по общепринятым 
формулам Еразр = &кр = Ш[1/(1 — ф‚)] И брзр = 5отр = 0,(1 + 1,4ф,). Здесь 5»„р — истинные напряжения 
разрушения при одноосном напряженном состоянии. Учтем, что текучесть в концентраторе наступает 
согласно (5). Заменим он = 002 и 002 на 002. В итоге получим для критических состояний (ат = в) 
относительные напряжения зарождения разрушения 07° /во,2 от концентратора: 
под = аРРоонолай аз — 

Отметим, что при очень малом радиусе вершины концентратора (р < 104, где 4, — средний диаметр 
зерна) распространение трещины ограничивается при неизменном @т. Как показано в [6], это зависит от а, 
стали. Формулы (8) и (9) привлекательны для инженера тем, что позволяют сравнивать чувствительность 
материалов к концентрации напряжений и использовать исходные данные механических характеристик стали. 

Можно получить несколько иную модель, учитывающую СНС в области вершины макроконцентратора. 
Автор [12] использует силовой критерий хрупкого разрушения в материале конструкции 01 > 5отр, ГДЕ 5отр — 
напряжение нормального отрыва в приближении к одноосному напряженному состоянию. Важно, что данная 
характеристика материала не зависит от температуры его испытания. Используем этот критерий несколько 
иначе — для момента возникновения несплошности в области вершины концентратора напряжений. Полагаем, 
Что 5отр Достижимо в случае постоянства энергетического равенства при СНС и одноосном напряженном 
состоянии, т. е. при потере пластической устойчивости и переходе от объемного напряженного состояния 
(когда о; мало) к одноосному напряженному состоянию с неизбежным возрастанием деформаций. 

Используем степенную аппроксимацию диаграммы деформирования в рамках деформационной теории 
пластичности. Приравняем истинные напряжения о; диаграммы деформирования в функции упругого решения 


задачи о; =Е а ) к 5отр. В результате получим: 
1 


тет 1 (@) 2 
9= 0% (00) =5ь = ов +14фк), (1) 
где в” —_ одноосные напряжения при упругом решении задачи. 


(у) 


Записывая (11) относительно о; 


(у) 


7 


‚ Получим выражение для упругого решения задачи о концентрации 


напряжений с другой стороны: 0; ° = ат ' он. Таким образом: 


а 
[ов +1 фк] | (12) 


1 


60,2 


Здесь он — Номинальная рабочая нагрузка на сварное соединение. В жестких точках сварных конструкций, 


“т ‘он = 


соединений или при перегрузках она может достигать значения предела текучести материала. Представим 
частный случай в условиях циклического нагружения с упрочнением материала для момента образования 
несплошности Богп сгаскК — зарождения разрушения. После преобразований (12) запишем упрощенную 


формулу: 


1+т 


абс = [есенею ®”, (3) 


60,2 
где &2° — значение теоретического коэффициента концентрации напряжений, при котором происходит 
зарождение разрушения в области вершины макроконцентратора; оу» — предел текучести материала в 
вершине концентратора, который может повышаться согласно (5) в случае СНС. 

В момент равенства ат = а2° номинальные напряжения должны быть эквиваленты напряжениям 
зарождения (разрушения) ой’. Для этой ситуации можно записать в относительных значениях 0 /бо2 и 
построить зависимости относительных значений напряжений образования разрушения от теоретического 
коэффициента концентрации &?< = ат: 
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=2т 
2-(«2°)"*1 ов фк+5) 


БС 
от +9 


БС — 


он = (14) 


К сожалению, не представляется возможным записать (14) относительно &?° при ой ". 002 = 1. Однако 0 
можно найти численно для конкретной нагрузки 07° или построить диаграммы опасности возникновения 
несплошностей в виде макротрещины в очаге концентратора. 

Аналитически получены зависимости для расчета критических теоретических коэффициентов концентрации 
напряжений, связанных с механическими характеристиками материала и заданной внешней статической 
нагрузкой при условии возникновения объемного напряженного состояния для концентраторов в сварных 
соединениях. Показано, что объемность напряженного состояния в очаге концентратора существенно влияет на 
величину критического теоретического коэффициента концентрации напряжений. Установлено, что сложное 
напряженное состояние в очаге концентратора может контролироваться геометрическими характеристиками 
самого концентратора и его расположением относительно внешнего поля напряжений. 

Экспериментальные данные сталей (таблица 1) получены из литературных источников. Это феррито- 
перлитные материалы, для которых формула (14) показала хорошее согласование. Зависимость (14) 
привлекательна тем, что в качестве исходных данных используются стандартные механические 
характеристики. Ее также можно применить для оценки опасности дефектов в зоне сплавления (рис. 1,2), где 
образуются закалочные структуры, свойства которых значительно отличаются от исходных характеристик 
свариваемой стали. 

Решения по (9) и (14) применяются чаще для самых хрупких сталей при невысоких &а?°. Однако в области 
высокой концентрации напряжений критические отношения 07° / 002 Для разных сталей будут весьма близки 
по значениям. Это можно объяснить тем, что при СНС в области вершины для всех сталей создается 
напряженное состояние, подобное опасному случаю. В области же малых концентраций напряжений решения 
по модели (9) показывают более высокие критические значения напряжений о#С. Отметим, что среднее 
значение 7‘ при 01° = оо2 составляет -8,5. Для разных по механическим характеристикам сталей оно 
меняется не так значительно, как по расчетным результатам модели (9). 

Рис. 5 построен для проверки результатов расчетов по физической адекватности модели (14). Это графики 
зависимостей критического состояния 0} /бо2 от а. 

БС 

он’ /от 
—= (13 
— 37ХНЗА 
—— Сталь 50 
—— Сталь 19 
— 22 
— ЗОХГСА 
—— Сталь 50 
—— 37ХНЗА 
== Ст 
—^ Сталь 19 


= 22К 
—— ЗО0ХГСА 


5 10 15 20 25 30 35 ос 


Рис. 5. Зависимости изменения относительного напряжения внешней нагрузки 07° /со2 от а? для критического состояния в 
вершине концентратора напряжений. Сплошные кривые рассчитаны по формуле 14; пунктирные — формула (10) 


Модель (9), не учитывающая объемность и повышение напряжений текучести в вершине концентратора, для 
стали 10 дает а = 91 при напряжениях 07° = 0,50%.2, что маловероятно. Если учитывать объемность в этой 
модели, то &2° = 33, что существенно меньше. При учете объемности напряженного состояния критическая 
пластическая деформация снижается до возникновения разрушения при статическом растяжении образцов. Это 
подтвердили эксперименты с различными сталями и сплавами°. Более того, на корродированных образцах 
исследователи получали меньшую пластичность сталей. При этом было показано, что в вершинах создается 
объемное напряженное состояние, и объемность препятствует свободному протеканию пластичности. 
Пластичность до разрушения этих образцов составляла -40 %, а предел текучести повышался до 27 % [21]. 
Можно заключить, что более глубокие концентраторы обусловливают большее повышение предела текучести. 


° Кроха В.А. Упрочнение металлов при холодной пластической деформации. Москва; 1980. 157 с. 
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Результаты косвенных опытов позволяют утверждать, что при гидростатической нагрузке текучесть не 
фиксируется, в этом случае материал разрушается только хрупко. Однако в очаге концентратора при 
увеличении внешней нагрузки объемное напряженное состояние не может оставаться достаточно стабильным 
для развития более высокого уровня объемности. Поэтому бифуркация неизбежна. Далее развитие может пойти 
по двум направлениям: 

— хрупкое разрушение или его зачатки в очаге; 

— потеря устойчивости напряженного состояния и, как следствие, переход к резкому возрастанию 
пластических деформаций. 

Вероятно, момент бифуркации следует считать критическим состоянием в очаге концентратора. Возможно, 
этот момент значительно зависит от глубины концентратора, его относительной длины и глубины к толщине 
детали, ориентации к внешнему полю напряжений, а также от механических характеристик и протяженности 
очага концентратора. Все это требует дополнительных исследований и более точных решений. 

Обсуждение и заключения. Дефекты и конструктивные особенности сварных соединений (форма разделки 
кромок, непроплавление, подрезы, поры, закалочные структуры и др.) снижают статическую и усталостную 
прочность. Для сварных соединений типичны концентраторы напряжений, имеющие длинную и неглубокую 
форму по отношению к толщине соединяемых деталей (подрезы) или длинную и глубокую (непроплавления 
и др.). Все они снижают усталостную прочность. Аналитические модели, полученные в представленной работе, 
позволяют сделать ряд заключений. 

1. Работа материала при высоких значениях теоретического коэффициента концентрации напряжений 
зависит от напряженного состояния и его жесткости, а также от механических и структурных характеристик. 
Последние могут отличаться от параметров исходного материала, так как вершины концентраторов могут 
находиться в зонах сплавления сварных соединений. 

2. Жесткое напряженное состояние в очаге концентратора напряжений может привести к повышению 
предела текучести. В этом случае прохождение сдвиговой деформации сдерживается, а наступление 
предельного состояния достигается при меньшем значении теоретического коэффициента концентрации 
напряжений и характеризуется несплошностью при постоянном значении внешней нагрузки. 

3. Сравнение высокопрочных сталей с пластичными и сталями обычной прочности при одном и том же 
уровне 07°/со2 дает следующие результаты. Во втором случае образование объемности напряженного 
состояния будет более существенно сказываться на изменении критического значения теоретического 
коэффициента концентрации напряжений, при котором происходит зарождение разрушения в области вершины 
макроконцентратора &2< (рис. 5). 

4. Сравнение моделе (14) и (10) позволяет сделать определенный вывод. С учетом жесткости напряженного 
состояния и повышения предела текучести при 017 /бо2 — 1 снизу (10) дает оценку &7° с запасом но прочности, 
если речь идет о ой / 00,2 = 0,8 (30ХГСА). Для менее прочной стали (14) запас по прочности при 0 И 00,2 = 0,5 
(37ХНЗА) ит. д. Таким образом, предпочтительнее учитывать объемность и использовать модель (14). 

5. Зависимости на диаграмме теоретических коэффициентов напряжений более плотно сгруппированы 
в (14), чем в (10) из-за очень похожих напряженных состояний в области вершин острых концентраторов 
напряжений. Поэтому механические характеристики сталей оказывают второстепенное влияние на &^-. 

6. При увеличении теоретического коэффициента концентрации напряжений и уменьшении &@?° при одном 
уровне о02’/оо2 внешней нагрузки снижается влияние структурного фактора в стали на возникновение 
начального разрушения (образование несплошности). Однако для дальнейшего разрушения существенное 
значение будет иметь сопротивляемость материала росту макротрещины. После ее возникновения срабатывают 
другие критерии механики и кинетики разрушения. 

Представленное исследование может служить предпосылкой к разработке аналитических моделей для 
оценки остаточного ресурса сварных соединений и конструкций, подверженных действию циклических 
нагрузок. Дальнейшие изыскания, предположительно, уточнят аналитические модели. Авторы проверят и 
апробируют результаты путем моделирования средствами САЕ и перейдут к оценке эффективных 
коэффициентов концентрации напряжений. 
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Аннотация 
Введение. Полиэтилен является наиболее широко применяемым материалом в различных областях народного 


хозяйства, а изделия из него обладают существенными достоинствами, такими как лёгкость, нерастворимость в 
органических растворах при вполне удовлетворительной прочности. Однако механизм его разрушения 
достаточно сложен и зависит от условий работы и контактирующих с ним веществ. Цель исследований 
заключалась в изучении устойчивости полиэтилена к нефти и сопутствующей воде при статическом и 
динамическом режиме и комнатной температуре. 

Материалы и методы. Исследования проводились на лабораторном устройстве для пропускания различных 
жидкостей (нефть, сопутствующая и дистиллированная вода) в полиэтиленовых трубах, собранном авторами 
статьи. При работе были применены методы статистической и динамической обработки данных. 

Результаты исследования. Графики изменения веса и объёма погруженных гранул во времени, построенные 
по результатам проведенных экспериментов, показали, что для полиэтилена характерна квадратичная 
зависимость и диффузия для трех исследуемых в данной работе жидкостей (дистиллированная вода, 
сопутствующая вода и нефть) описывается законом Фика. Это указывает на тот факт, что скорость диффузии 
жидкости через полиэтилен является наиболее важным фактором. 

Обсуждение и заключения. Эксперименты и полученные результаты показали, что скорость диффузии 
жидкости через полиэтилен является наиболее важным фактором. Погружение в нефть оказывает большее 
влияние, чем погружение в сопутствующую или дистиллированную воду из-за присутствия солей. Было 
обнаружено, что относительное изменение веса и толщины стенок полиэтиленовых труб, по которым проходит 
нефть, больше тех, по которым проходит сопутствующая вода, а микроскопические изображения срезов в 


образцах до и после испытаний подтвердили полученные результаты. 


Ключевые слова: полиэтилен, сопутствующая вода, твердые материалы, соли, дистиллированная вода. 
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Введение. Полиэтилен является одним из основных элементов семейства олефинов, который обычно 
характеризуется воскообразным внешним видом, химической инертностью и ухудшением физических свойств 
при высоких температурах. 

Поскольку полиэтилен является частично кристаллизованным материалом и не вступает в побочные 
химические реакции с какими-либо жидкостями, он не имеет растворителей при температуре окружающей 
среды и водопоглощение сравнительно невелико. 

Полиэтилен не растворяется в органических растворах, но он расширяется при температуре более 70 °С и 
может растворяться в гранулированных ароматических углеводородах, не подвержен влиянию нефти, жиров, 
керосина и других углеводородов, полученных из нефти. Он также устойчив к водным растворам кислот, солей 
и щелочей [1-5]. 

Полиэтилен окисляется азотной кислотой, что приводит к ухудшению характеристик, получаемых из него 
продуктов. Окисление происходит под воздействием ультрафиолетовых лучей при температуре окружающей 
среды, вызывая хрупкость и ухудшение механических и теплоизоляционных свойств [6-8]. Научный интерес 
представляет изучение взаимодействия полиэтилена с нефтью. Она представляет собой смесь углеводородов и 
неметаллических элементов, смешанных с ними (сера, азот, кислород), а также некоторых металлических 
элементов (железо, никель, мышьяк, ванадий). Наряду с ними в нефти присутствуют соли хлоридов земельных 
минералов, а также их сульфаты. Кроме того, нефть содержит воду в определённых процентных соотношениях, 
зависящих от ее типа и условий добычи. Этот процент может возрастать до 30 % по отношению к общему 
объёму нефти, поскольку она находится в форме насыщенных растворов хлорида натрия, сульфата магния и 
кальция. Двуокись углерода и газы сероводорода, растворяясь в воде, увеличивают её объём. Кислота и 
кислород, вступая в реакцию с железом и марганцем, растворенными в сопутствующей воде, образуют 
нерастворимые продукты. Сопутствующая вода содержит небольшой процент хлоридорганических солей, не 
превышающий 50 частей из миллиона [9, 10]. 
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Примеси, содержащиеся в сырой нефти (твердые материалы, соли, вода, химические добавки), создают 
серьезные риски для нефтеперерабатывающего оборудования, продолжительности его работы и 
эффективности. Известно, что полимерные материалы разрушаются совершенно иначе, нежели металлические, 
поэтому при производстве из них различных изделий требуется учитывать их физико-химические свойства и 
электрохимические реакции. 

Пластики не растворяются, как металлы, однако повреждаются или разрушаются вследствие набухания, то 
есть происходит так называемое охрупчивание, при котором они размягчаются, затвердевают и меняют цвет, 
что приводит к ухудшению их механических свойств. 

Полимерные материалы разрушаются по нескольким механизмам: 

— набухание структуры; 

— растворение; 

— разрушение химических связей под воздействием тепла, химических реакций или реакции свободных 
радикалов. Охрупчивание поверхностного слоя может происходить по какому-либо одному механизму или по 
комбинации механизмов [11-13]. 

Растворение и набухание без разрыва химических связей является основной причиной охрупчивания при 
контакте с примесями. Этот процесс является сложным, поскольку небольшие количества полимерного 
раствора проникают внутрь [14], образуя просачивающийся слой с измененными физическими 
свойствами (рис. 1). 


Растворитель 


Слой жидкости 


Гелевый слой 


Растворитель 
Твердый вспушенный слой 
Полимер а у 


Проницаемый слой 


Полимер 


Рис. 1. Слоистая структура процесса охрупчивания 


Целью исследования является изучение устойчивости полиэтилена к нефти и сопутствующей воде при 
комнатной температуре, как в состоянии покоя, так и в условиях движения. 

Материалы и методы. В таблице 1 показаны использованные в работе вещества, растворенные в нефти 
после её добычи. 


Таблица 1 
Вещества, растворенные и взвешенные в нефти после её добычи 
Растворимые вещества Взвешенные вещества 
Металлические Органические 
Твердые Газообразные 
материалы материалы 
Положительные Ее(ОН)з, Соединения, 
электролиты: Ме*, Ее5», полученные в 
Ее", К’, Ма", Ми" 5СО> результате добавления 
Отрицательные Растительные 0. Н5 м СО песок с глиняной деэмульгаторов, 
электролиты: Вх, и животные о крошкой, сульфатными бактерий, водных 
СГ, Е, СО;”, соединениями, растений, а также 
50: Н СО; карбонатами, щелочно- растительных и 
земельными минералами животных остатков 
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В таблице 2 показан состав сопутствующей ВОДЫ, использованной в экспериментах. 


Таблица 2 


Состав сопутствующей ВОДЫ, использованной в экспериментах 


Анализ Измеренное значение Единицы измерения 
Плотность 1,015 т/см? 
СГ 10100 
Ма* 2427 
Ме', 370 
Са”* 1270 
Полная жестокость СаСО} 4243 
Щелочность (СаСоз) 530 ррт 
50° 1093 
МОз 100 
МНу 13 
Н25 130 
РН 8,5 


В качестве чистящего средства применялся нефтяной спирт. Опыты проводили с использованием гранул и 


труб из полиэтилена высокого давления НОРЕ низкой плотности (таблица 3). 
Таблица 3 


Технические характеристики полиэтиленовых гранул и труб, использованных в экспериментах 


т Гранулы Трубы 
ип 
Полиэтилен высокого давления низкой плотности 
Плотность, г/ смз 0,914 0,917 

Твердость 44.6 467 
(по Шору) 

Отношение удлинения к срезу, % нет 4,42 

м 
инимальная температура | 109—122 110—112 
наступления размягчения, С 


Нефть и сопутствующая вода проходили по полиэтиленовым трубам в собранном автором лабораторном 
устройстве (рис. 2). Комплекс оборудования включает герметичную стеклянную тару; шкалу с 
чувствительностью 0,1 г, 9%; электронный микрометр с точностью до 1 микрометра; электронный микроскоп, 


подключенный к компьютеру. 


Полиэтиленовая труба 


Полиэтилен 


Нефть с 


Емкость 


Рис. 2. Лабораторное устройство для пропускания нефти и сопутствующей воды в полиэтиленовых трубах 
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В спокойном режиме при методе погружения последовательность действий при постановке эксперимента 
была следующей: 

— масса \/, и объем У, определяли для нескольких образцов гранул полиэтилена, характеристики которых 
приведены в таблице 3; 

— образцы погружали полностью в нефть, попутную и дистиллированную воду в герметичных контейнерах 
на определенные периоды времени; 

— далее образцы очищали; 

— в процессе эксперимента определяли вес образца \\/› и его объем У; 


— далее рассчитывали изменение веса А\У/=\/,-\ | и изменение объема ЛУ = У.-—У,, 


ДИ 
— После тестирования каждого образца рассчитывали: относительное изменение веса Е 100 и 
1 


ду 
относительное изменение объема т 100. 
1 


В динамическом режиме последовательность действия при постановке эксперимента была следующей: 

— нефть и сопутствующая вода протекали в полиэтиленовых трубах с постоянной скоростью при комнатной 
температуре, их вес М, и толщина стенки определялись для конкретных периодов времени; 

— После каждого опыта трубы очищали; 

— После испытания определяли их вес М”, и толщину стенок после испытания х,; 


— рассчитывали: изменение веса Д\/ = \/›\!:, изменение объема ДУ = \У›-У,, относительное изменение 
ДИ 
веса `` 100 и изменение толщины АХ после испытания. 
1 
Результаты исследования. Результаты испытаний изменения массы и объема частиц, погруженных в 
нефть, сопутствующую и дистиллированную воду, были сведены в таблицы 4—9. Графическая зависимость 
была найдена для относительного изменения веса и размера полиэтиленовых гранул, погруженных в нефть, 


сопутствующую и дистиллированную воду (рис. 3, 4). 


Таблица 4 
Изменение веса погруженных в нефть гранул 
№ 1 2 3 

Время 

погружения, У, г У» г ДУ, % Уьг У, г ДУ, % Уьг У» г ДУ, % 

дни 

10 27,022 27,673 2,240 26,9611 | 27,5542 | 2,200 25,384 | 25,921 2,261 

30 26,652 | 28,6376 | 7,450 26,4489 | 28,4589 | 7,600 26,577 | 28,613 7,661 

60 25,839 | 28,2497 | 9,330 23,486 25,7346 | 9,574 25,266 | 25,266 9,396 

90 26,863 | 29,6979 | 10,55 27,3432 | 30,1732 | 10,350 | 25,683 | 25,683 10,28 

Таблица 5 
Изменение объема погруженных в нефть гранул 
№ 1 2 3 
Время 
погружения, И У» м | ДУ, % и. а: АДУ, % Ты, ба ДУ, % 

Дни 

10 50 51,0 2,0 50 51,0 2,0 50 51,0 2,0 

30 50 52,1 4,2 50 52,4 4,8 50 52,3 4,6 

60 50 53,0 6,0 50 53,0 6,0 50 53,0 6,0 

90 50 53,7 7,4 50 53,7 7,4 50 53,9 7,8 
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Таблица 6 
Изменение веса гранул, погруженных в сопутствующую воду 
№ 1 2 3 
Время 
погружения, У, г И г ДУ, % У, г У г | ЛИ, % У, г И г ДУ, % 
дни 
10 55,045 | 55,430 0,700 58,765 58,765 0,70 28,514 28,69 0,63 
30 55,333 | 55,948 1,110 26,400 | 26,400 1,10 26,962 27,25 1,07 
60 45,831 | 55,629 1,450 27,332 | 27,332 1,60 26,614 26,99 1,44 
90 26,226 | 26,645 1,595 27,085 27,085 1,65 26,843 27,25 1,53 
Таблица 7 
Изменение объема частиц, погруженных в сопутствующую воду 
№ 1 2 3 
Время 
погружения, ие. Ур м | АУ, % А. У» м | АУ, % а ры ДУ, % 
дни 
10 50 50 0 50 50 0 50 50 0 
30 50 50,5 1 50 50,5 1 50 50,5 1 
60 50 51 2 50 51 2 50 51 2 
90 50 51,52 2,5 50 51,1 2,2 50 51,1 2.2 
Таблица 8 
Изменение веса гранул, погруженных в дистиллированную воду 
№ 1 2 3 
Время 
погружения, У, г И г ДУ, % У, г И г ДУ, % У, г И г ДУ, % 
дни 
10 56,003 | 56,003 0,000 26,761 26,762 0,00 26,723 26,723 0,00 
30 54,40 54,40 0,000 27,058 27,058 0,12 26,599 26,599 0,00 
60 55,48 55,489 | 0,000 26,303 26,303 0,00 26,084 26,084 0,0 
90 55,37 55,37 0,000 26,8300 | 26,8300 | 0,000 27,1665 | 27,1665 0,000 
Таблица 9 
Изменение объема гранул, погруженных в дистиллированную воду 
№ 1 2 3 
Время ду 
погружения, | Увы” | Уьы | АУ% | Увы | Узы | | Увы | Ум | А\% 
дни 
10 50 50 0 50 50 0 50 50 0 
30 50 50,1 0,2 50 50,1 0,2 50 50,1 0,2 
60 50 50 0 50 50 0 50 50 0 
90 50 50 0 50 50 0 50 50 0 
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Рис. 3. Относительное изменение веса погруженных гранул по времени 
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Рис. 4. Относительное изменение объема погруженных 


Рис.3 и4 показывают, что временная зависимость относительного изменения веса и размера гранул 
полиэтилена имеет характер квадратичной зависимости и следует закону Фика в диффузии в случаях трех 
исследуемых жидкостей (дистиллированная вода, сопутствующая вода и нефть), что указывает на то, что 
скорость диффузии жидкости через полиэтилен является наиболее важным фактором. 

Относительное изменение веса гранул, погруженных в дистиллированную воду, не упоминается, но оно 
едва заметно в сопутствующей воде и больше увеличивается в нефти (рис. 4). Это объясняется тем, что 
диффузия солей оказывает существенное влияние на относительный вес и объем гранул, погруженных в 
дистиллированную воду, с одной стороны, и гранул, погруженных в нефть и сопутствующую воду, с другой. 

Результаты испытаний изменения веса и толщины труб, в которых течет нефть и сопутствующая вода, 
систематизированы в таблицах 10 и 11. Была получена графическая зависимость относительного изменения 
веса и толщины труб, по которым проходит нефть с сопутствующим потоком воды, показанная на рис. 5, 6. 

Микрофотографии срезов образцов были сделаны до и после испытаний с увеличением в 510 раз, показаны 


на рис. 7. 
Таблица 10 
Изменение веса и толщины труб, по которым течет нефть 

Время погружения, дни 30 60 90 
№ трубы У, г Х,, мм У, % АХ ДУ, % АХ У, % АХ 
1 18,019 2,483 0,6 0,024 1,07 0,05 3,02 0,17 
18,6440 2,485 0,9 0,032 1,12 0,07 3,01 0,19 
3 18,6009 2,490 0,58 0,022 1,14 0,07 3,12 0,21 
4 18,6441 2,483 0,91 0,032 1,12 0,096 3,00 0,19 
5 18,5999 2,487 0,37 0,022 1,13 0,08 3,11 0,22 
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Таблица 11 
Изменение веса и толщины труб, по которым течет сопутствующая вода 
Время погружения, дней 30 60 90 
№ трубы У,: Х,, мм у, % АХ у, % АХ у, % АХ 

1 15,22 2,483 0,399 0,012 0,521 0,025 1,35 0,087 
2 19,134 2,484 0,404 0,015 0,478 0,029 1,42 0,086 
3 15,905 2,483 0,409 0,011 0,459 0,028 1,38 0,075 
4 18,6009 2,400 0,39 0,013 0,469 0,027 1,40 0,080 
5 18,019 2,483 0,38 0,012 0,498 0,027 1,39 0,079 
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Рис. 5. Относительное изменение веса труб, по которым течет нефть и сопутствующая вода, по времени 
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Рис. 6. Относительное изменение ТОЛЩИНЫ труб, в которых течет нефть и сопутствующая вода, в зависимости от времени 
Рис.би7 показывают, что относительное изменение веса и толщины стенок полиэтиленовых труб, по 


которым проходит нефть, больше тех, через которые проходит сопутствующая вода. Этот факт указывает на то, 


что полиэтилен больше подвержен влиянию нефти. 
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Рис. 7. Микрофотографии срезов образцов до и после испытаний с увеличением в 510 раз: а — сечение чистой гранулы; 
б — сечение гранулы нефти; в — сечение трубы с сопутствующей водой в нефти; г — сечение чистой воды; д — сечение 
трубы, погруженной в нефть; е — внутренняя поверхность трубы с сопутствующей водой; ж — сечение трубы, 
погруженной в нефть после 200 часов; з — внутренняя поверхность трубы, погруженной в нефть 


На рис. 7 показаны срезы образцов труб и гранул до и после испытаний, отражающие цвет и гладкость 
поверхности срезов, подвергнутых испытаниям, по сравнению с чистыми срезами, что подтверждает 
предыдущие результаты, полученные на основе графиков. 

Обсуждение и заключения. Эксперименты и полученные результаты показали, что скорость диффузии 
жидкости через полиэтилен является наиболее важным фактором при статическом и динамическом режимах. 
Погружение в нефть оказывает большее влияние, чем погружение в сопутствующую или дистиллированную 
воду из-за присутствия солей. 

Было обнаружено, что относительное изменение веса и толщины стенок полиэтиленовых труб, по которым 
проходит сопутствующая вода меньше, чем тех, по которым проходит нефть, а микроскопические изображения 
срезов в образцах до и после испытаний подтвердили полученные результаты. 
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Аннотация 

Введение. Электронные устройства, способные собирать данные по телеметрии индивидуума, открыли 
перспективы доклинического выявления признаков СОУТ-19. Известные решения предполагают анализ 
информации, которую сложно получить в моменте. Речь идет, например, о состоянии крови или ПЦР-тесте. Это 
существенно ограничивает возможности интеграции алгоритмов с наручными гаджетами. При этом сердечно- 
сосудистая система как объект наблюдения достаточно информативна, съем данных хорошо проработан. В 
статье описана задача детекции ковидных аномалий в ритмограммах. Цель работы — создание математической 
модели на базе алгоритмов машинного обучения для автоматизации процесса выявления ковидных аномалий в 
ритме сердца. Показана возможность интеграции полученных результатов с фитнесс-браслетами и умными 
часами. 

Материалы и методы. В работе задействовали открытый стек технологий: Рушоп, ЗсИ-еаги, Газ ебт. При 
оценке качества моделей для бинарной классификации использовалась метрика Ё\. Изучены 229 ритмограмм 
сердца (кардиоинтервалографий) пациентов с СОУТ-19. Наличие или отсутствие признаков аномалии 
определялось с учетом времени ритмограммы и интервалов между сердцебиениями. Графически показаны 
отклонения, которые могут свидетельствовать о заражении. По итогам разведочного анализа собран перечень 
признаков, указывающих на аномалию. 

Результаты исследования. В результате проделанной работы получена математическая модель, которая 
детектирует специфичные для СОУПТ-19 аномалии сердечного ритма с точностью 83 %. Выявлены и 
ранжированы основные признаки, определяющие прогностическую способность модели. Это текущее значение 
интервала между ударами сердца, производные в последующей и предыдущей точках измерения 
продолжительности сердцебиения, первая производная в текущей точке и отклонение от медианы текущего 
значения длительности АРА-интервала. Первый показатель в этом перечне признан наиболее значимым, 
последний —Й наименее. Для целей машинного обучения оценивался потенциал пяти алгоритмов: 
[зо]анопРогез, ГОВМОазяШег, КапаотРоге{СЛаз1Йег, Ехиа ТгеезСазвег, 5ОРОпеС]аз;5 УМ. 
Визуализированы нормальные и аномальные результаты наблюдений в изолирующих деревьях. Установлен 
параметр, который соответствует вероятности регулярного наблюдения за пределами нормы, и выбрано его 
значение — 0,11. С учетом данного показателя построен график для модели 5СПОпеСаз$5УМ. По набору 
данных с применением техники перекрестной проверки рассчитана метрика качества. Речь идет о ритмограмме 
с временным рядом наблюдений, снятых за один непрерывный интервал времени у одного человека. Описан 
пошаговый процесс получения усредненных значений метрики для каждой модели. При сравнении самый 
высокий показатель зафиксирован у модели ГОВМСИазуШег, наименьшие — у 5СООпеаз5УМ и 
[зо]анопРогез. 

Обсуждение и заключения. Полученная математическая модель занимает мало места в памяти мобильного 
устройства, то есть не предъявляет значимых требований к вычислительным ресурсам. Решение обладает 
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приемлемым качеством детекции для доклинического скрининга связанных с СОУТ-19 сердечно-сосудистых 
нарушений. Алгоритм обнаруживает аномалии в 83 % случаев. Для записи ритмограммы достаточно 4 минут. 
Предлагаемый сценарий использования интегрированного решения лаконичен и легко реализуем. Широкое 
использование разработки может способствовать выявлению СОУТО-19 на ранней стадии. 


Ключевые слова: СОУТШ-19, причины смерти ковид-положительных пациентов, осложнения в работе 
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Введение. Изучение влияния СОУТШ-19 на человека остается актуальной задачей. Так, в 2021-2022 гг. по 
данной теме опубликовано более 16 тыс. научных работ. Одна из основных причин смерти ковид- 
положительных пациентов — осложнения в работе сердечно-сосудистой системы (далее — ССС), вызванные 
воздействием коронавируса [1]. Для доклинической диагностики СОУПТ-19 в основном используются два 
метода: биохимический на основе полимеразной цепной реакции (ПЦР-тест) и анализ крови. Необходимые в 
данном случае контакты с медперсоналом (в том числе визиты в медицинские учреждения) затрудняют 
регулярный оперативный контроль и повышают нагрузку на систему здравоохранения. Таким образом, 
представляется актуальным применение современных технологий доклинического контроля ССС для раннего 
выявления признаков СОУГО-19. 

Регулярность контроля могут обеспечить носимые электронные устройства. Наиболее распространенные из 
них — фитнесс-браслеты и умные часы со встроенными датчиками частоты пульса и способностью выполнять 
измерения с высокой дискретностью [2]. Такой подход открывает возможности для анализа потоков данных на 
базе машинного обучения! [3]. 

Цель представленного исследования — создание обучаемой модели, способной выявлять ковидные 
аномалии, опираясь только на данные о ритме сердца. В ряде работ [4—6] рассматриваются подобные задачи, 
однако решения опираются на дополнительную информацию о состоянии крови и другие характеристики. Это 
существенно ограничивает возможности их интеграции с носимыми устройствами, т. к. в моменте невозможно 
ввести в модель результаты анализа крови или мазка для ПЦР-теста. Новизна предложенного решения состоит 
в том, что использовались только данные ритма сердца, которые можно с высокой частотой снимать удобным 
для человека способом и интерпретировать показатели в режиме реального времени. 

Материалы и методы 

Характеристика данных. В работе использовали 229 обезличенных ритмограмм (кардиоинтервалографий) 
пациентов с СОУТ-19. Сведения получены в 2021 году в рамках открытого всероссийского соревнования для 
профессионалов в сфере цифровой экономики «Цифровой прорыв». Фрагмент данных представлен в таблице 1. 


Таблица 1 
Фрагмент набора данных 


Номер Время ВК-интервал между Признак ковидной 
ритмограммы | в миллисекундах | сердцебиениями в миллисекундах аномалии* 
81 0 576 0 
81 568 568 0 
81 1140 572 0 


' Эндогенные аномалии кардиоритма у пациентов с СОУТ-19 / С. А. Пермяков [и др.] // Нелинейная динамика в когнитивных 
исследованиях — 2021 : тр. УП Всерос. конф. Нижний Новгород : Ин-т прикладной физики Российской академии наук, 2021. С. 109-110. 

Е 01а210$15 оЁ СОУТШ-19 апа 1$ сШшиса1 зресиат / Каззе Шшс. // Касе.сот : [сайт]. ОВЫ: 
Брз:/у\у м Каз е.сот/Ча!азе/ешуеш4аа4и/соу!1А19 (дата обращения : 10.09.2022). 
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Номер Время ВК-интервал между Признак ковидной 
ритмограммы | в миллисекундах | сердцебиениями в миллисекундах аномалии* 
176 44332 568 0 
176 44968 636 1 
176 45596 628 0 
*0 — аномалии нет, | — аномалия есть. 


На рис. 1 показана связь ритмограммы (КВ пиегуа1) с электрокардиограммой сердца (ЕСО). 
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Рис. 1. Сопоставление электрокардиограммы и ритмограммы сердца: по горизонтальной оси показано время в секундах, 


по вертикальной для ЭКГ — микровольты 


Во всех ритмограммах из этого набора есть промаркированные аномальные участки. На рис. 2 аномальные 


участки выделены красным пунктиром. По оси х показана продолжительность одного замера ритмограммы в 


миллисекундах, по у — интервал между соседними ударами сердца в миллисекундах. 
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Рис. 2. График ритмограммы № 69: красным пунктиром выделены аномальные участки, 
черными буллитами — аномальные точки 


Каждая ритмограмма представлена своим идентификатором. Продолжительность ритмограмм в 
исследуемом наборе данных различная: 4 минуты в среднем, 31 минута максимум. Каждый замер внутри одной 


ритмограммы имеет метку времени в миллисекундах от начала записи. Продолжительность КК-интервала 
также представлена в миллисекундах. Каждое конкретное значение в ритмограмме позволяет говорить о 
признаках аномалии (0 — нет, | — есть). 2,53 % наблюдений маркированы цифрой 1. Таким образом, набор 
данных имеет сильный дисбаланс классов, что типично для задач обнаружения аномалий. 

В разметке данных встречаются различные подходы к выделению аномальных участков. Как аномальные 


выделялись группы точек в окрестности характерного пика и падения продолжительности ритма сердца: 3-го, 
4-го, 6-го измерений (рис. 2). Не всегда количество точек в окрестности размечено одинаково — слева и справа 
от пика может быть разное количество аномальных точек. Кроме того, выявлены ритмограммы с 
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зашумленными показаниями. Так было при потере связи с гаджетом и замерах при установке или снятии 
прибора. 16 ритмограмм с некорректными данными пришлось исключить из рассмотрения, а разметку 
переделать: 

— у аномального участка выделяется лишь одна точка, характеризующая аномальный фрагмент; 

— аномальные точки обозначены черными буллитами (рис. 2). 

Извлечение признаков. В чистом виде представлен лишь один сигнал — значение интервалов между 
сердцебиениями. Поэтому для уточнения модели подготовили дополнительные признаки на основании 
имеющегося сигнала: отклонение от медианного значения и производные в соседних замерах ритма. Этот 
перечень признаков выбрали после разведочного анализа данных и визуального выявления паттерна в местах, 
соответствующих аномальным участкам. На рис. 2 их обозначили красной пунктирной линией. 

Результаты исследования 

Метрика для оценки качества детекции аномалий. Для оценки качества модели в задаче бинарной 
классификации ввиду дисбаланса классов задействовали метрику РЁ, [7] (1). Она позволяет оценить, насколько 
хорошо построенная модель детектирует редкий класс. В данном контексте под редким классом понимаются 


аномальные по продолжительности сердцебиения —= сердцебиения с аномальным ритмом: 
точность х полнота 
Е =2х 


ОЕ и 

Здесь: 

— точность — доля правильно детектированных моделью аномальных сердцебиений от общего количества 
сердцебиений, которые модель определила как аномальные; 

— полнота (или иными словами — чувствительность) — доля сердцебиений, которые модель верно 
детектировала как аномальные, от общего количества аномальных сердцебиений во всем наборе данных. 

Алгоритмы машинного обучения. В рамках исследования применены пять описанных ниже алгоритмов 
машинного обучения. 

1. зоайопРогезх — алгоритм с неконтролируемым самообучением на базе экстремально 
рандомизированных решающих деревьев [8]. 

2. Таз Отафеп Воозйпе Масбше Сазяйег (ГОВМ!азШНег) — алгоритм градиентного бустинга над 
решающими деревьями [9]. Для повышения скорости работы задействуются две техники: Ога еп(-Базе4 Опе- 
$14е ЗатрНие и Ехсшеуе Ееацие ВипаНп?°. 

3. ВапаотРогезСЛаз1Нег базируется на решающих деревьях и реализует многократный выбор случайного 
подмножества признаков. По ним строятся более простые оценщики —щ деревья решений. Результаты 
агрегируются для получения конечного предсказания [10]. 

4. ЕхцаТтеезСаз; ег аналогичен КапдотРогез аз Шег, однако в нем дополнительно реализован 
случайный выбор границы, по которой происходит ветвление узлов в деревьях решений [11]. 

5. $СРОпеС1аз$УМ“ — линейная версия Опе-СЙазз бирром Уесюг Масьше с использованием 
стохастического градиентного спуска. 

Тзо1вопРоге и $ОООпеСазз5УМ были выбраны ввиду их широкого использования в задачах детекции 
аномалий [12, 13]. ГОВМАазяШег, КапдотЕРогезСазНег и ЕхнаТтеезСазИег достаточно хорошо 
показывают себя в разных задачах, поэтому их тоже задействовали для сравнения результатов. 

Особенность алгоритмов 15091айопРогез и $СООпеСаз53УМ заключается в том, что они не требуют на 
входе четкой разметки аномальных наблюдений, в то время как для остальных использованных в исследовании 
алгоритмов она обязательна. 

[5о]анопРогезЕ базируется на предположении, что при построении изолирующих деревьев аномальные 
наблюдения можно изолировать (отделить) за меньшее количество операций, чем нормальные экземпляры 
наблюдений. Для каждого наблюдения алгоритм вычисляет значение оценки аномальности (апота]у зсоге) по 
формуле: 

_ Ето) 
5(х,п) =2 с®, (2) 


где А(х) — число ребер до экземпляра х в каждом изолирующем дереве решений; Ё (в (х)) — среднее значение А(х) 


на всем наборе изолирующих деревьев; с(п.) — нормализирующая константа для набора данных размером п (3). 


3 оМСВМ: А НзШу ЕЁйс<епЕ Стафеге Воозйпе Гесё\юп Тигее // ммм. писгозой.сот : [сайт]. УВТ: Вирз://Аммлу.пусгозой.сот/еп- 


и5/гезеагсЬ/\ир-сошепИир!юа4$/2017/1 1/12 еб т.раЕ (дата обращения: 10.09.2022). 
* Оше Опе-СМ!азз 5УМ /  $сщИ-еага Чеуе!орегз  (В$О  Шасепзе) //  зсЖи-еагп.оге: [сайт]. УВЫ: № рз://зсКИ- 
1еагп.ого/а ето и]ез/з А.В #юппе-опе-с1аз$-зут (дата обращения: 10.09.2022). 
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2(п-1) 


с(п) =2Нт- 1) - — (3) 
Н(®) = ш(®) +у. (4) 
В уравнении (4) у — постоянная Эйлера, равная 0,57721... 
Если наблюдение х имеет значение оценки аномальности $, близкое к 1, то оно считается аномальным. Если 
5 близко к 0,5, то наблюдение не имеет очевидных признаков аномальности. Если $ близко к 0, то наблюдение 
может считаться нормальным (рис. 3). 
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Рис. 3. Нормальные и аномальные наблюдения в изолирующих деревьях 


$ОРОпеС!аз$5УМ основывается на противоположном относительно [50]айопЕогезЕ подходе. Алгоритм 
определяет границы нормальных наблюдений и все новые наблюдения сопоставляет с границами этой нормы, 
чтобы выявить аномалию. 

Значимость признаков. Оценка степени влияния признаков на прогностическую способность модели 


представлена на рис. 4. 
ы ОО 
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Рис. 4. Диаграмма значимости признаков: х — текущее значение интервала; пех{_ 4 — производная в последующей точке 
измерения продолжительности сердцебиения; ргеу_АНР — производная в предыдущей точке измерения продолжительности 
сердцебиения; Ч! — первая производная в текущей точке; теФап_4е\аноп — отклонение текущего значения 
длительности АК-интервала от медианы 


Для расчета числовой оценки значимости использовался встроенный в ГОВМСЙазяШег механизм, 
возвращающий через свойство }{еайиге_ппроптапсеу_ обученной модели массив числовых оценок для каждого 
признака. Значимость в моделях на основе градиентного бустинга над решающими деревьями, как правило, 
рассчитывается на основе индекса Джини (Сии-ипригиу Ш4ех”) [14], используемого в процессе определения 
точек ветвления при обучении модели: 


ста) =1-УЕр?. (5) 


> КагаБфег Е. Сим ппригиу // 1еагпда!азс!.сот : [сайт]. ОВТ: Вирз://м\м\ 1еагп4а(азс1.сот/еоззагу/1-нпригИу/ (дата обращения: 10.09.2022). 
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Здесь Я — набор наблюдений, подходящих по условиям в рассматриваемой точке ветвления, 4 ЕР; К — 
количество классов, представленных во всем тренировочном наборе данных 0; р; — вероятность 
принадлежности наблюдений к классу 1 в рассматриваемой точке ветвления решающего дерева. 

Самыми значимыми оказались следующие признаки: текущее значение интервала (х), производная в 
последующей (пех! _Ч1!) и предыдущей (ргеу_А) точках измерения продолжительности сердцебиения (рис. 4). 
Полный перечень используемых признаков приведен в таблице 2: 


Таблица 2 

Перечень используемых признаков 

№ Признак Описание 

1 х КК-интервал в текущей точке измерения 

2 пех АН Первая производная в следующей точке 

3 ргеу_ЧЕ Первая производная в предыдущей точке 

4 ЧЕ Первая производная в текущей точке 

5 |шеФфап_аемайоп | Отклонение текущего значения длительности АК- 
интервала от медианы в рамках одной ритмограммы 


Сравнение моделей. Для результативности модели $СООпеС1аз55УМ важен подбор параметра пи, который 
соответствует вероятности обнаружения регулярного наблюдения за пределами границы нормы. Иными 
словами, пи определяет верхнюю границу доли ошибок при обучении модели и нижнюю границу доли опорных 
векторов. Для подбора пи с учетом природы имеющихся данных дополнительно оценивалась метрика качества 
при различных значениях указанного параметра (рис. 5). В итоге выбран пи, равный 0,11. 
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Рис. 5. Оценка параметра пи (по горизонтальной оси) для модели 5СДОпеС1а555УМ. 
На вертикальной оси — значения метрики РЁ! 


Для расчета метрики качества на различных моделях использовался весь набор данных с применением 
техники перекрестной проверки. Внутри одной ритмограммы мы имеем временной ряд наблюдений, снятых за 
один непрерывный промежуток времени у одного человека, поэтому следует рассматривать их как 
зависимые [15]. Для разделения данных на обучающие и тестовые наборы применялась следующая стратегия. 
Отобранный набор данных состоит из 213 ритмограмм, помеченных уникальным идентификатором (14). Это 
дает возможность выделить ритмограммы для обучения и тестирования моделей. Набор ритмограмм для теста 
можно случайным образом выбирать по идентификаторам. Ниже описан подход, примененный в 
представленной работе. 

Г. В цикле разделения данных выполняются пять действий. 

1. Фиксируется начальное число для генерации псевдослучайных чисел (зее4) — пр.гап4отп.зееа (а), где 
Гю14 — номер текущего разбиения данных. 

2. Генерируется 42 случайных целочисленных значения в диапазоне от 1 до 213. Так мы получаем 
случайные номера идентификаторов ритмограмм для тестового набора данных. 


$ $СООпеСаз5УМ — @4оситегаНоп. ЗсЖИ-еат  4еуеорез (ВЗР  П/сепзе) — зсЖИ-еагп.оге [сайт]. ПВГ: — Б@рз://зс П&И- 
1еагп.ого/зае/подиез/сепега{е/5Кеаго.Нпеаг_то4е1.5@ООпеС!аз;5УМ.В аи Юеага.Нпеаг_то4е1.5СООпеС1а55$УМ (дата обращения: 
10.09.2022). 
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3. Вносятся в отдельный список номера идентификаторов ритмограмм, которые остались после отбора 
идентификаторов для теста. Их задействуют для тренировочного набора. 

4. На ритмограммах из тренировочного набора обучаются модели, а на ритмограммах из тестового набора 
оцениваются метрики качества прогнозирования. 

5. Записывается значение метрики качества для каждой модели, посчитанной на тестовом наборе 
ритмограмм при текущем разбиении данных. 

П. Шаги 1-5 повторяются для каждого номера разбиения данных. 

Ш. Полученные значения метрики качества усредняются для каждой из моделей. 

Сравнительная оценка усредненной метрики качества прогнозирования для каждой модели приведена в 


таблице 3. 
Таблица 3 
Оценка метрики качества Ё\ 

Модель Метрика Е! * 
ГОВМС!азяШег 0,8328 
КапдотЕогеазяШЙег 0,7638 
ЕхцаТтеезСЛаз ег 0,7369 
$О0Опе а 5УМ 0,0169 
[зо]айопРоге$ < 1е-4 
*Среднее значение при выбранной стратегии 
перекрестной проверки на пяти разбиениях. 


Обсуждение и заключения. Разработана математическая модель обнаружения аномалий в ритме сердца с 
точностью 83 %. По метрике качества ЕЁ, лучшей оказалась модель на базе алгоритма ГОВМаз ег. 
Т5оайопЕоге$ и З@ООпеС!а555УМ на текущих данных показали слабые результаты. 

Предложенную модель можно реализовать в составе программной части носимых персональных смарт- 
устройств. Предлагаемый сценарий использования решения: 

— запись ритмограммы активируется на персональном носимом устройстве через пользовательский 
интерфейс; 

— по окончании запись подается в разработанную модель для анализа; 

— по итогам анализа данных математическая модель выдает уведомление о наличии или отсутствии 
аномалий на экране носимого устройства. 

Отметим, что для записи одной ритмограммы, видимо, достаточно в среднем 4 минут. За это время 
возможно обнаружение ковидных аномалий в ритме сердца. 

Модель занимает в памяти носимого устройства 493 килобайтов, что вполне подходит для практического 
применения. Решение опирается только на информацию о ритме сердца и не задействует факторы, недоступные 
для мобильных персональных гаджетов. 

Повышение точности детекции аномалий предполагает дополнительные изыскания. Их следует 
сфокусировать на разработке уникальных признаков, которые выявляются по исходному сигналу ритма сердца. 
Однако текущее решение уже дает возможность оперативной и легкой оценки вероятности СОУТ-19 на 
ранней стадии. Это наряду с выполнением рекомендаций медиков может дополнительно способствовать 
снижению риска смертности от негативного влияния коронавирусной инфекции на сердечно-сосудистую 
систему. 
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Отказоустойчивый кластер хранилища данных для 
аналитических запросов в банковской сфере 
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Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, 
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Аннотация 

Введение. Банковский сектор придает большое значение хранению данных, поскольку это критически важный 
аспект бизнес-операций. Объем данных в данной сфере неуклонно растет. С увеличением объемов данных, 
которые необходимо хранить, обрабатывать и анализировать, крайне важно выбрать подходящее решение для 
хранения данных и разработать необходимую архитектуру. Представленное исследование направлено на то, 
чтобы заполнить пробел в существующих знаниях СУБД, подходящих для банковского сектора, а также 
предложить способы для отказоустойчивого кластера хранения данных. Цель работы — анализ ключевых 
СУБД для аналитических запросов, определение приоритетов СУБД для банковского сектора и разработка 
отказоустойчивого кластера хранения данных. Для выполнения требований к производительности и 
масштабируемости предложено решение для хранения данных с отказоустойчивой архитектурой, отвечающее 
требованиям банковского сектора. 

Материалы и методы. Анализ предметной области позволил создать набор характеристик, которым должна 
соответствовать СУБД для аналитических запросов (ОГАР), выполнить сравнение некоторых популярных 
ОГАР СУБД и предложить отказоустойчивую кластерную конфигурацию, написанную на языке хи], 
поддерживаемую СУБД СПсКНоцзе. Автоматизация выполнена с помощью Апз1Ые Р1ауБооК$. Он интегрирован 
с системой управления версиями С\ИаБ и шаблонами Лца. Таким образом достигается быстрое развертывание 
конфигурации на всех нодах кластера. 

Результаты исследования. Для баз данных ОГАР были разработаны критерии, проведен сравнительный 
анализ нескольких популярных систем. В результате была предложена надежная кластерная конфигурация в 
банковской индустрии, которая удовлетворяет требованиям аналитических запросов. Для увеличения 
надежности и масштабируемости СУБД процесс развертывания был автоматизирован. Также приведены 
детальные схемы конфигурации кластера. 

Обсуждение и заключения. Составленные критерии для ОГАР СУБД позволяют определить необходимость 
данного решения в организации. Сравнение популярных СУБД может быть использовано организациями для 
минимизации затрат при выборе решения. Предлагаемая конфигурация кластера хранилища данных для 
аналитических запросов в банковской сфере позволит повысить надежность СУБД и удовлетворить требования 
к последующей масштабируемости. Автоматизация развертывания кластера путем механизма шаблонизации 
конфигурационных файлов в АпУчЫе Р1аубооК$ позволяет настроить готовый кластер на новых серверах за 
минуты. 


Ключевые слова: СУБД, ОГАР, хранилище данных, СПсКНоц$е, отказоустойчивый кластер. 
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Введение. Хранилище данных в банковской сфере является одним из ключевых бизнес-факторов. Для 
обеспечения безопасности информации о клиентах и транзакциях, необходимо провести меры защиты, 
распределения и создания резервных копий. Для оперативного анализа сотрудники должны иметь возможность 
делать оперативные аналитические запросы в хранилище данных, при этом не мешая работе других процессов 
внутри организации и не вызывая большую нагрузку на само хранилище. Базы данных (РаёаБазе) и хранилище 
данных (Раа У\/агерой$е) — это информационные системы, в которых производится хранение данных, но они 
используются и для решения различных задач. В статье описано, что делают такие системы, в чем основные 
различия между ними и почему их эффективное использование необходимо для развития бизнеса. 

Многие организации допускают ошибки в проектировании архитектуры баз и хранилищ данных, упуская из 
виду аспекты информационной безопасности, масштабируемости и отказоустойчивости. Актуальность этой 
проблемы обусловлена интенсивным развитием систем в банках, расширением сфер их применения и 
увеличением количества данных, нуждающихся в постоянном анализе. Для оперативного анализа большого 
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количества данных необходимо хранилище, которое должно соответствовать всем требованиям надежности и 
безопасности. 

Эффективные процессы принятия решений в бизнесе зависят от качественной информации. В современной 
конкурентной бизнес-среде требуется гибкий доступ к хранилищу данных, организованному таким образом, 
чтобы повысить производительность бизнеса, обеспечить быстрое, точное и актуальное понимание данных. 
Архитектура хранилища данных разработана для удовлетворения подобных требований и является основой 
этих процессов [1—5]. 

Цель работы — определение приоритетной СУБД для выполнения аналитических запросов в банковской 
сфере и проектирование отказоустойчивого кластера хранилища данных. Данное решение существенно 
повысит скорость выполнения аналитических запросов, решит проблемы с масштабируемостью и надежностью 
хранилища данных. 

Материалы и методы. База данных (Ра(аБазе) хранит информацию в режиме реального времени об одной 
конкретной части бизнеса, ее основная задача заключается в обработке ежедневных транзакций. Базы данных 
используют оперативную обработку транзакций (ОГТР) для быстрого удаления, вставки, замены и обновления 
большого количества коротких онлайн-транзакций. 

Хранилище данных (Ра{а \агенойзе) — это система, которая собирает данные из множества различных 
источников внутри организации для составления отчетов и анализа, используя оперативную аналитическую 
обработку (ОГ.АР) для быстрого анализа больших объемов данных. Данная система сосредоточена на чтении, а 
не на изменении исторических данных из множества различных источников, поэтому соблюдение требований 
АСШ (Аюпис, Сопз1\епь 150]ае4 ап Пига Ме) менее строгое. Хранилища данных выполняют сложные 
функции агрегирования, анализа и сравнения данных для поддержки принятия управленческих решений в 
компаниях. 

Хранилище в банковской сфере может содержать: 

— учетную информацию о клиентах (персональные данные, адреса, телефоны); 

— информацию о банковских продуктах и услугах (кредиты, депозиты, пластиковые карточки, мобильный 
банк ит. д.); 

— данные об операциях (включая карточные) в минимальной детализации за последние 3 года; 

— сведения о счетах, остатках на них ит.д. 

Для удовлетворения потребностей в онлайн-аналитической обработке запросов (ОГАР) существуют 
отдельные типы систем управления базами данных (СУБД) [3—6]. Каждая из систем имеет свои особенности в 
построении архитектуры. 

Для проведения эффективного анализа соответствия указанным требованиям хранилища должны: 

— обладать высокой вместимостью, способной вместить огромные объемы данных (миллиарды или 
триллионы строк); 

— быть организованы в виде широких таблиц с множеством столбцов; 

— выполнять запросы с небольшим количеством столбцов; 

— иметь высокую скорость выполнения запросов (в миллисекундах или секундах); 

— предусматривать большую часть запросов только на чтение; 

— поддерживать быструю массовую загрузку данных при их обновлении (более 1000 строк за раз) и 
добавлении, но без их изменения; 

— обладать высокой пропускной способностью для обработки одного запроса (до миллиардов строк); 

— обладать высокой надежностью; 

— обеспечивать безопасность и консистентность данных. 

Для ОГАР-сценария работы в банковской сфере предпочтительнее использовать колоночные 
аналитические СУБД, поскольку в них можно хранить много столбцов в таблице, что не будет сказываться 
на скорости чтения данных. Колоночные СУБД обеспечивают сильное сжатие данных в столбцах, так как в 
одной колонке таблицы данные, как правило, однотипные, чего не скажешь о строке. Они также позволяют 
на более маломощном оборудовании получить прирост скорости выполнения запросов в десятки раз. При 
этом, благодаря компрессии, данные будут занимать на диске в 5-10 раз меньше места, чем в случае с 
традиционными СУБД [7-11]. 

В ходе анализа требований выбраны следующие колоночные СУБД: СПсКНоцзе, Уегиса, Атахоп Вед. 

СИсКНоцзе является предпочтительным решением из-за следующих преимуществ: открытый исходный код; 
есть возможность определять некоторые или все структуры, которые будут храниться только в памяти; высокая 
скорость работы; хорошая степень сжатия данных; Юр и интерфейс командной строки; кластер можно 
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горизонтально масштабировать; высокая доступность; легкость установки и настройки. Установка 
осуществляется на серверах организации в изолированном сегменте, что отвечает требованиям безопасности 
для чувствительных данных в банковской сфере. Также СУБД включена в реестр отечественного программного 
обеспечения, поэтому позволяет внедрять данный программный продукт в государственные компании. 

Решение Ататоп Кедз ШИ предоставляется только в виде облачной службы. Для организаций из банковского 
сектора, которые не могут размещать свои данные в облаках по ряду причин, связанных с безопасностью, 
данный продукт теряет свою привлекательность. 

\Уегаса является альтернативным вариантом СПсКНоизе с платной лицензией для крупных кластеров и 
возможностью установки на локальные серверы компании. 

Ниже представлена реализация архитектуры распределенного хранилища данных. Для повышения 
отказоустойчивости и производительности предлагается реализация распределенного отказоустойчивого 
кластера СПсКНочзе с 3 шардами и 2 репликами. 

Шардинг (горизонтальное масштабирование) позволяет записывать и хранить части данных в 
распределенном кластере, обрабатывать и считывать их параллельно на всех нодах кластера, увеличивая 
пропускную способность данных. 

Репликация — копирование данных на несколько серверов, поэтому каждый бит данных можно найти на 
нескольких нодах. 

Масштабируемость определяется шардированием или сегментацией данных. Надежность хранилища 
данных определяется репликацией данных [12—16]. 

Шардирование и репликация полностью независимы, за них отвечают разные процессы. Необходимо 
локализовать небольшие наборы данных на одном шарде и обеспечить достаточно ровное распределение по 
разным шардам в кластере. Для этого в качестве ключа шардирования рекомендуется брать значение хэш- 
функции из поля в таблице. 

В зависимости от количества доступных ресурсов и серверов, предлагается реализовать эту конфигурацию 
на 3 или 6 нодах. Для производственной среды рекомендуется использовать кластер из 6 нод. Следует 
отметить, что репликация не зависит от механизмов шардирования и работает на уровне отдельных таблиц, а 
также, поскольку коэффициент репликации равен 2, то каждый шард представлен в 2 нодах [17-19]. Ниже 
описаны варианты конфигурации. 


Схема логической топологии выглядит следующим образом: 

З(Шард) х 2(Реплики) = Кластер СИсКВочзе из 6-ти нод. 

Вероятность безотказной работы системы с 2 репликами и 3 шардами на 6 нодах равна: 

в =[1-(-р)?. 

Вероятность безотказной работы — это объективная возможность того, что система без восстановлений 
проработает время # [7, 13]. 

Таким образом, таблица, содержащая 30 миллионов строк, будет распределена равномерно на 3 ноды 
кластера. Остальные 3 ноды будут хранить реплики данных. При выводе из строя одной из нод кластера, 
данные будут браться из другой доступной ноды, которая содержит ее реплику, тем самым достигается 
надежность [20]. Кластер из 6 нод изображен на рис. 1. 


поде 01 по4е 02 поде 03 


зВага 01 зВага 02 зВага 02 
терНса 01 терНса 01 терНса 02 


по4е 04 по4е 05 по4е 06 


звага 03 зВага 03 зВага 01 
герИса 01 герИса 02 герПса 02 


Рис. 1. Отказоустойчивый кластер из 6 нод 
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Для репликации данных и выполнения распределенных ООГ, запросов необходимо использовать +1 ноду с 
установленным ГооКеерег. Можно использовать также СПскНоизе Кеерег, совместимый с ГооКеерег, не 
требующий установки на отдельном сервере. 

Пример фрагмента конфигурационного файла представлен на рис.2, из которого видно, что шарде 
настроена репликация для 1-ой и 6-ой ноды. 

<уапдех> 


<гетоте_зегуег$> 
<сТи$фег_1> 


<5пага> 
<ме198т>1</ме19Н+> 
<1пфегпа\_гер11са{1оп>{гие</1п%егпа\_гер11са{1оп> 
<гер11са> 


<|1о$1>4{ поде1 }}</во$+> 
<рог{>9000</рог+> 
</гер\1са> 
<гер11са> 
<по$1>{{ подеб }}</по$+> 
<рог{>9000</рог+> 
</гер11са> 
</зПага> 


Рис. 2. Фрагмент конфигурационного файла для 6 нод 


Вариант конфигурации кластера из трех нод с цикличной репликацией изображен на рис. 3. 
поде 01 по4е 02 поде 03 


зВага 03 
терНса 01 


звага 01 


терНса 01 Е 


терНса 01 


звага 02 


зрага 03 звага 01 


терНса 02 терНса 02 терНса 02 


Рис. 3. Отказоустойчивый кластер из 3 нод 


Для такой реализации требуется два разных сегмента, расположенных на каждой ноде. Основная проблема 
возникает из-за того, что каждый шард имеет одинаковое имя таблицы, СИсКкНоцзе не может отличить один 
шард/реплику от другого, когда они расположены на одном сервере. 

Для решения данной проблемы необходимо: 

— поместить каждый шард в отдельную базу данных (схему); 

— установить деРаш_АабаБазе для каждого шарда; 

— установить параметр и\цегпа! герПсайоп каждого шарда в ие; 

— использовать пустой параметр базы данных в сценарии РОГ, распределенной таблицы. 

Для такой топологии в промышленной среде требуется 6 серверных нод, где каждый сервер хранит данные 
только одного сегмента, обходной путь для отдельной базы данных не требуется. Для экономии ресурсов в зоне 
разработки или тестирования можно использовать конфигурацию с 3 нодами. 

Автоматизация выполнена с помощью Ап е Р!аубооК$ и интегрирована с системой управления версиями 
СлИаБ. Таким образом достигается быстрое развертывание конфигурации на всех нодах кластера. При 
изменении конфигурации ее можно применить на всех нодах с помощью одной команды или развернуть новый 
кластер СУБД за несколько минут [21]. 


В.В. Сивов и др. Отказоустойчивый кластер хранилища данных для аналитических запросов в банковской сфере 


Результаты исследования. Отказоустойчивый кластер аналитической СУБД обеспечивает резервирование 
для важных компонентов системы, что позволяет продолжительно функционировать даже в случае 
возникновения ошибок в отдельных узлах кластера. Это делается за счет распределения нагрузки, репликации 
данных между узлами кластера и высокой надежности компонентов, используемых в кластере. Результатом 
является увеличение доступности и надежности аналитической СУБД, что важно для бизнеса, в котором 
аналитические запросы играют ключевую роль. Отказоустойчивая кластерная конфигурация хранилища 
данных для аналитических запросов в банковской сфере с учетом автоматизации процесса разворачивания 
позволяет повысить надежность аналитического хранилища данных и удовлетворить требования к 
масштабируемости. Разработанная задача по автоматизации развертывания кластера с использованием 
механизма шаблонизации конфигурационных файлов в Апзе Р!аурооКз$ позволяет настроить готовый кластер 
на новых серверах за несколько минут. В задачи шаблона входят операции по установке необходимых пакетов, 
созданию необходимой конфигурации и запуску кластера. 

Пример конфигурационных файлов для автоматического развертывания кластера СУБД приведен на рис. 4. 
Расширение ]2 говорит о том, что они созданы с помощью механизма шаблонов Липа. Специальные 
заполнители в шаблоне позволяют писать код, аналогичный синтаксису Руфоп. В шаблон передаются 
параметры для автоматической вставки в финальный документ, тем самым достигается автоматическая сборка 
в зоны разработки, тестирования и промышленной эксплуатации, которая не требует ручного изменения 
файлов конфигурации. 

у 1етр!а{е$ 

сйскпоизе_сопЯд.хт|.]2 
сИсквоцзе_Кеерег.хгт|.]2 
сИскпоц$е_1Чар_ацщВ.хпп!.]2 
сИсквоцзе_1Чар_изег_атес{огу.хт!.]2 
сИсквоц$е_тасго_п1.хт|.]2 
сИскпои$е_тасго_п2.хп|.]2 
сИскпоц$е_тасго_п3.хт|.]2 


сИскпои$е_тасго_п4.хг|.]2 


сИскпои$е_тасго_п5.хи|.]2 


сИсквоцзе_тасго_пб.хпп.]2 


сИскпоц$е_изег$.хпп|.]2 


си${ег.хт|.]2 


Рис. 4. Конфигурационные файлы 


Описание конфигурационных файлов: 

сПсКвойзе_сопНз.хп1.]2 — общая конфигурация кластера; 

сПсКвоизе_Кеерег.хп11.]2 — конфигурация 7ооКеерег, которая отвечает за синхронизацию нод и репликацию; 

сИскпопзе_1Чар_аи.хп1.]2 — конфигурация подключения к ГРАР для обеспечения безопасности данных; 

сПсКвопзе_1Чар_изег_(тесогу.хп11.]2 — конфигурация с ролевой моделью по группам доступа для 
обеспечения безопасности данных; 

сИскпоцзе_тасго_п1(6).хп1.]2 — файлы макросов (для каждой ноды свой); 

сИскпоцзе_изегз.хп.]2 — конфигурационный файл создания локальных пользователей, необходимых для 
администрирования; 

сшазбег.х111.]2 — конфигурационный файл кластера. 

Для проверки надежности данной конфигурации проведен эксперимент, в ходе которого были загружены 
данные в кластер СУБД с коэффициентом репликации, равным 2. Созданы схемы 4\ и таблицы 
сТаз(ег_(езЕ Чайа на каждой из нод кластера СУБД, а также создана распределенная таблица на кластере 
амВ.саз{ег_1е5_ Чаба_ 41зи1Ьщед. Строки таблицы 4\0.4е${_ 4аа_ 1516 ще4, распределенной по кластеру, 
равны 27 547 855. Ниже перечислены строки таблицы 4\В.стаз(ег {ез(_айа с каждой из нод кластера: 

9186544 строк — 1-ая нода; 

9182959 строк — 2-ая нода; 

9182959 строк — 3-я нода; 

9178352 строк — 4-ая нода; 

9178352 строк — 5-ая нода; 

9186544 строк — 6-ая нода. 

Как можно заметить, таблица распределилась на весь кластер. Согласно конфигурации, приведенной на 
рис. 1, коэффициент репликации равен 2, значит, каждый блок данных будет представлен на 2 нодах. Это 
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можно увидеть из количества строк на нодах: шестая нода хранит копию первой, третья — копию второй, 
пятая — копию четвертой. 

Отказоустойчивость данной конфигурации можно проверить попеременным выводом из строя нод в 
кластере. Для этого можно выключить ноду или остановить сервисы на одной из нод командой зузетсй юр 
сПсКвоцзе-зегуег. В ходе эксперимента была выполнена остановка сервисов СУБД на нодах кластера. 

При одновременном отключении 3, 4, 6-ой или 1, 2, 5-ой нод, которые содержат реплики, пользователи 
продолжали получать данные из таблицы А\В.саз{ег_1ез_ аа_ Ч1зи1Ьщеа, и количество строк было 
равным 27 547 855. При выводе из строя одной из нод, данные продолжали отображаться, а количество строк 
было равным 27 547 855. При одновременном отключении нод, которые содержат реплику и оригинальные 
данные, происходила потеря данных. Данную конфигурацию можно масштабировать на 12 нод, тогда 
коэффициент репликации будет равен 3, а коэффициент шардирования — 6. 

Обсуждение и заключения. Предлагаемое решение может повысить скорость выполнения аналитических 
запросов, решить проблемы с масштабируемостью и надежностью хранилища данных в организациях 
банковской сферы. Авторами выполнена автоматизация развертывания кластера путем шаблонов в Ап Ые 
Р1аурооКз$, которая позволяет настроить готовый кластер на новых серверах за минуты. Данную конфигурацию 
можно масштабировать, увеличив количество нод и добавив их в конфигурационные файлы. 

Обозначен набор характеристик, которым должна соответствовать ОГАР СУБД, выполнено сравнение 
СУБД, предложена отказоустойчивая кластерная конфигурация хранилища данных для аналитических запросов 
в банковской сфере, выполнена автоматизация процесса разворачивания конфигурации. Подобное решение 
применимо для развертывания на ЕгееВЗО, Глпах, тасОЗ. Приведены схемы конфигурации кластера. Данная 
конфигурация решит проблему, связанную с надежностью и масштабируемостью, которая часто встречается в 
организациях. 
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Аннотация 

Введение. В настоящее время концепции туманных и краевых вычислений используются широким кругом 
приложений самой различной направленности. Одной из ключевых проблем организации вычислений в 
группах мобильных устройств, составляющих краевой/туманный слой, является обеспечение выполнения 
миссии на основе наличия заряда батареи. В связи с этим к настоящему времени представлено немало 
разработок, направленных на энергосбережение систем устройств. Однако очень важный аспект остается за 
рамками рассмотрения проблемы ресурсосбережения, а именно — вопрос сбережения остаточного ресурса 
вычислительного устройства. Целью данного исследования является формализация задачи распределения 
нагрузки как двухкритериальной задачи оптимизации и выбор базового метода ее решения. 

Материалы и методы. В рамках данной статьи предлагается подход к ресурсосбережению на основе 
оценивания двух критериев устройств: ресурса батареи и остаточного ресурса вычислительного устройства. 
Остаточный ресурс вычислительного устройства может быть оценен при помощи значений вероятности 
безотказной работы устройства или как величина, обратная интенсивности отказов с учетом того, что при 
моделировании используется экспоненциальный закон распределения отказов. На основе этого 
сформулирована модель задачи двухкритериальной оптимизации с учетом динамики топологии сети в процессе 
выполнения пользовательской миссии. Динамика топологии отражена в модели как последовательность 
топологий, каждая из которых соответствует определенному отрезку времени функционирования системы. 
Результаты исследования. На основании представленной модели задачи двухкритериальной оптимизации 
предложен метод обеспечения ресурсосбережения в краевом и туманном слое сети, отражающий специфику 
динамических слоев сети, а также учитывающий важность критериев оценивания расхода ресурсов устройств. 
Проведен эксперимент, позволяющий оценить влияние способа распределения задач по сетевому кластеру на 
вероятность безотказной работы устройств и на средний остаточный ресурс. 

Обсуждение и заключения. Проведенный эксперимент демонстрирует целесообразность применения 
разработанного метода, поскольку распределение задач по исполняющим устройствам оказывает существенное 
влияние (до 25 % по итогам эксперимента) на средний остаточный ресурс вычислительного устройства. 


Ключевые слова: ресурсосбережение, планирование вычислений, туманные вычисления, краевые вычисления, 
оптимизация. 
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Введение. В настоящее время широкое применение нашли приложения, использующие краевой и туманный 
сегменты сети [1—5]. Распределенные вычисления, выполняемые на узлах этих сегментов, имеют существенное 
отличие от распределенных вычислений, выполняемых в облачных структурах, а именно — относительно 
высокую динамику топологии и, помимо этого, наличие определенных расстояний между используемыми 
узлами. Это приводит к необходимости учета данного параметра при моделировании времени 
информационных обменов между исполняемыми задачами. Существование динамики топологии объясняется 
тем, что устройства края сети — это, как правило, пользовательские устройства (включая разного рода датчики, 
смартфоны, ноутбуки и т.д.), которые принадлежат определенным людям и могут быть перемещены, 
выключены, доступ к данным/ресурсам может быть запрещен. Все они при этом находятся в зависимости от 
наличия заряда батареи (источника питания) [6, 7]. Туманный слой сети менее динамичен, однако и для него 
характерно все перечисленное в той или иной степени: наличие определенного расстояния между узлами с 
присутствием транзитных участков сети, возможная мобильность самих узлов. Например, в случае обработки 
данных группами мобильных устройств. Интересным примером в этом смысле являются группы 
низкоорбитальных спутников, которые к настоящему времени используются для обработки данных, но при 
этом узлы смещаются относительно поверхности Земли с достаточно высокой скоростью (время доступа — 
порядка нескольких минут) и, следовательно, возникает вопрос о маршрутизации и передаче результатов 
обработки на наземные станции [8]. 
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В рамках вопросов ресурсосбережения используемых устройств в динамичных слоях сети был опубликован 
ряд работ [9-11], в которых внимание акцентируется на энергосбережении. Как правило, модель предполагает 
наличие зависимости потребляемой энергии устройства от его загруженности [12]. При этом вне фокуса 
исследований остался такой важный фактор как сбережение остаточного вычислительного ресурса устройства [13]. 

Остаточный вычислительный ресурс устройства — это величина, тесно связанная с такими надежностными 
характеристиками, как вероятность безотказной работы и гамма-процентная наработка на отказ. Также средний 
остаточный ресурс, с учетом использования экспоненциального закона распределения интенсивности отказов, 
характеризуется следующим выражением: 


1 [°.®) 
В(®) = =. ] Р(Хах, (1) 


где Р(>х) — вероятность безотказной работы (ВБР) объекта в течение времени х. 

Для экспоненциального закона распределения ресурса ВБР определяется как: 

Ро) =е*, (2) 
где /\, — интенсивность отказов, /,>0. 

Тогда согласно формуле (1), средний остаточный ресурс определяется по следующей формуле: 

вФ=В=А. (3) 
А, 

При этом имеются работы, посвященные связи величин ВБР, гамма-процентной наработки на отказ и 
среднего остаточного ресурса от температуры вычислительного элемента устройства, которая, в свою очередь, 
описана как функция от загруженности устройства [14—16]. Задача ресурсосбережения устройства также может 
быть поставлена как задача максимизации среднего остаточного ресурса устройства. 

Таким образом, проблема ресурсосбережения рассматривается в современных работах либо как задача 
минимизации расхода электроэнергии устройств, либо как задача сбережения остаточного вычислительного 
ресурса. 

Однако в современных условиях целесообразно говорить о двухкритериальной оптимизации расхода 
ресурсов устройств по следующим причинам: 

— электроэнергия — восполняемый, но критичный для выполнения миссии устройства ресурс, будь то 
датчик информационно-управляющей системы или управляемое мобильное устройство; 

— остаточный вычислительный ресурс определяет длительность целесообразной эксплуатации устройства, а 
время эксплуатации, соответственно, складывается из миссий, выполненных устройством. 

Поэтому целью данного исследования являлась разработка метода ресурсосбережения для устройств 
динамичных слоев сети, в рамках которого бы оценивание проводилось на основании двух критериев — 
расхода электроэнергии и расхода вычислительного ресурса устройства. 

Материалы и методы. Определим понятие «миссия группы устройств» как заранее определенную 
последовательность решения группой устройств комплекса вычислительных задач, связанных 
информационными обменами. Миссия группы устройств может быть описана ациклическим графом, вершины 
которого будут взвешены трудоемкостями задач, а дуги, соответственно, будут определять ограничения 
следования и объемы данных, передаваемых задачами. Пусть миссия описывается графом И = {м%, 3» Г}, где 
и; — трудоемкость задачи, <; — доля задачи, выполненная к моменту смены топологий, / — матрица объемов 
данных, передаваемых между задачами. Группа устройств в динамическом слое сети описывается набором 
графов, каждый из которых определяет топологию сети в момент времени Е. Будем считать, что время 
выполнения миссии включает некоторый временной интервал с дискретными моментами времени [15 1, ...1%1, 
каждому из которых соответствует граф топологии С‚={<ру,е,, Ку>, Н}, 1=1...т, ру — производительность ]-го 
узла 1-й топологии, е, — энергоресурс ]-го узла 1-й топологии, К, — средний остаточный вычислительный 
ресурс ]-го узла 1-й топологии, и1 — число устройств в группе, Н — матрица сетевых соединений между узлами. 
Также примем допущение о том, что в каждый последующий момент времени #„., для графов Снги С; 
справедливо следующее: по меньшей мере одна из вершин, принадлежащая С’, будет содержаться в С.+, что 
необходимо для продолжения выполнения миссии. 

Также будем считать, что обеспечивается, в случае изменения графа топологии в момент времени 1.1, 
возможность переназначения задач, которые в 1-Й топологии выполнялись на узлах, переставших существовать 
в {+/ момент времени, на любые из узлов топологии 1+/, имеющиеся в наличии. При этом сделаем допущение 
о том, что перенесенные на новые узлы задачи начинают выполняться на том же месте, где было прервано их 
исполнение. 
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Оценивание остаточного вычислительного ресурса будет реализовано следующим образом: для оценки ВБР 
устройства туманного слоя будем использовать выражения, предложенные в [17, 18]: 


АР 
аа, (4) 
где р— загрузка вычислителя (в долях). 
У’. 
ВЕ, (5) 
Р И : 1 опотани 


где и, — трудоемкость выполняемой работы. 


Соответственно может быть оценена ВБР (0 а также средний остаточный ресурс: 


0-е 2 РИеопзтаии 


сб ША. -1.2*2Л0 
Р=е“ =е"®"" =е | (6) 
ГДе Е.опзнайи — ограничение по времени, за которое должна быть выполнена задача узлом, чтобы уложиться в 
ограничение по времени выполнения миссии т 
При оценивании остаточного энергетического ресурса также будем полагать, что расход энергии 
пропорционален загруженности вычислительного узла и больше энергетических затрат не имеется: 


Ем = Ех — ёр , (7) 


где Ею —щ исходный уровень энергии в наличии; Е остаточный уровень; & — коэффициент, 


Лт 
выражающий зависимость между нагрузкой вычислительного элемента и его энергозатратами. 

Результатом работы двухкритериального метода обеспечения ресурсосбережения будет последовательность 
распределений нерешенных еще задач миссии за время [4, &1| по топологии С; с учетом того, что часть задач, 
находящихся в процессе решения на узлах, наследуемых из топологии С’;/, остаются на своих местах и не 
переносятся: А={А;}; 


А =. № = (8) 


Матрица А описывает распределение нерешенных задач, для которых 5,<1. При этом количество узлов 
соответствует новой топологии С’;. Узлы, унаследованные от топологии С’, перенумеровываются и занимают 
индексы, начиная с 1. 

Поскольку выполнение миссии может включать более одной топологии, проведем декомпозицию задачи 
следующим образом: воспользуемся «жадной» стратегией и будем полагать, что удовлетворяющее решение 
(минимизация ресурсопотребления, оцениваемое двумя выбранными критериями) может быть получено, если 
для каждой топологии мы будем выбирать наилучшее решение. 

То есть модель задачи будет иметь следующий вид: \8,,5 <Ь1:>0 получить такое А» что В(&,) 


КЕ, ) —> шах,е (1 


„)— тах ‚ при тах(А„ ) <1,. 

Рассмотрим возможную приоритетность сформулированных целевых функций. Для каждой отдельно взятой 
топологии необходимо минимизировать как расход электроэнергии, так и расход вычислительного ресурса. В 
некоторых частных случаях может быть достигнуто совпадение целей оптимизации для двух критериев. 
Например, задача распределяется на узел таким образом, что обеспечит ему минимальную загруженность 
вычислительного элемента и при этом минимальный расход энергопотребления. 

Однако, такая ситуация возможна не всегда. Например, имеется узел, загруженный менее прочих 
вычислениями, но при этом его энергоресурс практически исчерпан. Назначая на него задачу, мы ставим 
выполнение миссии под угрозу (либо необходима смена топологий по мере окончания запаса электроэнергии), 
а назначая задачу на узел более загруженный, мы ухудшаем значение критерия остаточного вычислительного 
ресурса. Противоположная ситуация также может иметь место: узел с большим остаточным энергоресурсом 
может иметь низкий остаточный вычислительный ресурс как результат участия в предыдущих миссиях 
устройств. 

Поэтому целесообразно провести ранжирование целевой функции (ЦФ) по важности, на основе чего 
выбрать дальнейший метод решения. 

Для миссий, критически важных к выполнению, ведущую роль играет ресурс батареи устройства [19—21]. 
Соответственно, ЦФ, определяющая остаточный вычислительный ресурс, может быть сведена к 
дополнительному ограничению. Невыполненные задачи миссии будут дораспределяться в зависимости от 
остаточного заряда батареи. 
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Следует отметить, что и для миссий, не обладающих критичностью выполнения, справедливо сказанное 
выше. В случае, если устройство лишится заряда батареи, не выполнив миссию, оно может быть утеряно. Тогда 
оптимизация вычислительного ресурса станет бессмысленной. 

С другой стороны, сводя остаточный вычислительный ресурс к ограничению, решение может быть не 
получено вовсе по причине отсутствия доступных узлов, характеристики которых бы ему удовлетворяли. 

Результаты исследования. Многокритериальные задачи оптимизации, как правило, сводят к решению 
одной и более однокритериальных задач. В данном случае более целесообразно оставить одну ЦФ, а вторую 
представить в качестве ограничения. 

Для рассматриваемой системы выполнение миссии критично для каждого узла, поэтому преобразуем ЦФ по 
энергопотреблению в ограничение, выполнение которого отсечет неприемлемые варианты по распределению 
задач по узлам. Далее в рамках полученного множества узлов распределение будет происходить с 
минимизацией расхода остаточного вычислительного ресурса. 

Опишем основные этапы двухкритериального метода обеспечения ресурсосбережения, как основного 
результата проведенного исследования: 

— на имеющейся С;-топологии выбрать узлы, обладающие достаточным для завершения миссии без смены 
топологии энергоресурсом; 

— произвести распределение задач миссии по выбранным узлам таким образом, чтобы для незанятых узлов, 
к которым будут прикреплены нерешенные еще задачи, минимизировать расход вычислительного ресурса. 

Допустимо ранжирование узлов, полученных на первом шаге метода, по убыванию наличия энергоресурса. 
Но тогда, распределяя задачи в первую очередь по узлам с максимальным количеством энергоресурса, может 
возникать ухудшение значений критерия остаточного вычислительного ресурса. 

Проиллюстрируем зависимость значений ВБР устройств от выбора узла для решения задачи (в данном 
эксперименте сравниваем состояния узлов при передаче поступающих данных и при их обработке). 

Топология представлена на рис. 1. 
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Рис. 1. Экспериментальный фрагмент топологии сети 


Эксперимент проводится для 10 временных отрезков, каждый из которых по 100 часов при следующих 
значениях параметров: \\тесеуе (объем получаемых задачей данных) — 500 ед., \зепа (объем отправляемых 
задачей данных) — 100 ед., \/ргосез$ (трудоемкость обработки данных) — 150 ед., Т4ес! (время, отводимое на 
решение задачи) — 50 ед. 

Далее приведены графики ВБР-узлов в зависимости от того, обрабатывает ли узел данные или только 
передает их (рис. 2—5). 
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Рис. 2. Сравнение значений ВБР для узла 1 при обработке данных (Рр1) и при передаче данных (Ри) 
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Рис. 3. Сравнение значений ВБР для узла 2 при обработке данных (Рр2) и при передаче данных (Ри?) 
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Рис. 4. Сравнение значений ВБР для узла 3 при обработке данных (Рр3) и при передаче данных (РИ) 
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Рис. 5. Сравнение значений ВБР для узла 4 при обработке данных (Рр4) и при передаче данных (Ри4) 


Далее на рис. 6 представлено сравнение различий в значениях среднего остаточного вычислительного 
ресурса. 
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Рис. 6. Сравнение значений среднего остаточного ресурса для узлов 1-4 при транзите данных (Кетат_6гапзй) и при 


обработке данных (Кетат_регогт) 


Обсуждение и заключения. В статье предложен  двухкритериальный метод обеспечения 
ресурсосбережения в краевом и туманном слоях сети, который отличается от ранее представленных в 
предметной области организации распределенных вычислений в динамических слоях сети тем, что учитывает 
одновременно состояние батареи устройства (что критично для выполнения миссии) и остаточный ресурс 
вычислительного устройства (что влияет на длительность целесообразного использования устройства за 
пределами выполняемой миссии). 

В основу разработанного метода положена сформулированная задача двухкритериальной оптимизации, 
которую затем предлагается свести к однокритериальной, переведя ЦФ остатка энергоресурса в ограничение. 

Отсекая узлы, заранее не подходящие по уровню заряда батареи, далее распределение задач происходит по 
критерию остаточного вычислительного ресурса. 

Проведенный эксперимент позволяет сравнить результаты распределения задачи по узлам при наличии 
информационных взаимодействий между задачами. Наглядно показано, что способ распределения задач 
существенно влияет на величину остаточного ресурса вычислительного устройства (до 25 %). Таким образом, 
возможны улучшения показателей ресурсопотребления устройств, участвующих в вычислениях в краевом и 
туманном слоях сети. 
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Аннотация 

Введение. Существенная нелинейность химических реакций кислот и оснований определяет алгоритмы 
управления в режиме подкисления или подщелачивания, то есть периодической дозации минимального объема 
кислоты или щелочи. Такое регулирование может быть малоэффективным, т. е. допускает недостаточную или 
избыточную концентрацию контролируемого вещества. Статья посвящена проблеме точного регулирования 
водородного показателя в мини-биореакторах. Предлагается использовать цифровую модель системы 
управления кислотностью для подбора концентраций доливаемых растворов, определения методики 
регулирования и повышения точности. Цель работы — обеспечение требуемой точности регулирования рН 
в ш уйто мини-модели желудочно-кишечного тракта статического типа. 

Материалы и методы. Исходная структурная схема модели включает накопители и потоки. Она представляет 
собой базу для основных дифференциальных уравнений, характеризующих изменение объема и кислотности. 
Для корректировки показаний кислотности результирующего раствора по температуре создана статическая 
модель, основанная на полиномиальной аппроксимации экспериментальных данных методом наименьших 
квадратов. В прикладном пакете МаЧаб ЗипиПпК исследованы структурные элементы математической модели. 
На реальной системе ш уйго моделирования искусственного желудочно-кишечного тракта домашней птицы 
определены переходные характеристики для подтверждения адекватности математической модели. 
Результаты исследования. В рамках данной работы авторы создали и проанализировали нелинейную 
математическую модель изменения рН в биореакторе с учетом внешних управляющих воздействий. Потоки 
раствора соляной кислоты, раствора щелочи и слива раствора из реактора представлены как элементы 
дифференциального уравнения, описывающего накопление жидкости в реакторе. Для повышения точности 
решение доработали с учетом температурной зависимости водородного показателя. Предложена 
математическая модель дозирования на основе регулятора с каналами щелочи и кислоты. Полученные данные 
позволили генерировать объединенную модель процесса регулирования рН в биореакторе. Адекватность 
решения подтвердили опытным путем. В виде структурных схем показаны модели: регулятора РН, 
регулирования объема содержимого в реакторе и химических реакций. Сравниваются переходные процессы 
математической модели и реальной системы управления. Установлено, что переходные характеристики 
математической модели и реальной системы идентичны по времени регулирования. Относительная 
погрешность регулирования реальной системы составила 0,35 %, а математической модели — 0,1%, что 


соответствует требуемой точности регулирования -0,1 рН. В виде графиков показано влияние исследованных 
потоков на реакцию нейтрализации. 

Обсуждение и заключения. Предлагаемая математическая модель позволит подобрать оптимальные методы и 
алгоритмы регулирования кислотности, что ускорит создание регулятора нелинейного процесса регулирования 
водородного показателя. В будущем эти наработки можно интегрировать в комплексную цифровую модель 
всего искусственного желудочно-кишечного тракта домашней птицы для оптимизации алгоритмов управления 
(дозации, перемешивания, периодичности и т. п.), а также аппроксимации к объектам ш У1у0. 
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Введение. С начала 1980-х годов исследователи физиологии и функциональных особенностей внутренних 
органов человека пришли к пониманию, что фокус экспериментов будет смещаться от практик ш У1\о к Ш 
уго [1]. Актуальность разработки систем аппаратной имитации физиологических процессов обусловлена: 

— снижением затрат и времени на проведение исследований: 

— возможностью единовременного проведения однотипных экспериментов; 

— акцентом на этичность исследований, что предполагает ограничение экспериментов на живых существах. 

К началу ХХПГвека сформировалось научное направление по разработке методик и оборудования для 
моделирования и исследований т уго [1]. 

В настоящее время известны статические и динамические ш уго модели желудочно-кишечного 
тракта (ЖКТ). В первом случае речь идет о моделировании процессов пищеварения при использовании жидких 
питательных сред [2]. Динамические модели, в отличие от статических, позволяют моделировать 
перистальтические движения в ЖКТ и используются при работе с пищевыми матрицами в виде суспензий, 
включающих твердые частицы и волокна. 

Цель данной работы — обеспечение требуемой точности регулирования рН в ш \о мини-модели ЖКТ 
статического типа. 

Для достижения цели решаются следующие задачи: 

— разработка нелинейной математической модели изменения рН в биореакторе с учетом внешних 
управляющих воздействий; 

— исследование математической модели методами имитационного моделирования; 

— доработка математической модели для повышения ее точности путем учета температурной зависимости 
водородного показателя; 

— создание математической модели процесса дозирования на основе двухканального конкурентного 
регулятора (по каналам щелочи и кислоты); 

— разработка объединенной математической модели процесса регулирования рН в биореакторе; 

— проведение экспериментов для подтверждения адекватности объединенной математической модели и 
определения достигнутой точности регулирования. 

Материалы и методы. С точки зрения моделируемого объема модели ЖКТ подразделяются на микро- (до 
50 мл), мини- (от50 до 400 мл) и макросистемы (свыше 400 мл) [3]. Отдельно следует упомянуть 
микрофлюидные системы [4], но они выходят за рамки данной работы. 

Патентный анализ позволяет назвать первое зарегистрированное изобретение в данной сфере — «т уйго 
модель ш У\о пищеварительного тракта»'. Эта динамическая система предшествовала появлению систем типа 
ТИМ? [5]. Идея развивалась в моделях ООМ [6], ТИМ? [5], НО$ [7] и др. Статические модели (ОТООТ [8], 
У1МОТ [9], 5НИМЕ [10], АВСОЕ. [11] и др.) ушли от сложного моделирования перистальтики и сфокусировались 
на более точном моделировании собственно процессов пищеварения. Первая причина — техническая 
сложность и дороговизна реализации систем с динамическим сжатием содержимого. Вторая — развитие 
методов математического моделирования и микроконтроллерных систем управления, позволивших мониторить 
среду в пробирке (реакторе) с высокой точностью и стабильностью [4]. 

В данной работе рассматривается решение частной задачи математического моделирования процесса 
регулирования рН для ш УЙго мини-модели желудочно-кишечного тракта (ЖКТ) статического типа [12]. 
Точность управления кислотностью химуса определяется жесткими требованиями к условиям экспериментов п 
уго. Точность автоматического регулирования рН в биореакторе не может быть ниже +0,1 единицы.” Это 
принципиально для правдоподобия условий процессов пищеварения и соблюдения требований к 
повторяемости результатов опытов. При этом рН может быть от слабощелочной (7,9 рН, ротовая полость) до 
сильно кислотной (1,3 рН, желудок). 


Для моделирования реакции нейтрализации воспользуемся методом накопителей и потоков [13]. Его база — 
это законы сохранения (накопители), а также математическое описание закономерностей, изменяющих 
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количество накапливаемого признака (потоки). Структурная схема данной модели с двумя накопителями 
представлена на рис. 1. 


(1) Объем 
кислоты 5 
1 
Емкость (3) Сброс 
(3) Объем 
щелочи 
Реактор 
(4) Поток 
Маон о, (6) Процессы 
Раствор 
(5) Поток пищеварения 
НС1 


Рис. 1. Структурная схема модели с двумя накопителями 


Реактор как физический объект накапливает некоторое количество: 
— жидкости (01); 
— ионов водорода (О). 
В этой емкости идет контролируемый процесс пищеварения с изменением рН и объема жидкости. 
На первый накопитель действуют три потока: 
— 1-й — объем доливаемой кислоты; 
— 2-й — объем доливаемой щелочи; 
— 3-й — объем, выливаемый из реактора. 
На второй накопитель также влияют три потока: 
— 4-й — подача насосом гидроксида натрия [2]; 
— 5-й — поступление соляной кислоты [2]; 
— 6-й — изменение рН в процессе пищеварения. 
Последний поток представляет собой возмущающее воздействие и не исследуется в данной работе. 
Результаты исследования. Рассмотрим математическую реализацию модели для описания физико- 
химических процессов. 
Объем жидкости в реакторе: 
У=5. 1. (1) 
Изменение объема жидкости в реакторе: 
аУ =5 -ай. (2) 
Для удобства назовем: 
— поток раствора соляной кислоты — 41 = С1; 
— поток раствора щелочи (гидроксид натрия) — 42 = Со; 
— поток слива раствора из реактора — 42 = Сз. 
Таким образом, дифференциальное уравнение, описывающее накопление жидкости в реакторе: 
оао (3) 
% 5 
Для определения динамики изменения концентрации раствора в реакторе необходимо найти дифференциал 
количества вещества: 


УЕУ.С, (4) 
где У — это объем вещества, л; С — молярная концентрация, моль/л. 
Из (4) и (1) получаем: 
аа ь: | (5) 
& Ш ф @ 


На изменение концентрации раствора влияют потоки щелочи и кислоты с определенной концентрацией. 
Значит: 


ай ас 
5-м. | - (@-Сыи-,-Сл). (6) 
Из (6) выведем изменение концентрации и упростим выражение: 
ас 1 
=— (С, (Сион ие (п 
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Представим систему дифференциальных уравнений, описывающих динамику изменения объема и 
концентрации раствора в биореакторе: 


Ца +6,-С,), 
@ 5 
(8) 
ас _ 1 
и? - яр (С (мон -С+ „(С — Сна )). 


Зная основы химии и принимая во внимание (8), можно определить концентрацию катионов Н ([Н*]). 
Концентрации ионов водорода ([Н*]) и гидроксид-ионов ([ОН ]) в дистиллированной воде при 25 °С равны и 
составляют по 10” молыл. Это следует из определения ионного произведения воды [11]: 


к = [Н*|х [ОН ] =Ю “тор /Р (а 25 °С). (9) 
Также стоит учитывать, что: 
Т=[ОН`]-[Н"] = [Ма] СГ]. (10) 
Теперь найдем [Н"] (обозначим, как Си) через квадратный корень из предыдущих выражений с учетом 
концентрации раствора по уравнению (4) [14]: 


О и 
С, =УК,Т=0, (11) 
с, =-5 т 


Здесь ионы 1= 4. 

Для создания системы регулирования с обратной связью по кислотности в виде водородного 
показателя (рН) используем отрицательный десятичный логарифм из концентрации, рассчитанной по 
уравнениям (7), и получим кислотность раствора при нормальных условиях, то есть при 25 °С. 

Однако содержимое биореактора не всегда находится в нормальных условиях. Например, для создания 
искусственной среды, свойственной ш У1уо среде желудочно-кишечного тракта домашней птицы, в реакторе 
поддерживается температура 42 °С. Поэтому на высокоточных измерителях рН получили матрицу изменения 
рН буферных растворов в зависимости от температуры. Показания фиксировали по измерительной системе с 
термокомпенсацией для применяемого электрода [13]. С этой целью использовали данные эксперимента и 
сведения из открытых источников°. 

В итоге построили статическую модель. Ее получили полиномиальной аппроксимацией методом 
наименьших квадратов. Зависимость рН от температуры визуализирована на рис. 2 [14—16]. 


Реальный рН > с 


Рис. 2. Поверхность, описывающая зависимость рН от температуры (12) 


3 Ашотанс Тетрегацие Сотрепзаноп ш РН Меазигетеи/ Нопфа// Вопфа.сот: [сайт]. ВТ: Вирз:/Амум.Вофа.сот/езр/мацег- 


диаНеу/зирронАесьщса]-Нрз/БепсН-пеегз/ащюотайс-{етрегаиге-сотрепзаноп-ш-рЬ- 
шеазигетеп И? ии _зоигсе=ив\у&ини_тедат=301&хани_сатраеи=иб\-ге@ тес! (дата обращения: 25.12.2022). 
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Е(®) = р00+ р10-х+ р01- у+ р20-х? + р. х. у+ р02- у? + 2) 
+р30-х* + р21-х’-у+р12-х- у”. 


Коэффициенты уравнения представлены в таблице 1. 


Таблица 1 
Коэффициенты уравнения (11) 
Коэффициент Значение 
200 0,154 
р10 —0,00511 
ро1 0,899 
р20 —2,998 10° 
РИ 0,00333 
р02 0,0130 
р30 -3,226. 10” 
р21 1,588 . 10° 
р12 —4,8. 10“ 


Для определения переходных процессов в системе построим математическую модель в прикладном пакете 
МаЧаь ЗитиНак (рис. 3): 


ОН 
ОН 
Реальный РН У РН(25С) 
Н 
Модель системы Модель Н 
управления дозацией изменения 


объема 


РН Реальный рН 


Статическая модель 
зависимости рН от 
температуры 


Модель 
химических реакций 


Рис. 3. Общая структура математической модели в МаНаб ЗипаПаК 


Рассмотрим созданную систему по блокам. Первый — блок двухканального регулятора по кислотности с 
обратной связью. В цепь обратной связи введены шумы для моделирования случайных событий реальной 


системы измерения кислотности и погрешности измерительных приборов. 


ила сы 
Требуемый рН ПП 
нь 


Контроллер Н 


+ ий Ошибка считывания рН 
Реальный рН 
ь Со 


Рис. 4. Структура регулятора рН, блок «Модель системы управления дозацией» 


Второй блок — это система контроля объема раствора в реакторе. Функция Уоате свапее то4е 
моделирует включение откачивания. Так обеспечивается постепенный сброс раствора до необходимого объема 


при его приближении к максимуму. 


Д.Ю. Донской. и др. Математическая модель системы управления рН вт ъйто модели желудочно-кишечного тракта 


Начальный 
объем 


Контроллер сброса 


Уо 


у ощ 


х 1/5 (реактор) 


Объем в высоту 
столба жидкости 


Графики 


и 


У 


Рис. 5. Структура модели регулирования объема содержимого в реакторе, блок «Модель изменения объема» 


Третий блок — основной. Он реализует расчет химической реакции кислоты и щелочи при добавлении в 
реактор и вычисление рН конечного раствора при нормальных условиях. Функция МеибаПтайоп определяет 
концентрацию ионов водорода в растворе [12]. Блок Сопсепиайоп ш рН преобразует концентрацию в 
водородный показатель рН. 


-(Расчет реакции) 


Концентрация рН 
при н.у. 


Нео 


РН (25 С) 


Рис. 6. Структура модели, описывающей химические реакции, блок «Модель химических реакций» 


Функция статистической модели в четвертом блоке (Тетрегаеиге соггесНоп) позволяет учесть температурные 


изменения раствора. 


Температура 
раствора 


Температурная коррекция 
водородного показателя 


Реальный РН 
раствора 


Рис. 7. Структура модели, описывающей изменение кислотности от температуры, блок 


«Статическая модель зависимости РН от температуры» 
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Для определения адекватности переходных характеристик результирующей математической модели 
проводились эксперименты на имитационной модели. Реальная модель представлена системой 
автоматизированного управления искусственным ЖКТ домашней птицы с одним блоком управления и 
биореактором. Оборудование представлено на рис. 8. 


а) 6) 


Рис. 8. Оборудование для проведения исследований: а — система управления параметрами среды в биореакторе; 


б — элементы биореактора. 1— управляемые перемешивающие устройства; 2— нагреватели; 
3— порты для датчиков и системы дозации жидкостей 


Опишем методику проведения эксперимента. В математической модели учитываются характеристики 
реальной системы: 

— заданная кислотность; 

— параметры реактора (площадь основания и высота); 

— концентрации доливаемых веществ (кислоты и щелочи); 

— начальная концентрация раствора; 

— начальный объем раствора; 

— температура раствора на протяжении всего эксперимента; 

— дискретность процесса регулирования [14]. 

Математическая модель представляет собой идеальную систему: смешивание растворов в реакторе 
происходит мгновенно и равномерно. В реальных условиях это невозможно, поэтому для равномерного 
растворения добавленной кислоты или щелочи необходимо принудительное перемешивание. В имитационной 
модели данная проблема решается магнитной мешалкой с управлением $САПА“. 

К сожалению, такой тип перемешивающего устройства наводит на рН-электрод электромагнитные помехи. 
Для исключения помех устройство перемешивания кратковременно отключалось. Дискретность измерений и 
работы регулятора составила 20 секунд, что обусловлено экспериментально установленным временем 
нормализации показаний рН электрода в реакторе (по документации электрода возможна дискретизация до 1- 
2 минут) ”°. В целях обеспечения безопасности регулятор реальной системы ограничен в максимальном объеме 
дозируемой кислоты или щелочи, что учитывалось в модели. 

На рис. 9 сравниваются переходный процесс математической модели и реальной системы управления. РН- 
электрод — это чувствительный элемент, на который воздействуют внешние электромагнитные помехи. 


4 От англ. Зирегу1зогу Сопто! Апа Рава Асаш$#юоп — диспетчерское управление и сбор данных. Программный пакет для создания или 
поддержки работы систем сбора, обработки, отображения и хранения информации об объекте мониторинга или управления. 

° РЬ-4502< РН шеег саНьаноп поез / Шопе01,Ы0е / И опе01510°. Му @есиотсз Фу ргодесз: [сайт]. ОВ: 
Брз://Фоп201 .6102/2019/04/26/рь-4502с-рЬ-теег-саШгаНоп-побе5/ (дата обращения: 30.11.2022). 

° Н7Отм рН Сисай ПаазВее / АЧаз Заепийс МС // НШез.аМаз-зслепайс.сот: [сайт]. ОВГ:  Бирз://ез.аЧаз- 
зс1епайс.сош/рН_ЕЙО_РаазВееЕ.рАЕ (дата обращения: 30.11.2022). 
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Рис. 9. Переходные характеристики математической модели и реальной системы при регулировании от 6,2 до 2,15 рН 


Один из экспериментов проводился при заданной кислотности раствора 2,15+0,05 рН. Кислотность 
исходного раствора — 6,2 рН. Температура раствора в реакторе — 22,25 °С. Для изменения кислотности 


использовались растворы 


0,1 моль/литр МаОН и 0,1 моль/литр НС. 


дискретном пропорциональном регулировании. 


Алгоритм дозации заключался в 


На рис. 10 представлены интегралы по потокам из дифференциального уравнения (7). 


6:10-8 


4-10 


Интегралы по потокам 
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50 


100 
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200 
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| — юн 

--° @>[Н] 

- С! [поток р-ра] 
---_ О›[поток р-ра] 


300 


350 


400 


450 — Время, с 


Рис. 10. Влияние каждого потока на реакцию нейтрализации в ходе экспериментов 


Обсуждение и заключения. Анализ результатов моделирования показал, что потоки с концентрацией 


раствора в реакторе не вносят весомого вклада в регулирование всей системы, поэтому ими можно пренебречь. 
Тогда система дифференциальных уравнений цифровой модели (8) приобретет следующий вид: 


а<_ 1 


а 1 
в = < (6, +С(, —6. }, 
(С, } Стон —С, "Сна ). 


@& 51 


(13) 


По результатам экспериментов с учетом концентрации раствора в реакторе относительная погрешность 


установившейся кислотности раствора для математической модели составила 0,1 %, а для реальной системы с 
тем же алгоритмом — 0,35 %. 


Итак, в математической модели учтем погрешность показаний электрода и пренебрежем перемешиванием 


растворов. Абсолютная погрешность результатов регулирования находится в Допустимых пределах и не 


превышает 0,1 РН от заданной величины. Это означает, что разработанная математическая модель адекватно 


описывает кривую реакции нейтрализации и в будущем позволит подобрать оптимальные алгоритмы 
регулирования, что было бы трудоемким процессом в реальной системе. 
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Разработанная модель решает поставленную задачу. Она представляет собой математический инструмент 
для создания методик и алгоритмов регулирования РН с учетом температуры содержимого мини-реактора, 
начальной концентрации основного раствора и концентраций доливаемых растворов. Такой подход позволяет 
совершенствовать существующие алгоритмы управления. 
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